Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,529,180 T→C pseudogene (1279/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,1800TC97.3% 81.8 / NA 37pseudogene (1279/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (16/20);  minor base A (0/1);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGAC  >  NC_000913/1529014‑1529476
                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                        
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)
                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                        
GGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGAC  >  NC_000913/1529014‑1529476

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