Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,595,660 A→G pseudogene (2328/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,595,6600AG98.6% 278.3 / NA 73pseudogene (2328/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (45/27);  total (46/27)
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

AGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGAGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTAATCTCCTGC  >  NC_000913/1595513‑1595806
                                                                                                                                                   |                                                                                                                                                  
agcagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGCTTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGt                                                                                                                                                  >  1:927379/1‑150 (MQ=255)
agcagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTg                                                                                                                                                 >  2:1217426/1‑151 (MQ=255)
agcagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTg                                                                                                                                                 <  1:1209826/151‑1 (MQ=255)
agcagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTg                                                                                                                                                 >  2:378571/1‑151 (MQ=255)
agcagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTg                                                                                                                                                 <  1:48868/151‑1 (MQ=255)
 gcagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTg                                                                                                                                                 >  1:91179/1‑150 (MQ=255)
  cagGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGaa                                                                                                                                               <  1:244329/151‑1 (MQ=255)
     gAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGt                                                                                                                                            <  1:926538/151‑1 (MQ=255)
     gAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGt                                                                                                                                            <  2:165045/151‑1 (MQ=255)
        tatCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACt                                                                                                                                         >  2:1057101/1‑151 (MQ=255)
        tatCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACt                                                                                                                                         >  2:1201916/1‑151 (MQ=255)
            cggaTCGAAAATGATGGTACCTCCATTCATGGCCAGTTCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGAATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCt                                                                                                                                     >  2:869567/4‑151 (MQ=255)
            cGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCt                                                                                                                                     >  2:697526/1‑151 (MQ=255)
            cGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCt                                                                                                                                     >  1:592691/1‑151 (MQ=255)
                tCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTc                                                                                                                                  >  2:63983/1‑150 (MQ=255)
                tCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTc                                                                                                                                  >  1:673886/1‑150 (MQ=255)
                 cGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCa                                                                                                                                 <  2:1018690/150‑1 (MQ=255)
                 cGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAg                                                                                                                                <  2:184266/151‑1 (MQ=255)
                  gAAAATGATGGTACCTCCATTCATGGCCAGTCCGTCAATGTGTTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGTCTGCAACATCGCGTGAGTAAGCGCTGCGGTGTTGTTTCGTTCCAGGGTGAAGGTACTGTCTTTc                                                                                                                                  <  2:1114921/148‑1 (MQ=255)
                  gAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGt                                                                                                                               >  1:1196564/1‑151 (MQ=255)
                   aaaaTGATGGTACCTCCATTCATGTCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGTGACTGTCTTTCAGtt                                                                                                                              >  2:914115/1‑151 (MQ=255)
                      atgatgGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTggg                                                                                                                           <  2:805731/151‑1 (MQ=255)
                         atgGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCa                                                                                                                         >  2:526451/1‑150 (MQ=255)
                          tgGTACCTCCATTCATGGCCAGACCACCAATGGATTGCTCTCCCACTTTTACCGATGTGGTATTTTCACTGCCAGACTGCAACCTCGCGTGTGTAAGCGCAGCGGTGATGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTCGGCAac                                                                                                                       >  2:679212/1‑151 (MQ=255)
                             tACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACAcc                                                                                                                    >  2:553220/1‑151 (MQ=255)
                              aCCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACAcc                                                                                                                    <  2:31047/150‑1 (MQ=255)
                              aCCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACGTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACAcc                                                                                                                    <  2:1044510/150‑1 (MQ=255)
                                   cATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGa                                                                                                               >  2:651601/1‑150 (MQ=255)
                                    aTTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAAt                                                                                                             <  1:219889/151‑1 (MQ=255)
                                      tCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAAtt                                                                                                            <  1:1159039/150‑1 (MQ=255)
                                         tGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGt                                                                                                        >  1:875497/1‑151 (MQ=255)
                                             cAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCt                                                                                                    >  1:661379/1‑151 (MQ=255)
                                             cAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCt                                                                                                    >  2:1183722/1‑151 (MQ=255)
                                                 ccaccaATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGtt                                                                                                 <  1:91910/150‑1 (MQ=255)
                                                 ccaccaATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGtt                                                                                                 <  1:1065554/150‑1 (MQ=255)
                                                     caATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATg                                                                                            >  2:680844/1‑151 (MQ=255)
                                                         ggATTGCTCTCGAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTCCTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGt                                                                                         >  2:832191/1‑150 (MQ=255)
                                                           aTTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTaaa                                                                                      <  2:91182/151‑1 (MQ=255)
                                                           aTTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTaa                                                                                       >  1:721045/1‑150 (MQ=255)
                                                                  tcCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCCGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCaaa                                                                                >  1:1237663/1‑150 (MQ=255)
                                                                     aaCTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATc                                                                            >  1:121411/1‑151 (MQ=255)
                                                                          ttACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTc                                                                       >  1:555641/1‑151 (MQ=255)
                                                                           tACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCg                                                                      >  1:619858/1‑151 (MQ=255)
                                                                              ctctgtgGTATTTTCACTGTCTTAGTGCAACATCGCGTGGGTAAGCTCAGCGGTGTTCTCGCGTTGCAGGGTGAAGGGACTGTCTTTCAGTTGGGCAACACCGGCGACTTCAGTCCCTATTGCATGGGGGAACCCAAACATCTTGTc                                                                       >  2:1127501/4‑147 (MQ=255)
                                                                                aTGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCggagg                                                                  >  1:1240496/1‑150 (MQ=255)
                                                                                      tATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGAc                                                             >  2:294862/1‑149 (MQ=255)
                                                                                       aTTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGc                                                          >  1:565998/1‑151 (MQ=255)
                                                                                         tttCACTGTCAGAGTGCAACATCGCGAGAGGAAGCTCAGCGGTGTTGTCGCGGTCCAGATTGAAGGTACTGTCTTTCAGTTTGGCAACACCGCCGAATTCAGTCCCTGTTTCAGTGTTAGCGGCAAACATCTTGTCGGAGGATGAAAGc                                                          >  2:620185/1‑149 (MQ=255)
                                                                                               tGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCg                                                  >  2:978900/1‑151 (MQ=255)
                                                                                               tGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCg                                                  >  2:857509/1‑151 (MQ=255)
                                                                                                  cAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCa                                                <  2:668166/150‑1 (MQ=255)
                                                                                                     aCTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCa                                             >  1:610623/1‑150 (MQ=255)
                                                                                                         cAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGcc                                         >  1:595574/1‑150 (MQ=255)
                                                                                                          aaCATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCAt                                       >  1:952459/1‑151 (MQ=255)
                                                                                                             aTCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGc                                    >  1:83321/1‑151 (MQ=255)
                                                                                                               cgcgTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCt                                  >  1:299469/1‑151 (MQ=255)
                                                                                                                gcgTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCtt                                 <  2:172611/151‑1 (MQ=255)
                                                                                                                gcgTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAAAACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCtt                                 >  1:558300/1‑151 (MQ=255)
                                                                                                                 cgTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCttt                                <  2:1123691/151‑1 (MQ=255)
                                                                                                                  gTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTg                               <  2:1185477/151‑1 (MQ=255)
                                                                                                                      gTAAGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAgcg                           >  2:442624/1‑151 (MQ=255)
                                                                                                                         aGCGCAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAgctc                          <  2:43939/149‑3 (MQ=255)
                                                                                                                           cgcAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTg                      <  2:592669/151‑1 (MQ=255)
                                                                                                                            gcAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTgg                     <  1:566112/151‑1 (MQ=255)
                                                                                                                            gcAGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTgg                     <  1:63984/151‑1 (MQ=255)
                                                                                                                              aGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTc                   >  1:851786/1‑151 (MQ=255)
                                                                                                                              aGCGGTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTc                   <  2:1008552/151‑1 (MQ=255)
                                                                                                                                 ggTGTTGTCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGc                <  2:1194062/151‑1 (MQ=255)
                                                                                                                                       gtCGCGTTCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTa           >  2:286510/1‑150 (MQ=255)
                                                                                                                                             ttCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTAATCTCCt    >  1:630347/1‑151 (MQ=255)
                                                                                                                                              tCCAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTAATCTCCt    <  1:625246/150‑1 (MQ=255)
                                                                                                                                               ccAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTAATCTCCTGc  >  1:240359/1‑151 (MQ=255)
                                                                                                                                               ccAGGGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTAATCTCCTGc  >  2:973608/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCAGGAATATCCGTATCGAAAATGATGGTACCTCCATTCATGGCCAGTCCACCAATGGATTGCTCTCCAACTTTTACCGATGTGGTATTTTCACTGTCAGACTGCAACATCGCGTGAGTAAGCGCAGCGGTGTTGTCGCGTTCCAGAGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGTAATCTCCTGC  >  NC_000913/1595513‑1595806

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