Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,084,357 Δ1 bp coding (1555/1977 nt) speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,084,3570G.99.2% 575.1 / ‑2.2 120coding (1555/1977 nt)speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base . (61/57);  minor base C (0/1);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 120 total aligned reads displayed.

TGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACGTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCTCGTCGATAATCGGAC  >  NC_000913/3084211‑3084500
                                                                                                                                                  |                                                                                                                                               
tGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCg                                                                                                                                            >  2:1637183/1‑151 (MQ=255)
 gcgcCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCgg                                                                                                                                           >  1:243211/1‑151 (MQ=255)
    ccGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCAc                                                                                                                                        <  1:1442896/151‑1 (MQ=255)
    ccGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCAc                                                                                                                                        <  1:577385/151‑1 (MQ=255)
     cGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACt                                                                                                                                       >  1:1248081/1‑151 (MQ=255)
     cGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACt                                                                                                                                       >  2:1804396/1‑151 (MQ=255)
        ccATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTg                                                                                                                                     >  1:1820534/1‑150 (MQ=255)
        ccATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGa                                                                                                                                    <  1:1598407/151‑1 (MQ=255)
        ccATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGa                                                                                                                                    <  2:597765/151‑1 (MQ=255)
            aaaGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATcc                                                                                                                                 <  1:61125/150‑1 (MQ=255)
            aaaGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATcc                                                                                                                                 <  2:1089414/150‑1 (MQ=255)
             aaGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAg                                                                                                                               <  1:1672654/151‑1 (MQ=255)
                   ccGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTc                                                                                                                         >  1:1885485/1‑151 (MQ=255)
                   ccGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTc                                                                                                                         >  1:1112727/1‑151 (MQ=255)
                     gAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCa                                                                                                                       <  1:1040180/151‑1 (MQ=255)
                      aGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTcc                                                                                                                        >  2:1510369/1‑149 (MQ=255)
                      aGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTcc                                                                                                                        >  1:720605/1‑149 (MQ=255)
                        cATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCg                                                                                                                    <  1:180365/151‑1 (MQ=255)
                        cATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCg                                                                                                                    >  2:1811755/1‑151 (MQ=255)
                         aTCGGCGGATTCTCTGGATCGTACTCCGGCATTTGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCgg                                                                                                                   >  1:1095648/1‑151 (MQ=255)
                            ggcggATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAg                                                                                                                <  1:1518500/151‑1 (MQ=255)
                            ggcggATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAg                                                                                                                <  1:1824933/151‑1 (MQ=255)
                             gcggATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCGCAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGc                                                                                                                  >  2:264052/1‑148 (MQ=255)
                             gcggATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACGTTCCGGCACTTGATCCAGCCCTTCCAGCGGc                                                                                                                  >  2:614209/1‑149 (MQ=255)
                               ggATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAAc                                                                                                             >  2:525680/1‑151 (MQ=255)
                               ggATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAAc                                                                                                             >  2:431665/1‑151 (MQ=255)
                                gATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAAcc                                                                                                            <  1:1691001/151‑1 (MQ=255)
                                   tctcTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCggg                                                                                                         >  2:1605022/1‑151 (MQ=255)
                                   tctcTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCggg                                                                                                         >  1:98671/1‑151 (MQ=255)
                                      ctGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGa                                                                                                        <  1:1814379/149‑1 (MQ=255)
                                        ggATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGaa                                                                                                       <  1:1737890/148‑1 (MQ=255)
                                          aTCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACt                                                                                                  >  2:1533202/1‑151 (MQ=255)
                                                 tCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGa                                                                                            >  1:1550838/1‑150 (MQ=255)
                                                 tCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGAt                                                                                           >  1:1575777/1‑151 (MQ=255)
                                                  ccGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATc                                                                                          <  1:189029/151‑1 (MQ=255)
                                                    ggcattggcattGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATccc                                                                                        <  1:280788/151‑1 (MQ=255)
                                                     gcattggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATcccc                                                                                       >  1:516195/1‑151 (MQ=255)
                                                       attggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATccccc                                                                                      >  1:1933133/1‑150 (MQ=255)
                                                         tggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCAt                                                                                    <  1:755568/150‑1 (MQ=255)
                                                         tggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCAt                                                                                    <  1:139148/150‑1 (MQ=255)
                                                         tggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCAt                                                                                    <  2:456105/150‑1 (MQ=255)
                                                         tggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCAt                                                                                    <  2:1833028/150‑1 (MQ=255)
                                                          ggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGc                                                                                  >  2:708801/1‑151 (MQ=255)
                                                          ggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGc                                                                                  >  2:1356351/1‑151 (MQ=255)
                                                          ggcattgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGc                                                                                  <  1:1840652/151‑1 (MQ=255)
                                                               tgTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTcc                                                                              <  1:223412/150‑1 (MQ=255)
                                                                gTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCg                                                                             <  1:657522/150‑1 (MQ=255)
                                                                   gTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCAt                                                                         >  2:1232558/1‑151 (MQ=255)
                                                                     ggCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATc                                                                        <  2:1775771/150‑1 (MQ=255)
                                                                     ggCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATc                                                                        <  1:1746357/150‑1 (MQ=255)
                                                                     ggCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCg                                                                       <  1:1637454/151‑1 (MQ=255)
                                                                            ccGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGaa                                                                <  2:681176/151‑1 (MQ=255)
                                                                             cGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAAc                                                               >  1:79057/1‑151 (MQ=255)
                                                                                cACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGc                                                            <  2:1103092/151‑1 (MQ=255)
                                                                                cACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGc                                                            <  2:612052/151‑1 (MQ=255)
                                                                                cACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGc                                                            <  1:1705579/151‑1 (MQ=255)
                                                                                 aCCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCg                                                           >  1:1127118/1‑151 (MQ=255)
                                                                                  ccATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCg                                                           >  2:455278/1‑150 (MQ=255)
                                                                                   cATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCgag                                                         <  2:843936/151‑1 (MQ=255)
                                                                                   cATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCg                                                           >  1:1181361/1‑149 (MQ=255)
                                                                                    aTCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCg                                                           >  1:1291281/1‑148 (MQ=255)
                                                                                      cAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCAGCATCGACTGGAACAGCGAGAAg                                                      <  1:669838/151‑1 (MQ=255)
                                                                                       aaTATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGt                                                     >  1:598082/1‑151 (MQ=255)
                                                                                        atatAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGtt                                                    >  1:432598/1‑151 (MQ=255)
                                                                                        atatAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGtt                                                    >  1:534135/1‑151 (MQ=255)
                                                                                          ataGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGATGTTg                                                   >  2:1921746/1‑150 (MQ=255)
                                                                                          ataGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGa                                                  <  2:863021/151‑1 (MQ=255)
                                                                                          ataGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGa                                                  <  1:1223073/151‑1 (MQ=255)
                                                                                                  cGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACAtt                                          >  2:198898/1‑151 (MQ=255)
                                                                                                   gATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACAttt                                         <  2:565590/151‑1 (MQ=255)
                                                                                                    aTAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACAtttt                                        <  1:1357091/151‑1 (MQ=255)
                                                                                                    aTAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACAtttt                                        >  1:210249/1‑151 (MQ=255)
                                                                                                      aGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGt                                      >  2:827216/1‑151 (MQ=255)
                                                                                                       ggaCCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTACCTGCACTTGAGCCAGCACTTCCAGCGGCAGACACGTGAACAACTGGTCGACCCCCCATGCGTCCGGCATCGTCTGGAACAGCGAGAATTTGGCGTACATTTTGt                                      <  1:517250/148‑1 (MQ=255)
                                                                                                       gCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTc                                     >  1:1033604/1‑151 (MQ=255)
                                                                                                        cACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTcc                                    >  2:1651497/1‑151 (MQ=255)
                                                                                                         aCCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCg                                   >  2:766594/1‑151 (MQ=255)
                                                                                                           cGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCGCTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCCACTGGAACAGCGAGAAGTTGACGTACATTTTGtcctcc                                 <  1:566268/151‑1 (MQ=255)
                                                                                                           cGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGcc                                 >  1:1740767/1‑151 (MQ=255)
                                                                                                              cAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATa                              <  2:571659/151‑1 (MQ=255)
                                                                                                               agagTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATAc                             <  1:1871003/151‑1 (MQ=255)
                                                                                                               agagTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATAc                             <  1:1648630/151‑1 (MQ=255)
                                                                                                               agagTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATAc                             >  1:910110/1‑151 (MQ=255)
                                                                                                               agagTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAACCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATAc                             <  2:1921815/151‑1 (MQ=255)
                                                                                                                gagTCACAGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACg                            <  1:1621176/151‑1 (MQ=255)
                                                                                                                    cacaGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTcc                        >  2:382913/1‑151 (MQ=255)
                                                                                                                     acaGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCt                       >  2:608052/1‑151 (MQ=255)
                                                                                                                     acaGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCt                       >  2:1307377/1‑151 (MQ=255)
                                                                                                                     acaGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCt                       >  1:888342/1‑151 (MQ=255)
                                                                                                                       aGGTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGc                     <  1:1716760/151‑1 (MQ=255)
                                                                                                                        ggTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCa                    <  2:1778203/151‑1 (MQ=255)
                                                                                                                         gTAATATCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGc                     >  2:92634/1‑149 (MQ=255)
                                                                                                                              atCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCtcgt              <  1:1230901/151‑1 (MQ=255)
                                                                                                                               tCCAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCtcgtc             >  1:87705/1‑151 (MQ=255)
                                                                                                                                ccAGCAGCACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCTTGCAGCtcgtc             <  2:712615/150‑1 (MQ=255)
                                                                                                                                 cagcagCACAGCTCGCCTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCtcgtcg            <  1:1177693/150‑1 (MQ=255)
                                                                                                                                  agcagcACAGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGGGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCTCGTCGAt          >  1:1414539/1‑151 (MQ=255)
                                                                                                                                         caGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCTCGTCGATAATCGGa   <  2:1399022/151‑1 (MQ=255)
                                                                                                                                          aGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCTCGTCGATAATCGGAc  <  1:1534593/151‑1 (MQ=255)
                                                                                                                                          aGCGCGACTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCTCGTCGATAATCGGAc  >  1:424424/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGCGCCGACCATAAAGAAACCGAGCATCGGCGGATTCTCTGGATCGTACTCCGGCATTGGCATTGTCGTGGCAATACCGTCACCATCAATATAGTGGTCGATAGCACCGTCAGAGTCACAGGTAATATCCAGCAGCACAGCGCGACGTTCCGGCACTTGATCCAGCCCTTCCAGCGGCAGAACCGGGAACAACTGGTCGATCCCCCATGCGTCCGGCATCGACTGGAACAGCGAGAAGTTGACGTACATTTTGTCCGCCATACGTTCCTGCAGCTCGTCGATAATCGGAC  >  NC_000913/3084211‑3084500

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