Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 258,954 T→A 7.7% V163D (GTC→GAC)  PP_0207 → nitrate ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947258,9540TA 7.7% 347.0 / 15.1 115V163D (GTC→GAC) PP_0207nitrate ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (51/55);  new base A (4/5);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 115 total aligned reads displayed.

ATCAGTGTGCTGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGTGGCTTCG  >  NC_002947/258806‑259100
                                                                                                                                                    |                                                                                                                                                   
atcaGTGTGCTGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCACGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCCCGGCATGCTGTTGCGCGCCGTc                                                                                                                                                    >  2:252625/1‑150 (MQ=255)
atcaGTGTGCTGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTc                                                                                                                                                    >  1:2174925/1‑150 (MQ=255)
atcaGTGTGCTGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTc                                                                                                                                                    >  2:1606823/1‑150 (MQ=255)
       tgCTGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCgcgc                                                                                                                                             >  1:1973764/1‑150 (MQ=255)
          tGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGGCGACGGCGCGCaac                                                                                                                                          <  2:1056920/150‑3 (MQ=255)
          tGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAgg                                                                                                                                          >  2:1322589/1‑150 (MQ=255)
          tGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAgg                                                                                                                                          <  2:1902647/150‑1 (MQ=255)
          tGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAgg                                                                                                                                          <  1:598260/150‑1 (MQ=255)
           gTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGGCGACGGCGCGCaaca                                                                                                                                         <  2:1808151/150‑4 (MQ=255)
             cGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCt                                                                                                                                       <  2:1381082/150‑1 (MQ=255)
             cGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCt                                                                                                                                       <  2:542645/150‑1 (MQ=255)
              gggCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTg                                                                                                                                      <  2:103194/150‑1 (MQ=255)
               ggCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTgg                                                                                                                                     <  2:2161843/150‑1 (MQ=255)
               ggCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTgg                                                                                                                                     <  2:474356/150‑1 (MQ=255)
                           ggCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAg                                                                                                                         <  2:2082985/150‑1 (MQ=255)
                                    agcGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCgg                                                                                                                >  1:1668622/3‑150 (MQ=255)
                                    aaCGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCgg                                                                                                                >  2:2153107/1‑150 (MQ=255)
                                    aaCGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCgg                                                                                                                <  1:2119428/150‑1 (MQ=255)
                                     aCGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGa                                                                                                               >  1:670138/1‑150 (MQ=255)
                                     aCGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGa                                                                                                               <  2:777008/150‑1 (MQ=255)
                                     aCGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCGATTTGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGa                                                                                                               <  1:1833473/150‑1 (MQ=255)
                                        gCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGatca                                                                                                            <  1:2007199/150‑1 (MQ=255)
                                          aTCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGatcatc                                                                                                          >  1:1595034/1‑150 (MQ=255)
                                               cgtGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCa                                                                                                     <  2:1617850/150‑1 (MQ=255)
                                                gtGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTgg                                                                                                                                     >  1:2403660/1‑117 (MQ=255)
                                                gtGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTgg                                                                                                                                     <  2:2403660/117‑1 (MQ=255)
                                                gtGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAg                                                                                                    <  2:670138/150‑1 (MQ=255)
                                                gtGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAg                                                                                                    <  1:1606823/150‑1 (MQ=255)
                                                 tGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGAccc                                                                                                                                <  2:1959667/121‑1 (MQ=255)
                                                     ctccgccTCCAGCGTGCAGTCGCTGGCGGAGCTTAAGGCCAAGACCTTTTCCGTGCCGTTCGCCTCCTCAGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGcc                                                                                               <  2:1706200/148‑1 (MQ=255)
                                                       ccgccTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCa                                                                                                     >  1:458695/1‑142 (MQ=255)
                                                       ccgccTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCa                                                                                                     <  2:458695/142‑1 (MQ=255)
                                                         gccTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGaa                                                                                           <  1:2147782/150‑1 (MQ=255)
                                                            tCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAGATc                                                                                        >  2:1592372/1‑150 (MQ=255)
                                                            tCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATc                                                                                        >  2:1699831/1‑150 (MQ=255)
                                                                          gCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCt                                                                         >  1:1639635/1‑150 (MQ=255)
                                                                          gCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCt                                                                         <  1:2153107/150‑1 (MQ=255)
                                                                           cTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTg                                                                        >  1:160935/1‑150 (MQ=255)
                                                                            tGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTgc                                                                       <  1:253759/150‑1 (MQ=255)
                                                                              gCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAg                                                                                                                         >  1:18115/1‑99 (MQ=255)
                                                                              gCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAg                                                                                                                         <  2:18115/99‑1 (MQ=255)
                                                                               cGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGACGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCgca                                                                    <  1:163224/150‑1 (MQ=255)
                                                                                   gCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCaa                                                                <  1:433910/150‑1 (MQ=255)
                                                                                    cTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCAACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAAc                                                               <  1:1123227/150‑1 (MQ=255)
                                                                                       aaGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTgg                                                                                                                                     <  2:642217/78‑1 (MQ=255)
                                                                                       aaGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTgg                                                                                                                                     >  1:642217/1‑78 (MQ=255)
                                                                                       aaGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGc                                                            >  1:1651740/1‑150 (MQ=255)
                                                                                        aGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGGCGACGGCGCGCaaca                                                                                                                                         <  1:191371/73‑4 (MQ=37)
                                                                                        aGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCa                                                           >  2:1987027/1‑150 (MQ=255)
                                                                                          ggCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATc                                                         <  2:1639635/150‑1 (MQ=255)
                                                                                          ggCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATc                                                         <  2:160935/150‑1 (MQ=255)
                                                                                             aaGACCATTTCCGTTCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGaa                                                      <  1:2176399/150‑1 (MQ=255)
                                                                                             aaGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGaa                                                      <  1:597308/150‑1 (MQ=255)
                                                                                               gACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGaa                                                      >  2:187328/1‑148 (MQ=255)
                                                                                               gACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGaa                                                      <  1:187328/148‑1 (MQ=255)
                                                                                                aCCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCgg                                                                                                                <  2:832586/90‑1 (MQ=255)
                                                                                                aCCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGcc                                                   >  2:1629588/1‑150 (MQ=255)
                                                                                                 ccATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGccc                                                  <  1:131802/150‑1 (MQ=255)
                                                                                                   aTTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGc                                                                                                                              >  2:1048898/1‑73 (MQ=255)
                                                                                                   aTTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGc                                                                                                                              <  1:1048898/73‑1 (MQ=255)
                                                                                                   aTTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATc                                                         >  2:466625/1‑141 (MQ=255)
                                                                                                   aTTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATc                                                         <  1:466625/141‑1 (MQ=255)
                                                                                                         gTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGAc                                          <  2:1668622/150‑1 (MQ=255)
                                                                                                          tGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACt                                         >  2:2407463/1‑150 (MQ=255)
                                                                                                               ttCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTg                                    >  2:1115416/1‑150 (MQ=255)
                                                                                                                   ccTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGcc                                                   >  1:2191267/1‑131 (MQ=255)
                                                                                                                   ccTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGcc                                                   <  2:2191267/131‑1 (MQ=255)
                                                                                                                    cTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGtt                               <  2:1888157/150‑1 (MQ=255)
                                                                                                                     tCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAg                                                     <  2:454538/127‑1 (MQ=255)
                                                                                                                     tCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAg                                                     >  1:454538/1‑127 (MQ=255)
                                                                                                                     tCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAAACGCGCGGCATTCCCAGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTc                              >  2:1268930/1‑150 (MQ=255)
                                                                                                                        aCTGCCCACGGCATGCTGTTGCGCGCCGTCGCCTCGCAGGGCTGGGACCCGCATAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCAGGCTCCGCCCTGCGCAGCAACCACATCGAAGGTCACGCCGACTTCGTGCCTTTCGcc                           >  1:754471/1‑150 (MQ=255)
                                                                                                                        aCTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGcc                           >  2:1423279/1‑150 (MQ=255)
                                                                                                                         cTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCg                          >  2:1922725/1‑150 (MQ=255)
                                                                                                                             ccACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAg                        <  2:1358165/148‑1 (MQ=255)
                                                                                                                             ccACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAg                        >  1:1358165/1‑148 (MQ=255)
                                                                                                                              cACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGatca                                                                                                            <  1:1152362/64‑1 (MQ=255)
                                                                                                                              cACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGatca                                                                                                            >  2:1152362/1‑64 (MQ=255)
                                                                                                                              cACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTa                     <  2:1973764/150‑1 (MQ=255)
                                                                                                                                 ggCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACt                                         <  2:1913776/127‑1 (MQ=255)
                                                                                                                                 ggCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACt                                         >  1:1913776/1‑127 (MQ=255)
                                                                                                                                 ggCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTc                  <  1:1922725/150‑1 (MQ=255)
                                                                                                                                 ggCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTc                  <  1:1987027/150‑1 (MQ=255)
                                                                                                                                   cATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGaa                                                                                                                          >  1:669847/1‑45 (MQ=255)
                                                                                                                                   cATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGaa                                                                                                                          <  2:669847/45‑1 (MQ=255)
                                                                                                                                     cagccctGCGCGGCGACGGCGCGCAAGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCa              >  1:705846/7‑150 (MQ=255)
                                                                                                                                        tGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAAcc           >  1:363912/1‑150 (MQ=255)
                                                                                                                                         gTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGt                       <  1:2151853/137‑1 (MQ=255)
                                                                                                                                         gTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGt                       >  2:2151853/1‑137 (MQ=255)
                                                                                                                                         gTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCg          >  1:555038/1‑150 (MQ=255)
                                                                                                                                         gTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCg          >  1:1516137/1‑150 (MQ=255)
                                                                                                                                         cctGCGCGGCGACGGCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCAAGCCGACTTCGTGCCGTTCGCCGAGTGAGTCCCCAACCg          >  2:522579/3‑150 (MQ=255)
                                                                                                                                          ttGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGt         >  2:380163/1‑150 (MQ=255)
                                                                                                                                          ttGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGt         <  1:252625/150‑1 (MQ=255)
                                                                                                                                          ttGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGt         >  1:879182/1‑150 (MQ=255)
                                                                                                                                          ttGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGt         >  1:905261/1‑150 (MQ=255)
                                                                                                                                          ctGCGCGGCGACGGCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGcc                                             >  2:1094607/2‑114 (MQ=255)
                                                                                                                                          ctGCGCGGCGACGGCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGcc                                             <  1:1094607/113‑1 (MQ=255)
                                                                                                                                           tGCGCGGCGACGGCGCGCAAGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGTg        <  2:705846/150‑1 (MQ=255)
                                                                                                                                           tGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGTg        <  1:248969/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCAGTGTGCTGTCGGGCAGCGTGCATGGCAGCGGCAACGGCATCGTCGTGCCCGCCGCCTCCAGCGTGCAGTCGCTGGCGGAGCTCAAGGGCAAGACCATTTCCGTGCCGTTCGCCTCCACTGCCCACGGCATGCTGTTGCGCGCCGTCGCCGCGCAGGGCTGGGACCCGCAGAAGGATGTGCGGATCATCGCCCAGGCGCCGGAAATCGCCGGCTCCGCCCTGCGCAGCAACCGCATCGAAGCCCACGCCGACTTCGTGCCGTTCGCCGAGTTATTCCCCAACCGTGGCTTCG  >  NC_002947/258806‑259100

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