Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,294 G→C 6.2% D292H (GAT→CAT)  ygiC → putative enzyme with ATPase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2940GC 6.2% 377.2 / 8.6 112D292H (GAT→CAT) ygiCputative enzyme with ATPase activity
Reads supporting (aligned to +/- strand):  ref base G (52/52);  new base C (5/2);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 117 total aligned reads displayed.

AGATGTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTA  >  NC_000913/3181145‑3181443
                                                                                                                                                     |                                                                                                                                                     
agaTGTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatg                                                                                                                                                       >  1:151597/1‑150 (MQ=255)
 gaTGTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatga                                                                                                                                                      >  1:468346/1‑150 (MQ=255)
    gTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                                                                                   <  2:1323424/150‑1 (MQ=255)
     ttCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcat                                                                                                                                                  >  1:1628245/1‑150 (MQ=255)
               aGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                                                                        <  1:532381/150‑4 (MQ=255)
                 cTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataag                                                                                                                                      <  1:1016499/150‑6 (MQ=255)
                  tGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                                                                     <  2:466411/150‑1 (MQ=255)
                   ggaggaTGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaa                                                                                                                                    >  1:634538/1‑150 (MQ=255)
                         tGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                                                                              <  1:642987/150‑1 (MQ=255)
                         tGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                                                                              <  2:1181781/150‑1 (MQ=255)
                           cAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTgg                                                                                                                            >  2:1193637/1‑150 (MQ=255)
                           cAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTgg                                                                                                                            >  2:679795/1‑150 (MQ=255)
                             ggCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                                                                          >  1:548990/1‑150 (MQ=255)
                             ggCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                                                                          >  1:1570126/1‑150 (MQ=255)
                               cGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTaa                                                                                                                        >  1:687735/1‑150 (MQ=255)
                                gTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTaaa                                                                                                                       >  2:706934/1‑150 (MQ=255)
                                   cGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCg                                                                                                                    >  1:1010329/1‑150 (MQ=255)
                                   cGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCg                                                                                                                    >  1:13837/1‑150 (MQ=255)
                                      tggctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGATAAACCGAtc                                                                                                                 <  2:1129277/150‑1 (MQ=255)
                                       ggctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtct                                                                                                                <  1:577896/150‑1 (MQ=255)
                                        gctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtctt                                                                                                               >  1:546387/1‑150 (MQ=255)
                                         ctggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtcttc                                                                                                              <  1:206646/150‑1 (MQ=255)
                                           ggAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTtctc                                                                                                            >  1:607517/1‑150 (MQ=255)
                                            gAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                                                                           <  1:706934/150‑1 (MQ=255)
                                             aaCCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTccc                                                                                                          <  2:1010329/150‑1 (MQ=255)
                                                cGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                                                                       >  1:705268/1‑150 (MQ=255)
                                                 ggCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGa                                                                                                      <  2:1272466/150‑1 (MQ=255)
                                                       gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                                                                                <  2:40043/150‑1 (MQ=255)
                                                        aaGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                                                                               >  1:106424/1‑150 (MQ=255)
                                                          gagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaaa                                                                                             >  2:243593/1‑150 (MQ=255)
                                                            gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                                                                           >  1:1094508/1‑150 (MQ=255)
                                                              aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtg                                                                                         <  2:1518948/150‑1 (MQ=255)
                                                                tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                                                                       <  2:565663/150‑1 (MQ=255)
                                                                tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                                                                       <  2:1513157/150‑1 (MQ=255)
                                                                tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTc                                                                                       <  2:228907/150‑1 (MQ=255)
                                                                 aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                                                                      >  2:306786/1‑150 (MQ=255)
                                                                 aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                                                                      <  1:1113478/150‑1 (MQ=255)
                                                                 aTCTCCAACAAGGCACTTCTACAGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGAGGAAGATGATCATCCGCAAATGGAAAAATATGTGCTTAAACCGATCTTCTCGCGAGAAGGCGCAAACGTGTCg                                                                                      >  1:1015619/1‑150 (MQ=255)
                                                                    tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATc                                                                                   >  2:1160834/1‑150 (MQ=255)
                                                                     ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCa                                                                                  >  1:459553/1‑150 (MQ=255)
                                                                       aacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAtt                                                                                <  2:619842/150‑1 (MQ=255)
                                                                        acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                                                                               <  2:151597/150‑1 (MQ=255)
                                                                           aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                                                                            >  1:1484247/1‑150 (MQ=255)
                                                                            ggCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCGGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGaa                                                                           >  1:359902/1‑150 (MQ=255)
                                                                              cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                                                                         <  2:1484247/150‑1 (MQ=255)
                                                                               aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACgg                                                                        >  1:514429/1‑150 (MQ=255)
                                                                                cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                                                                       <  2:634538/150‑1 (MQ=255)
                                                                                     aCCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                                                                  >  1:708291/1‑150 (MQ=255)
                                                                                         cTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt                                                              >  1:1615921/1‑150 (MQ=255)
                                                                                          tACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg                                                             <  2:1628245/150‑1 (MQ=255)
                                                                                              gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagc                                                         >  2:1056327/1‑150 (MQ=255)
                                                                                                gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcag                                                       >  2:321862/1‑150 (MQ=255)
                                                                                                  gagaTGTTCCCGAATCACCCGAACCTGCTGCGCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg                                                     <  1:1098698/150‑1 (MQ=255)
                                                                                                  gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg                                                     <  1:1104286/150‑1 (MQ=255)
                                                                                                  gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg                                                     <  2:708291/150‑1 (MQ=255)
                                                                                                   agaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCgg                                                    >  1:395166/1‑150 (MQ=255)
                                                                                                      tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAg                                                 >  1:236628/1‑150 (MQ=255)
                                                                                                        ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAggg                                               >  1:459069/1‑149 (MQ=255)
                                                                                                        ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGt                                               >  1:95588/1‑150 (MQ=255)
                                                                                                         tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTc                                              >  2:1373770/1‑150 (MQ=255)
                                                                                                         tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTc                                              >  2:82699/1‑150 (MQ=255)
                                                                                                         tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTc                                              >  2:88725/1‑150 (MQ=255)
                                                                                                         tCCCGAATCACCCGAACCTGCTGCCCGCATATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGGTAAACCGATCTTCTCCCGTGAAGGCGCAAACGGGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTc                                              <  2:1620500/150‑1 (MQ=255)
                                                                                                          cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTcc                                             >  2:723523/1‑150 (MQ=255)
                                                                                                          cccGAATCACCAGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCACGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTcc                                             >  2:587076/1‑150 (MQ=255)
                                                                                                           ccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCg                                            >  1:1429432/1‑150 (MQ=255)
                                                                                                            cGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGt                                           >  1:739243/1‑150 (MQ=255)
                                                                                                              aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTAt                                         <  1:1193637/150‑1 (MQ=255)
                                                                                                               aTCACCCGACCCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATg                                        >  2:215376/1‑150 (MQ=255)
                                                                                                               aTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATg                                        <  1:29824/150‑1 (MQ=255)
                                                                                                                tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGCTCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATgg                                       >  2:814024/1‑150 (MQ=255)
                                                                                                                 cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATTATCATCCGCAAATGGAAAAAAATGGGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGACCGGCAAACCCATTGAAGCAGCGGAAGGTCCGCATGGc                                      <  2:748204/150‑1 (MQ=255)
                                                                                                                 cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGc                                      <  1:639466/150‑1 (MQ=255)
                                                                                                                  aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCg                                     >  1:547931/1‑150 (MQ=255)
                                                                                                                   cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg                                                     >  1:1364134/1‑133 (MQ=255)
                                                                                                                   cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCga                                    <  2:1094508/150‑1 (MQ=255)
                                                                                                                       aaCCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCgaagaa                                <  1:1567228/150‑1 (MQ=255)
                                                                                                                          ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAggg                             >  1:218712/1‑150 (MQ=255)
                                                                                                                            gctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGgat                           <  2:546387/150‑1 (MQ=255)
                                                                                                                             ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGgatg                          <  2:172687/150‑1 (MQ=255)
                                                                                                                               gcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGGTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGgatgat                        >  1:633596/1‑150 (MQ=255)
                                                                                                                               gcCCGCTTATTTTGCGGAAGATGATAATCCGCATATGGAAAAATATGTGGTTTAACCGATCTTCTCCCGTGAAGGCGCCAACGTGTCGATGATTGAGAACGGCAAAACAAGTGAGGCCGGGGAAGGTGCGGATGGGGAAGCCGGgatgct                        >  2:698302/1‑148 (MQ=255)
                                                                                                                                  cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGt                     >  1:139228/1‑150 (MQ=255)
                                                                                                                                  cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGt                     >  2:653710/1‑150 (MQ=255)
                                                                                                                                    cTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTc                   <  2:514429/150‑1 (MQ=255)
                                                                                                                                     ttATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCTTTGATAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTct                  <  2:1570126/150‑2 (MQ=255)
                                                                                                                                      ccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAg                 >  1:515693/4‑150 (MQ=255)
                                                                                                                                       aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGc                <  1:1160834/150‑1 (MQ=255)
                                                                                                                                       aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGc                >  2:1365799/1‑150 (MQ=255)
                                                                                                                                        atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCa               >  1:1234531/2‑150 (MQ=255)
                                                                                                                                        atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCa               >  1:1076379/2‑150 (MQ=255)
                                                                                                                                         tttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCaa              <  2:95588/150‑1 (MQ=255)
                                                                                                                                            gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAgg                              >  1:244513/1‑131 (MQ=255)
                                                                                                                                            gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAgg                              <  2:244513/131‑1 (MQ=255)
                                                                                                                                            gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTc           <  1:1513668/150‑1 (MQ=255)
                                                                                                                                            gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGGCGATCATTGAGAACGGCAAAACCATTGACGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTc           <  2:1220801/150‑1 (MQ=255)
                                                                                                                                                 aGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCAccc      >  1:1117509/1‑150 (MQ=255)
                                                                                                                                                    tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGtt   >  1:1062826/1‑150 (MQ=255)
                                                                                                                                                     gatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTa  >  1:1505976/1‑150 (MQ=255)
                                                                                                                                                     gatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTa  >  2:1335860/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGATGTTCTCAACCAAGCTGGAGGATGCAGGCGTACGCTGGCTGGAACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGGAAGGTCCGTATGGCGAAGAAGGGATGATTGTTCAGCAATTCCACCCGTTA  >  NC_000913/3181145‑3181443

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