Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,529,185 C→T pseudogene (1284/2037 nt) rhsE → pseudogene, Rhs family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,529,1850CT97.4% 82.3 / ‑6.0 38pseudogene (1284/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (15/22);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATG  >  NC_000913/1529014‑1529482
                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         
ggtgTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTa                                                                                                                                                                           >  1:598226/1‑300 (MQ=14)
      tcgggCTGAACTGGTGGGGCTGTTGGACGGGCTGAAGGAAGAAATCCGGGCAGCCCGCGTGGGCAGTGAGAGCCGGGCGTGGCTTGCGCAGTGCGGCCTGACGGTTGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGGGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTg                                                                                                                                                                    <  1:214896/298‑1 (MQ=11)
       ccGGCTGAACAGGTGCGGTTGCTGCACAGTCTGGAGGAAGAAATCCGCGCAGACCGCGTGAGCAGTGAAAGCCGGGAGTGGCTTGCGCAGCGCGGGCTGACGGTCGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTCACCACTGCGACCATCTGGGACTGCTGCCTGCGCTTTTCAGCGAATACGCTAATACGGCTTGGTGCGCGGAGTATGATGAAAGGCGCAACCAGCCTCATTAGGCGAACACGCATCACCTCCACCAGCCGTACAGTCTg                                                                                                                                                                    >  2:213183/1‑300 (MQ=11)
         ggCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGTAAAGCTCATCTTTATCAATGCGACCACCGGGGACGGCCGACGGCGCTTGTCCGCGGAGACCGAAATAAGGCGCGGCGCGGGGAGTATGATGAATGGGGCGACCATCATAAAGAGGAGGACTGGCGTCACATCCACCAACCGTTACGTccgagg                                                                                                                                                                 >  2:214894/1‑296 (MQ=11)
         ggCTGAACGGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAGAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGcc                                                                                                                                                                  <  2:346850/300‑1 (MQ=9)
             gAAGTGGTGCGGCGGCGGAGCAGGCTAGAGGAAGAAAGCCGGGCAGCCGGCGTGGGCAGTGAAAGGCGGGCGTGGGTTACGGAGGACGGACGGACGATGGAGCAAATGGACAGAGAGGTGGAACCGGGATACAAGCCGGGGGGAAAAGCTCATCTTTATCACGGCGACCACCGGGGACTGCGGCTGGCGCTGATCAGCGGAGACGGCAATACGGCGTGGCGCGGGGAGTAGGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAggg                                                                                                                                                              <  2:385722/300‑1 (MQ=11)
              aaCTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgc                                                                                                                                                             <  1:507847/300‑1 (MQ=14)
                              ggACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGTCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATTATTGCGACCACTGGGGATTGCCGCCGGCGATTATCAGAGAAGACGGCAATACGGCGGGGCGCGGGGAGGACGATGAAGGGGGCCACCAGCTTAATGAGGAGAACCAGCAACACCTGAACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                                                            >  2:212333/1‑301 (MQ=11)
                               gcaaGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAACGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                                                            <  2:323183/297‑1 (MQ=12)
                                    gCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                      >  1:432742/1‑301 (MQ=18)
                                    gCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCggg                                                                                                                                       <  1:212335/300‑1 (MQ=18)
                                     cTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTTAAAGCCGCGCGTGGCTTGCGCAGTGCGGGCTGACGGTGAAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGTCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                     <  1:570099/301‑1 (MQ=255)
                                      tGGGGGAAGAAATCCGGCCAGACCGCGTGAGCAGTGAACGCCGTGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGCCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                     <  2:598220/300‑1 (MQ=255)
                                      gggaggaAGAAATCCGGGCAGACAGAGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACGGGCCAGAAAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                     <  2:43207/299‑1 (MQ=11)
                                       ggaggaAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGTAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                      >  2:128894/1‑298 (MQ=16)
                                               aaaTCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACtata                                                                                                                            <  1:110088/300‑1 (MQ=18)
                                                        cAGACCGCGTTGGCAGTGAAAGCCGGGAGTGGCTTGGGCAGTGCGGGCTGACGGTGGAGAAACTGGCCAGACAGGTGGAGCCGGAATACACGCCGGCGCGAAAGGCTCAGCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATAAGGAGTCGGGGCTGTACTATAACCGTCAcc                                                                                                                   <  2:109119/300‑1 (MQ=11)
                                                           ctcttccgatctgTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCa                                                                                                               <  1:512650/289‑1 (MQ=255)
                                                                gTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTCCCTTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTTACCGGCACTAc                                                                                                           >  1:115402/1‑300 (MQ=255)
                                                                gTGAGCAGGGAAAGCCGGGCGGGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTAc                                                                                                           <  2:33673/300‑1 (MQ=255)
                                                                  gAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGGACTATAACCGTCACCGGCACTACGaa                                                                                                        >  1:194590/1‑300 (MQ=255)
                                                                  ctctttcccaaaacgacgcggttccgagctgTGCGGGAGGAGGGTGGAGCAACTGGCCAGGCAGGGGGGGCCGGAATACACACCGGCGCGAAAAGCTCATATTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAt                                                                                                        <  2:570483/271‑1 (MQ=11)
                                                                      aGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCGGACAGGTGGAGCCGGAATACACACCGGCTCGAAAAGCTCTTCTTTATCCCTGAGAACACCGGGGACTGCTGCTGGCGCTTACCAGCGAAGAGGGCAATACGGCGCGGCGCGGGGAGTACGATGTATGGGGCAAACAGCTTAATGGGGAGAATCCGCACATCCTGCACCAGTCGCACAGTCCTCCAGGGCAGCAGCATGACGAGGAGTAGGGGACTTATTATAACTGCCACCGGGAcatcgagacgg                                                                                                    >  2:96939/1‑291 (MQ=2)
                                                                       gTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCCGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCAATatcggaagag                                                                                                   >  2:512645/1‑292 (MQ=255)
                                                                         gAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACCGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCCCGAAAAGCCCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATTAGCGAAGACGGTATTACGGCCTGGCGCGGGGAGTATGACGAACGAGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGCACGGTCTGCCAGGCCAGCCTCTTGATGAGGAGTCGGGGGCGTACTGTCACCGTCACCGGTACTCCGATGCGCTg                                                                                                  >  2:21789/1‑300 (MQ=11)
                                                                                   tacccgacgctcttccgatcttGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGta                                                                                       <  2:224654/280‑1 (MQ=11)
                                                                                    gTGGCGTGCGCAGTGCGGGCGGACGGTTGAGCAACTGGCCAGACAGGTGGAGCCGGAAGACACACCGGCGCGAAAAGCTCATCTTTATCGCTGCGACCACCGGGGACTGCCGATGTCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGTAGTATGATGAATGGGTCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCACCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGta                                                                                       <  1:636271/300‑1 (MQ=11)
                                                                                                gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAtagatcggaagagcatactgctgaactccat                                                                          >  1:570488/1‑271 (MQ=11)
                                                                                                   cGGGGTGACGGTAGAGCAACGGGCCAGAAAGGGTGAGCAGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCTCCGGGGACTGCCGCGGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATAATGAATGGGGCAACCAGCTTAATGTGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGAccc                                                                        <  2:115400/300‑1 (MQ=14)
                                                                                                        tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGtaagatcggaagagcacacgt                                                                    >  1:224656/1‑280 (MQ=11)
                                                                                                                 aGCAACAGACCAGTCATGTGTGGCCGGAAGACACACCGGCGCGAAAAGCTCATCTTTATCACGGCGACCACCGGGGACGGCCCCTGGCGCTTAACAGCTAAGACGGCAGTACGGCGTTGCGCGGGGAGTATGATGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCAGCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGGCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTGGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGagag                                                          <  2:307978/300‑1 (MQ=14)
                                                                                                                   cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAgg                                                         <  1:128896/299‑1 (MQ=255)
                                                                                                                        ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCGCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatg                                                   <  1:271280/300‑1 (MQ=255)
                                                                                                                           aaGAAAGGTGGGGCCGGAATCCACACCGGCGCGAAAAGCTCATCTGTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAat                                               <  2:194588/300‑1 (MQ=255)
                                                                                                                                                     cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCAt                     >  1:228678/1‑301 (MQ=255)
                                                                                                                                                                gctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAATTGTCACCGGCACTACGACCCGTTGCAGGGGCGGTATATCACCTCGGACCCGATTGGGTTGCGGGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTgttcg          >  2:53801/3‑297 (MQ=17)
                                                                                                                                                                   cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGAc        >  1:23285/1‑300 (MQ=255)
                                                                                                                                                                        ttATCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGAACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGGATATGTATCAGTATCCGTTGAAACCCATACAAGTGATAGACCCAATg  <  1:21789/301‑1 (MQ=255)
                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         
GGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATG  >  NC_000913/1529014‑1529482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: