Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 870,234 A→C 7.5% A64A (GCT→GCG)  PP_5456 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947870,2340AC 7.5% 251.6 / 13.6 93A64A (GCT→GCG) PP_5456hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/42);  new base C (5/2);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

GGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGTT  >  NC_002947/870087‑870383
                                                                                                                                                   |                                                                                                                                                     
ggCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAgc                                                                                                                                                     >  1:1060049/1‑150 (MQ=255)
 gCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCAATTCTTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGGTTGCGCGAgcg                                                                                                                                                    >  1:1677867/1‑150 (MQ=255)
   gAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCa                                                                                                                                                  >  1:2299841/1‑150 (MQ=255)
        cGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGCCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAAccc                                                                                                                                             <  2:911807/150‑1 (MQ=255)
              ccGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCtt                                                                                                                                       >  2:2142851/1‑150 (MQ=255)
               cGATCAGGATACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGACATTCGTGGTTTTTATTGAAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGAGAGCGCAAACGAGGGCttt                                                                                                                                      >  1:940055/1‑149 (MQ=255)
                   cAGGAACCGTCCTTAGAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTcggcg                                                                                                                                  <  2:1812188/150‑1 (MQ=255)
                    atgataCGTTCTTATCATCGGCTGTATGTTATAGCGATGTTGCGAGCAGAAATGTAGCTATTCGTGGTTTTTATTTCACCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCTAACGTCGGTTTGCGCGAGCGCAATCCCGGGCTTctgcgg                                                                                                                                 <  2:20498/145‑1 (MQ=255)
                    aGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTcggcgg                                                                                                                                 <  1:1269104/150‑1 (MQ=255)
                      gagaCGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTGGCCATTAGTGGTTTTTATTTCAGCGTTTCTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCgggg                                                                                                                               >  1:288621/1‑149 (MQ=255)
                           gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGa                                                                                                                                                       <  1:962445/121‑1 (MQ=255)
                           gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGa                                                                                                                                                       >  2:962445/1‑121 (MQ=255)
                           gTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTc                                                                                                                          <  2:1121719/150‑1 (MQ=255)
                              cTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTcgcg                                                                                                                       <  1:2142851/150‑1 (MQ=255)
                                     tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                           >  1:121200/1‑123 (MQ=255)
                                     tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCgg                                                                                                                                           <  2:121200/123‑1 (MQ=255)
                                     tCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagc                                                                                                                >  2:939137/1‑150 (MQ=255)
                                        gggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCccagccag                                                                                                             <  2:1849623/150‑1 (MQ=255)
                                         ggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                            <  2:1027950/118‑1 (MQ=255)
                                         ggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCg                                                                                                                                            >  1:1027950/1‑118 (MQ=255)
                                         ggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccagg                                                                                                            >  1:887180/1‑150 (MQ=255)
                                         ggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccagg                                                                                                            <  2:1060049/150‑1 (MQ=255)
                                         ggggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccagg                                                                                                            <  2:1402389/150‑1 (MQ=255)
                                           ggAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCTCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggcc                                                                                                          >  1:971574/1‑150 (MQ=255)
                                            gagaTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggcca                                                                                                         >  2:1966794/1‑150 (MQ=255)
                                               aTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCggg                                                                                                                                          >  1:2201905/1‑114 (MQ=255)
                                               aTTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCggg                                                                                                                                          <  2:2201905/114‑1 (MQ=255)
                                                  ataGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGaa                                                                                                   <  1:926751/150‑1 (MQ=255)
                                                       gATGTCGCGGGCAGGAATGTCGCCACTAGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTAAATTTCg                                                                                              >  2:8232/1‑150 (MQ=255)
                                                       gATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCg                                                                                              >  2:117760/1‑150 (MQ=255)
                                                        aTGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGc                                                                                             <  2:2056901/150‑1 (MQ=255)
                                                             gcgGGCAGAAATGTCGCCATTCGTGGTTTTTATTGAAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGGATTTCGCTGGTc                                                                                        >  1:846794/1‑150 (MQ=255)
                                                                ggCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGc                                                                                     <  1:1273096/150‑1 (MQ=255)
                                                                  cAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAc                                                                                   >  1:1200332/1‑150 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCTCGCGCAAAccgcc                                                                                                                                           <  2:1901137/93‑4 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCTCGCGCAAAccgcc                                                                                                                                           >  1:1901137/1‑90 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCa                                                                                                                              >  2:811297/1‑106 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCa                                                                                                                              <  1:811297/106‑1 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAcc                                                                                  >  2:386968/1‑150 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAcc                                                                                  >  1:770143/1‑150 (MQ=255)
                                                                   aGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGAAAGGAAGCGATGAGACCATAGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCAcc                                                                                  >  2:672252/1‑150 (MQ=255)
                                                                         gTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCg                                                                            >  2:1659360/1‑150 (MQ=255)
                                                                           cGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTTCTAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGAt                                                                          <  2:654295/150‑1 (MQ=255)
                                                                             ccATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCTCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGc                                                                        >  1:302680/1‑150 (MQ=255)
                                                                             ccATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggc                                                                                                           <  1:385152/115‑1 (MQ=255)
                                                                             ccATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAgccaggc                                                                                                           >  2:385152/1‑115 (MQ=255)
                                                                             ccATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGc                                                                        >  2:2023179/1‑150 (MQ=255)
                                                                                   gTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  >  2:1993848/1‑150 (MQ=255)
                                                                                            aTTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTgg                                                          <  2:333244/149‑1 (MQ=255)
                                                                                            aTTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTgg                                                          >  1:333244/1‑149 (MQ=255)
                                                                                             ttGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCt                                                        >  1:1475769/1‑150 (MQ=255)
                                                                                                cAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCTCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCa                                                     <  2:971574/150‑1 (MQ=255)
                                                                                                cAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCcca                                                                                                                  <  2:474755/89‑1 (MQ=255)
                                                                                                cAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCcca                                                                                                                  >  1:474755/1‑89 (MQ=255)
                                                                                                  gCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCAtt                                                   >  1:1161194/1‑150 (MQ=255)
                                                                                                     tttGTCAGGCAGCGAGGAGGCTCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCt                                                <  2:302680/150‑1 (MQ=255)
                                                                                                     tttGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCt                                                <  1:939137/150‑1 (MQ=255)
                                                                                                     tttGGCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCATGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCACGTGGTGCCTGGCTGCATTGCt                                                <  2:361903/150‑1 (MQ=255)
                                                                                                        gTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGa                                             <  2:2299841/150‑1 (MQ=255)
                                                                                                         tCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCATGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCACCGTGGTGCCTGGCTGCATTGCTTGaa                                            <  2:1677867/150‑1 (MQ=255)
                                                                                                           aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAg                                                                               <  2:1489691/113‑1 (MQ=255)
                                                                                                           aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAg                                                                               >  1:1489691/1‑113 (MQ=255)
                                                                                                           aGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAgg                                          >  2:1054843/1‑150 (MQ=255)
                                                                                                             gCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTAGCCCCAGGTCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGcc                                        >  1:143791/1‑150 (MQ=255)
                                                                                                              cAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCTGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGTTGCCTGGCTGCATTGCTTGAAGGCCa                                       <  2:940055/150‑1 (MQ=255)
                                                                                                              cAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                       <  2:1161194/150‑1 (MQ=255)
                                                                                                              cAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCa                                       <  1:758228/150‑1 (MQ=255)
                                                                                                               aGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAATCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                      <  2:692539/150‑1 (MQ=255)
                                                                                                               aGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                      >  1:2059611/1‑150 (MQ=255)
                                                                                                               aGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAg                                      >  1:539744/1‑150 (MQ=255)
                                                                                                                gCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAgg                                     >  2:1252692/1‑150 (MQ=255)
                                                                                                                   aggaggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCACCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCAt                                  >  1:145971/1‑150 (MQ=255)
                                                                                                                      aggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTg                               <  2:539744/150‑1 (MQ=255)
                                                                                                                       ggCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTTCTTGAAGGCCAGGCATTTGc                              >  2:544084/1‑150 (MQ=255)
                                                                                                                         cgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggt                                                                <  1:465217/114‑1 (MQ=255)
                                                                                                                         cgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtggt                                                                >  2:465217/1‑114 (MQ=255)
                                                                                                                         cgcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCa                            <  1:1993848/150‑1 (MQ=255)
                                                                                                                          gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTg                               <  1:921641/146‑1 (MQ=255)
                                                                                                                          gcTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTg                               >  2:921641/1‑146 (MQ=255)
                                                                                                                                   aaCGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCaa                  <  1:1054843/150‑1 (MQ=255)
                                                                                                                                      ggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                            >  2:232515/1‑121 (MQ=255)
                                                                                                                                      ggcggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                            <  1:232515/121‑1 (MQ=255)
                                                                                                                                      ccgggTTTGCGCTCGCGCAAACCGGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  <  2:301075/96‑1 (MQ=255)
                                                                                                                                      ccgggTTTGCGCTCGCGCAAACCGGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCgtg                                                                  >  1:301075/4‑99 (MQ=255)
                                                                                                                                      ccgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGc               >  2:1196294/4‑150 (MQ=255)
                                                                                                                                       cgggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCg              >  1:1953985/3‑150 (MQ=255)
                                                                                                                                        gggTTTGCGCTCGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             >  2:1031391/2‑150 (MQ=255)
                                                                                                                                        cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  1:86826/150‑1 (MQ=255)
                                                                                                                                        cggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCgg             <  1:1729476/150‑1 (MQ=255)
                                                                                                                                         ggTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCggg            <  1:386968/150‑1 (MQ=255)
                                                                                                                                           tttGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                            <  2:1930282/116‑1 (MQ=255)
                                                                                                                                           tttGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGaa                                            >  1:1930282/1‑116 (MQ=255)
                                                                                                                                            ttGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCAt         >  1:1537916/1‑150 (MQ=255)
                                                                                                                                                  gAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                           >  1:891852/1‑110 (MQ=255)
                                                                                                                                                  gAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAg                                           <  2:891852/110‑1 (MQ=255)
                                                                                                                                                   aGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGcagt  <  1:672252/150‑4 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCGAGTGCGGGGCCCGATCAGGAGACGTCCTTATAATCGGGGGGAGATTATAGCGATGTCGCGGGCAGAAATGTCGCCATTCGTGGTTTTTATTGCAGCGTTTGTCAGGCAGCGAGGAGGCGCTGGTGCCAACGGCGGTTTGCGCGAGCGCAAACCCGGGCTTCGGCGGGCAGCTCGCGCCCCAGCCAGGCCAGCTGAATTTCGCTGATCTGCACCCAGTCGATGCCCCGTGGTGCCTGGCTGCATTGCTTGAAGGCCAGGCATTTGCCATCCGCGTCAAGGCGGGCATAAGCGTT  >  NC_002947/870087‑870383

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