Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,393 C→T 11.3% pseudogene (1002/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,3930CT11.3% 100.7 / 4.5 51pseudogene (1002/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base C (24/20);  new base T (3/3);  total (27/24)
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.52e-01

GCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACG  >  NC_000913/1466155‑1466630
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
gggATCTGGATGGCAACGGCGGCGGTCAAGGTATTGATATCACCGGTGGCAGGGCTACGGTGGATAACCAGGGCACCATGACCGTCCCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCTAGAGCACCATCAGCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCGCGGCAAACAACAACGGCAAAACCTCCGTTGACGGc                                                                                                                                                                                                                                                <  2:724311/237‑1 (MQ=255)
            tCAGGGGGGGCGTCCCCGGTTTTGATATCCCCGGGTACAGCGGGACCGGGGATAACAAGGGCCCCAGGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCAc                                                                                                                                                                                                                                     <  1:675557/238‑1 (MQ=17)
              aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGATTAATGGTGATGACGCCACGGCAAAAAACAACGGCAAAACCACCGTTGACCGCAAGGATTCCACCgg                                                                                                                                                                                                                                  >  2:842494/1‑239 (MQ=255)
                 ggcggcggTCAAGGTATTGAAATGTCCGGTGAACGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGATAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa                                                                                                                                                                                                                                <  2:854765/238‑1 (MQ=255)
                                  tGATATCACCGGTGACAGCGCGACGGTGGATAACCAGGGCACCATGACCGTCACCGACCCGCAGTCCATGGGGATCCAGATCGACGGCGACACGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGTGATGCAGCCACGGCAAACAACACCGGCAAAACCACCATTGAAGGCACGGATTCCACCGCTACGGCAATAAACggacaac                                                                                                                                                                                                              >  1:259169/1‑234 (MQ=14)
                                      atCACCCGTGACCGGGCGACGGGGGATAACCAGGGGACCATGCCCGTCCCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCCCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACgg                                                                                                                                                                                                          <  2:369237/239‑1 (MQ=18)
                                           cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCCAACACCAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGAACGGAAATCACCGGAAATAACGGCAAAg                                                                                                                                                                                                     >  2:858654/1‑239 (MQ=255)
                                           cGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGgaagg                                                                                                                                                                                                     <  2:862592/239‑3 (MQ=255)
                                              tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCAGGAGTCCATGGGTATCCAGATCGACGGTGACACGGCCATCGTCACTAACGCAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGTAGGACGCCACGGCAAACAACAACGCCAAAACAACCGTTGACGGCAAGGATTCCACCGGCACGGAACTCAACGGTAATAACaaaaaaagaa                                                                                                                                                                                                  >  1:735053/1‑235 (MQ=17)
                                                    cgcgACGGTGGATAACAAGTGCACCATGCCCGTCACCGATCCGGGGTCCATAGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAAGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGGGATGACGCCACGGCAAAAAACAACGGAAAAACCCCCGGTGACGGCAGGGCTTCACCCGGCACGGAAATTAACGGAAATAACCGGAAAggaaataaga                                                                                                                                                                                            >  1:619982/1‑230 (MQ=255)
                                                     ggggtttgtgaaaaCAAGGGCACCATGACCGTCACCGATCCTGCGGTCCTGGGTATCCAGAGCGACGGTAGCAAGGCCAAGGCAAATAACGGAGGCGAGAGCACCAAACACAACCGTGGCAGCGGCCTCCAGATTAAGGGTGATGATGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGa                                                                                                                                                                                          <  1:704708/227‑1 (MQ=17)
                                                         cgggtgATAAAAAGGGTACCATGACCTTACCCGTTCCGGAGTCCATGGGTTTCCAGATCGACGGGTAAAAGGCCATCGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                                                                                                                                                                                       <  1:426388/233‑1 (MQ=18)
                                                         cGGTGGATAACAAGGGCACCATGACCGTCACCGGTCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACgg                                                                                                                                                                                       <  2:97730/238‑1 (MQ=255)
                                                          ggTTGATAACCAGGGCCCCCTCCCCGTCACCGACCCGGAGTTCATGGGTCTCCAGATCGGAGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGTTGATGGCGGCGCGGGAAGTAACAACGGGAAAGCCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCg                                                                                                                                                                                     <  1:282641/239‑1 (MQ=25)
                                                                        gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATAAACGGAAATAACGGAAAAGGGATTCAGGACGCGGACCTgggttcacgc                                                                                                                                                                       >  2:843491/1‑231 (MQ=17)
                                                                            cATGACCGGCCCCTATCCGGGGGCCTTGGGTACCCAGATCGGCGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGGTATTCAGGACGGCGATCTGGATGTCAgcggcg                                                                                                                                                                   <  1:842494/239‑1 (MQ=18)
                                                                              tGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcggc                                                                                                                                                                 >  2:494290/1‑239 (MQ=255)
                                                                                  cgtgaaggATGCGGTGTCAATGGGGAACAAGAGCGACGGGGAAAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGGCCCCTGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                                                                                                                                                              <  2:842806/231‑1 (MQ=18)
                                                                                      aCCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCAACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCCATCTGGATGTCAGCCGCGGCGGCCACg                                                                                                                                                          >  2:364180/1‑238 (MQ=255)
                                                                                               gAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAACAGTTATTCAGCAAGGCGTCCCGAATGTCagcagcccgggtacgggtattaaaa                                                                                                                                                 >  2:724751/1‑219 (MQ=18)
                                                                                               gAGTCCATGGGTATCCAGATCGACGGGGACAAGGCCATCGTCAATAACGCAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGACACGGCAAACACCACCGGCAAACCCACCCTTGACGGCAAGGATTCCACCGGTCCGAAAATAAAATGCAATAACGGAAAGGTTATTTGAGCACGCGCCGTGGATGTCAGCGGCGGGCGGCACGGGAtgagca                                                                                                                                                 >  1:587283/1‑233 (MQ=255)
                                                                                                    cATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGCAATCAACGGTAATAACGCAAAAGTTATTACGGACGGCGATCTGGATGTCACCGGCGCCGGTCCCGGCATTGagattccc                                                                                                                                            >  2:468841/1‑234 (MQ=18)
                                                                                                     aTGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGCTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGGATTGATTTAACCg                                                                                                                                           >  2:807278/1‑238 (MQ=255)
                                                                                                         gTATCCAGATCGGTGGTGAAAAGGCCATCGGCAAAAACGAAGGCGAGAACCCGCACACCAACGGTGACCGCGGCACCCAGATTAAGGGTGATGACGCCACGGTAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGa                                                                                                                                       <  2:786074/238‑1 (MQ=17)
                                                                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCCGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGCCGGCGCCTTGGATGTAACCGACGGGGGGCACGGGATTGAAATCAACGGCGACaag                                                                                                                                   >  1:854081/1‑237 (MQ=18)
                                                                                                            tCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTAACGGGATTGATATAACCGGccaacg                                                                                                                                    >  2:856889/1‑233 (MQ=255)
                                                                                                                  tccgtggTGAAAAGGGCATGCTCAAAAACGGAGGCGTGAGGTGCAACACCAACCGTGACCGCGGCCTTCAGATTACTGGTGATAACGCCACGGTAAACAACAACGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCaa                                                                                                                               <  2:422871/232‑1 (MQ=14)
                                                                                                                  tCGACGGTGACAAGGCCAACGACAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCaa                                                                                                                               <  2:600731/237‑1 (MQ=255)
                                                                                                                   cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGCCACGGTATTGGTATACCCGGCGACAGCGCAACgg                                                                                                                            >  2:296980/1‑239 (MQ=255)
                                                                                                                              aGGACATCGTCAATAAAGAAGGCGAGAGCACCATCACCAACGGTTGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                 <  2:735053/239‑1 (MQ=255)
                                                                                                                                 ccgtCGTCAATAACGAAGGCGAGAGGACCAACAACAACCGTTGCACCGGCCCCCGGATTAATGGGGATGAAGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAggg                                                                                                               <  2:825209/235‑1 (MQ=17)
                                                                                                                                   aTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATACCCGGCGCACGCCCAACGGTGGGAAAAAAgggaa                                                                                                             >  2:261446/1‑236 (MQ=17)
                                                                                                                                     cGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGGAATAACGGAAAAGTTATTCAGGACGGGGATCTAGATGTCAGCGGCGGCGGGCACGGGATTGGTATAAACGGCGCAAGCGCAAAGGGGGGAAAAAAggggagcc                                                                                                          >  1:716532/1‑234 (MQ=14)
                                                                                                                                              cgaagggccgagCACCAACACCAACGGTGGACCCGGCCCCCAGATTAATGGTGATGACGCCACTGAAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGGCCGTCa                                                                                                 <  2:755040/230‑1 (MQ=14)
                                                                                                                                                 aGGCGAGGGCACCTTCCCCAACGGTCGAACCGGCACCCAGATTAAAGGAGAAGGCGCCACCGCAAACAACAACGGCAAAACCATCGTTGACGGGAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCAGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCAcc                                                                                               <  1:364180/238‑1 (MQ=14)
                                                                                                                                                  ggCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAGCGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCACCGGCAATAACGGGAAGGTTATTCAGGACGGCGATCTGGAGGCCAGCGGCGGCGGGCCAGGTATTGAAATAACCGGCGCCAGCGACACGGGGGGTAACAAGGGCACCATAAACGACAcca                                                                                              >  1:346800/1‑237 (MQ=11)
                                                                                                                                                        agCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGCGACAGCGCAAAGGGGGGTAACAAGGGAACCATGACCGTCACCGATCCGGa                                                                                       >  2:830415/1‑239 (MQ=255)
                                                                                                                                                             cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGGGGATCACAAGGGCACCATGACCGTCACCGATCCCGAGccc                                                                                   >  2:417879/1‑238 (MQ=255)
                                                                                                                                                                     aCGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATc                                                                          <  2:854081/239‑1 (MQ=255)
                                                                                                                                                                             accggcaccCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGGACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGAACGGGATTGAGATCACCGCCCACACGGCAACGGGGGAGAACAAGGGCACCAGGACCGACCCCGATCCGCGGGCCATCggcaccccacttaa                                                                  >  1:715201/1‑227 (MQ=14)
                                                                                                                                                                                        agttAATGGTGAAGAGGACACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAgg                                                       <  2:776380/237‑1 (MQ=255)
                                                                                                                                                                                         aTTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCgcacaagc                                                      >  2:813074/1‑232 (MQ=255)
                                                                                                                                                                                           tagtggggttgaCGCCACGGTAAACAACAATGGCAAAACCACCGTGGACGGTAAGGATTCCAGCGGTACGGAAAACAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATTGCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCgg                                                    <  1:193843/230‑1 (MQ=14)
                                                                                                                                                                                                   atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTTACGGGGCCCcggcgggtttttaa                                            >  2:161692/1‑231 (MQ=255)
                                                                                                                                                                                                                aaaaaaCAACCCAAAAACCCCCGTGGCCCGCAAGGATTGCACCGGTTCGGAAATCAAACGTAAAAACGCAAAATATATTCAGGACGGCGAGCTGGATTGCGTCGGAGGCGGTCACGGTATGTATAGCACCGGCTGCAGCGCAACGGTGGATAACAAGGGCGCCATGACCGTAACGGATCCGGAGGCCATCGGTAACCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgag                                <  1:468841/238‑1 (MQ=255)
                                                                                                                                                                                                                 aacaacaacGGCAAACCCTCCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCgaga                               <  2:183817/238‑1 (MQ=255)
                                                                                                                                                                                                                     acaccGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAAGAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGcc                          <  2:587283/239‑1 (MQ=255)
                                                                                                                                                                                                                         tggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCCTGACCGTCCCCGGACCGGGGGCCCTGGGGAGCCAGGACGCCGGGGACCAAGCCGATGTTAAAAAAGaaagggaaaaccacaaac                      >  1:41117/2‑223 (MQ=11)
                                                                                                                                                                                                                         tggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAAGGTATTGATATCACCGGTGACAGCGCGACGGTGGGTAACAAGGGGACCATGAACGTAACCGTTCCGGATTCCATTGGTGTCCCAGTAGACGGAGAACAGGCCGTTGTTaacaaacaaacagcaaacgcccta                       >  1:183817/2‑219 (MQ=11)
                                                                                                                                                                                                                         cGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGCCGGCGGTCACGGAATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAAGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAAAAACGAAGGCGAGAACGCCATc                       >  2:259291/1‑238 (MQ=255)
                                                                                                                                                                                                                           gCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGGCCATGGGTATCCAGATTGCCGGGGACCAGACCGTCGGTaacaaaaaaagaggaagggacaaccca                    >  1:394406/1‑217 (MQ=11)
                                                                                                                                                                                                                                 cccacggtGACGGCAAAGATTCCCACGGTAAGGAAACCAACGGGAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGt              <  2:394406/232‑1 (MQ=17)
                                                                                                                                                                                                                                              cAAGGGTTCCCCCGGTACGGAAAACAACGGTAATAACGGAAAAGTTATTCAGGACGGGCATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg  <  2:411935/238‑1 (MQ=255)
                                                                                                                                                                                                                                              cAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACg  <  2:716532/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
GCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACG  >  NC_000913/1466155‑1466630

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: