Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,041,173 C→A 5.1% intergenic (+55/‑202) yedZ → / → zinT inner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,1730CA 5.1% 365.6 / 6.4 127intergenic (+55/‑202)yedZ/zinTinner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (72/47);  new base A (2/6);  total (74/53)
Fisher's exact test for biased strand distribution p-value = 6.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

GCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTAT  >  NC_000913/2041025‑2041321
                                                                                                                                                    |                                                                                                                                                    
gctggctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGCTCGCCCt                                                                                                                                                     >  2:103587/1‑150 (MQ=255)
gctggctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCt                                                                                                                                                     <  2:871674/150‑1 (MQ=255)
   ggctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATATCTTCATATTGCCGACAAAGTAAGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATc                                                                                                                                                  >  2:187650/1‑150 (MQ=255)
    gctgTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATcc                                                                                                                                                 <  1:1267950/150‑1 (MQ=255)
         aCTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCAt                                                                                                                                            >  2:1358495/1‑150 (MQ=255)
         aCTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCag                                                                                                                                            <  1:1784178/150‑2 (MQ=255)
         aCTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTAAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCag                                                                                                                                            <  1:1897300/150‑2 (MQ=255)
          cTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATa                                                                                                                                           <  2:1318358/150‑1 (MQ=255)
          cTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGCagg                                                                                                                                           <  1:1140997/150‑3 (MQ=255)
           tGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCTCGCGCCACATCCGGCaggg                                                                                                                                          <  1:31169/150‑4 (MQ=255)
           tGCTTTTAGCCTTACGGTATAAGAAGTTGCGATCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTACCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATac                                                                                                                                          <  2:1338693/150‑1 (MQ=255)
            gCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATaca                                                                                                                                         >  2:1469418/1‑150 (MQ=255)
              ttttAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACACAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATAcaca                                                                                                                                       <  1:87567/150‑1 (MQ=255)
              ttttAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGGACGCGTTATCCCTGCCGGATGTGGAGCGCGCGCCCTAGCCGGCATAcaca                                                                                                                                       >  2:1378231/1‑150 (MQ=255)
              tattAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATAcaca                                                                                                                                       <  1:269626/148‑1 (MQ=255)
               tttAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAg                                                                                                                                      >  1:1970421/1‑150 (MQ=255)
                   gCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCTGCATACACAGATCg                                                                                                                                  >  1:269776/1‑150 (MQ=255)
                        aCGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCa                                                                                                                             <  2:1706512/150‑1 (MQ=255)
                         cGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGAAAGCAGGTTCACAACAAATTATCTGAATAACTTCATATTGCCGACAAAGTACGCGTGATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg                                                                                                                            >  2:1714768/1‑150 (MQ=255)
                          ggTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAgg                                                                                                                           >  1:1441342/1‑150 (MQ=255)
                           gTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGt                                                                                                                          >  1:1332559/1‑150 (MQ=255)
                             ataAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTa                                                                                                                        >  1:915776/1‑150 (MQ=255)
                              taAGAAGTTGCGTACCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGACGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGAACGCATCAGGatat                                                                                                                       >  2:1033444/1‑150 (MQ=255)
                               aagaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCCCAACAAATTATCTGTATAACTTCATATTGCCGCCAAAGTTCGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCCGCATACACAGATCGCATCAGCTTATg                                                                                                                      >  2:942430/1‑150 (MQ=255)
                                agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGATAACAAAAAATTATCTGTATAAATTGATATTGCATACAAAGAACGCGTTATCCCTGCCGGATGTGGCGCTAGCGCCCTATCCGGCAAACACAGATCGCATCAGGTTATGt                                                                                                                     >  2:1583289/1‑150 (MQ=255)
                                 gaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGtt                                                                                                                    >  2:1638253/1‑150 (MQ=255)
                                      tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCgat                                                                                                               >  2:899411/1‑148 (MQ=255)
                                        cGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTg                                                                                                             <  1:595798/150‑1 (MQ=255)
                                        cGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTg                                                                                                             >  1:755414/1‑150 (MQ=255)
                                        cGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGGTTCGTTTg                                                                                                             >  1:789817/1‑150 (MQ=255)
                                        cGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAAATATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCTCGAGCGCCCTATCCGGCATACACAGATCGCAACAGGTTATGTTTCGTTTg                                                                                                             >  2:669849/1‑150 (MQ=255)
                                           tCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtca                                                                                                          <  2:938855/150‑1 (MQ=255)
                                            cccTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCTGATGTGGCGCGAGCGCCCTATCCGGCATACACAGCTCGCATCAGGTTATGTTACGTTTGtcat                                                                                                         >  1:90954/1‑150 (MQ=255)
                                                   tAACCGGTTACGCATGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTc                                                                                                  <  2:1994653/150‑1 (MQ=255)
                                                   tAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCAGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTGTCGTTTGTCATCAACCTc                                                                                                  >  2:1628884/1‑150 (MQ=255)
                                                      ccGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaaa                                                                                               >  2:1423109/1‑150 (MQ=255)
                                                       cGGTTACGCAAGCAGGTTCACAACAAATCATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAt                                                                                              >  2:1884245/1‑150 (MQ=255)
                                                        ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTACATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGAGGCTCGAGCGCCATATCCGACATACACAGAACCCATCAGGATATCTTTCTTTTGTCATCAACCTCaaagc                                                                                             >  1:360330/1‑148 (MQ=255)
                                                               cAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCAACAACCTCAAATCACCaaaa                                                                                      >  1:1504966/1‑150 (MQ=255)
                                                                 tgCAGTTTCACAACAATTTATCTGTATATCTTCATATTGCCGACAATGTACGAGTTATCACTGCCGGATGTGGCGCGAGCTCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaata                                                                                    <  1:128750/149‑1 (MQ=255)
                                                                  tcAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAaataa                                                                                   <  1:1878341/149‑1 (MQ=255)
                                                                     ggTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTACCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAattt                                                                                >  2:649266/1‑150 (MQ=255)
                                                                       ttCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGTTGTGGCGCGAGCGCCCTATCCGGCATAAACAGATCGCATCATGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttat                                                                              <  1:1948695/150‑1 (MQ=255)
                                                                       ttCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAatttat                                                                              <  2:1441342/150‑1 (MQ=255)
                                                                         cacaACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTACGTTTGTCATCAACCTCAAATCACCAAAATAatttattt                                                                            >  1:977719/1‑150 (MQ=255)
                                                                          acaacaAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttatttt                                                                           >  2:518397/1‑150 (MQ=255)
                                                                           caacaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAAtttatttta                                                                          >  1:810471/1‑150 (MQ=255)
                                                                              caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTtattttattt                                                                       >  1:504752/1‑150 (MQ=255)
                                                                              caaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTtattttattt                                                                       >  1:1244547/1‑150 (MQ=255)
                                                                                 ataaTCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGGGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCaa                                                                    <  2:848333/147‑1 (MQ=255)
                                                                                 aTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCaa                                                                    <  2:842137/150‑1 (MQ=255)
                                                                                     tCTGTATAACTTCATATTTCCTACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCTTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAag                                                                <  1:1263363/150‑1 (MQ=255)
                                                                                      cTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAaga                                                               <  2:1932894/150‑1 (MQ=255)
                                                                                       tGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagat                                                              <  2:1332559/150‑1 (MQ=255)
                                                                                        gTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagata                                                             >  2:1497935/1‑150 (MQ=255)
                                                                                         tataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagatag                                                            >  2:1425025/1‑150 (MQ=255)
                                                                                          ataACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAagataga                                                           >  1:38947/1‑150 (MQ=255)
                                                                                              cTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATaat                                                       >  1:1521377/1‑149 (MQ=255)
                                                                                               ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGt                                                      >  1:1174427/1‑150 (MQ=255)
                                                                                                 cATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAg                                                    <  2:784429/150‑1 (MQ=255)
                                                                                                   tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCCCCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGaa                                                  >  2:2039358/1‑150 (MQ=255)
                                                                                                    tttGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGTCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAAc                                                 <  1:1865308/149‑1 (MQ=255)
                                                                                                    atTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGaaa                                                 >  2:1921994/1‑149 (MQ=255)
                                                                                                      tGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGAGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTGATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTg                                               >  2:347648/1‑150 (MQ=255)
                                                                                                      tGCCGACAAAGTACGCGTAATCCCTGCCGGATGTTGCGCGAGCGCCATATCCGGCATACACAGATCGCATCAGGTTATGTTTAGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTg                                               <  1:1404320/150‑1 (MQ=255)
                                                                                                       gCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGAATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACAAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTga                                              >  1:31399/1‑150 (MQ=255)
                                                                                                        ccGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCCTACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCATTAAGATAGATAAGTAGAACTgag                                             >  1:1899521/1‑150 (MQ=255)
                                                                                                        ccGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTgag                                             >  1:563639/1‑150 (MQ=255)
                                                                                                        ccGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGAGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTgag                                             <  1:1023936/150‑1 (MQ=255)
                                                                                                           aCAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCACCCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGACCTGAGaaa                                          >  1:2033773/1‑150 (MQ=255)
                                                                                                               aGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCa                                      >  1:1781451/1‑150 (MQ=255)
                                                                                                               aGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCa                                      <  1:1077552/150‑1 (MQ=255)
                                                                                                               aGTACGCGTTATCCCTACCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCa                                      <  1:649266/150‑1 (MQ=255)
                                                                                                                 tACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATATGATAGATAAGTAGAACTGAGAAAGCCATg                                    >  2:732548/1‑150 (MQ=255)
                                                                                                                      gTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCtctc                               <  1:1371736/150‑1 (MQ=255)
                                                                                                                       ttATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCg                              >  1:1673435/1‑150 (MQ=255)
                                                                                                                          tCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGttt                           <  2:1837636/150‑1 (MQ=255)
                                                                                                                          tCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCAGCACCCTCAAATCACCAAAATAATTTATTTTAGTTCAATAAGATAGATAAGTTGATCTGAGAAAGCCATGCTCTCGttt                           >  1:392049/1‑150 (MQ=255)
                                                                                                                           cccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTc                          >  2:1941445/1‑150 (MQ=255)
                                                                                                                            ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACTCAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTcc                         <  2:678368/150‑1 (MQ=255)
                                                                                                                            ccTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTcc                         <  1:1333237/150‑1 (MQ=255)
                                                                                                                                ccGGATGTTGCGCGATCGCCATATCAGGCATACACAGATCGCATCAGGTTTTGTTTCGTTTGTCATCAACCTCAAATCACCAAAATATTTTATTTTTTTTCAATAAGATAGATAAGTAGTACTGAGAAAGCCATGCTCTCGTTTCCTaat                     <  1:1033444/150‑2 (MQ=255)
                                                                                                                                ccGGATGTGGCGCGAGCGCCCTATCGGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAag                     <  2:428690/150‑1 (MQ=255)
                                                                                                                                ccGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAag                     <  1:1921312/150‑1 (MQ=255)
                                                                                                                                 cGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGAATAAGGTTATTTTTCGTTTGTCATCAACCTCAAAACTCCAAAATAATTTATTTAAAATAAATAAGATAGATAAGTAGAACTGAGATAGCCATGCTCTCGTTTCCTAcga                    >  2:1379344/1‑150 (MQ=255)
                                                                                                                                  ggATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAagag                   >  1:1586146/1‑150 (MQ=255)
                                                                                                                                  ggATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAagag                   >  1:110293/1‑150 (MQ=255)
                                                                                                                                        ggagcgAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTAAACTAAGAgttgttt             >  2:1433007/4‑149 (MQ=255)
                                                                                                                                        ggCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATGAGCCTCATATCACCAAAATAATTTATTTTATTTCATTATGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGttgttg             >  1:2058625/1‑150 (MQ=255)
                                                                                                                                         gcgcGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGc            >  1:2043847/1‑150 (MQ=255)
                                                                                                                                          cgcgAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCa           >  1:55346/1‑150 (MQ=255)
                                                                                                                                              aGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCAtttt       >  1:1586732/1‑150 (MQ=255)
                                                                                                                                                cgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAATATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTTAGAGTTGTTGCATTTTGc     <  2:507539/150‑1 (MQ=255)
                                                                                                                                                cgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGc     <  1:2082816/150‑1 (MQ=255)
                                                                                                                                                cgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGc     <  1:1951678/150‑1 (MQ=255)
                                                                                                                                                cgcCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGATGTTGCATTTTga     >  1:1410698/1‑149 (MQ=255)
                                                                                                                                                  cccTATCCGGCATACACAGCTCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGCTAAGTAGAACTGAGAAAGCCATACTCTCGTTTCCTAAGAGTTGATGCATTTTGAta   >  1:111465/1‑150 (MQ=255)
                                                                                                                                                  cccTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGATTGTCATCAGCCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCta   >  2:1818735/1‑150 (MQ=255)
                                                                                                                                                   ccTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCtat  >  2:752357/1‑150 (MQ=255)
                                                                                                                                                   ccTATCCGGCATACACAGAACGCATCAGGTGATGTTGCGTTTGTCATCAAACTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCtat  <  2:169250/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTAT  >  NC_000913/2041025‑2041321

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