Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,635,624 C→T 5.8% L331L (CTG→CTA) ‡ PP_2305 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,635,6240CT 5.8% 198.4 / 2.4 69L331L (CTG→CTA) ‡PP_2305hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

CCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTGGTTGATCAAATCCA  >  NC_002947/2635482‑2635772
                                                                                                                                               |                                                                                                                                                     
ccTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGa                                                                                                                                                <  1:1676421/150‑1 (MQ=255)
ccTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGa                                                                                                                                                >  2:522484/1‑150 (MQ=255)
  tCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATg                                                                                                                                              <  2:1692/150‑1 (MQ=255)
           cTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTCACGAAAAAGCAAAACTCAACTATCCGCCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAACTCCGgc                                                                                                                                     >  1:228184/1‑150 (MQ=255)
                tCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCtt                                                                                                                                >  1:1815570/1‑150 (MQ=255)
                tCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCCCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCACCCCAGGGCGATGAGCTCCGGCGCCtt                                                                                                                                <  1:724314/150‑1 (MQ=255)
                   ttttCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTa                                                                                                                             <  2:1095577/150‑1 (MQ=255)
                      tCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGaa                                                                                                                          <  1:1811303/150‑1 (MQ=255)
                         ccAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCTGGGGCGagagat                                                                                                                                           <  2:245474/130‑6 (MQ=255)
                         ccAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCTGGGGCGagagat                                                                                                                                           >  1:245474/1‑125 (MQ=255)
                          cAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGt                                                                                                                              >  2:1376660/1‑142 (MQ=255)
                          cAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGt                                                                                                                              <  1:1376660/142‑1 (MQ=255)
                             ccGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGaaa                                                                                                                  <  2:1815570/150‑1 (MQ=255)
                              cgcgAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGaaaa                                                                                                                 <  1:1865294/150‑1 (MQ=255)
                                        gTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCacac                                                                                                       >  1:246637/1‑150 (MQ=255)
                                               tCGAAGTTCACAGTGCCTTCAGGTGTGGCGCTAACCAGGTGCTAACGCTTAACCCCTTCATTACTATGCAGCGTATCGCGTAAAATCTCTCGCCGCAGAGCTATGAGTTCCGGCGCCTTATAGAACTAGAAAAACCGACACACCCGCATg                                                                                                >  2:127344/1‑150 (MQ=255)
                                               tCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTg                                                                                                >  1:1662163/1‑150 (MQ=255)
                                                cGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTgg                                                                                               <  1:1490459/150‑1 (MQ=255)
                                                  aaGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGt                                                                                             >  2:791463/1‑150 (MQ=255)
                                                      tCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCtt                                                                                                 <  1:796231/142‑1 (MQ=255)
                                                      tCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCtt                                                                                                 >  2:796231/1‑142 (MQ=255)
                                                          aGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGAt                                                                                     >  1:1154913/1‑150 (MQ=255)
                                                          aGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGAt                                                                                     <  2:378733/150‑1 (MQ=255)
                                                           gTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAAACCGGCACACCCGCTTGGGTACGTTGAt                                                                                     <  2:222694/150‑1 (MQ=255)
                                                           gTGCCTTAAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATc                                                                                    >  1:1155625/1‑150 (MQ=255)
                                                            tGCTTTCAGGGGTGGCGCAAACCAGGTGCCACGGCTTCACCCCTTCATAACTGTCCAGCGCATCGCGTCCAATCTCCCGCCCCAGCGAGATGAGCTCCGGCGCCGTGTAGAAATCAAAAAACCGGCACACCCGCTTGGGTACGTTGATCa                                                                                   <  2:246637/150‑1 (MQ=255)
                                                              ccTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGa                                                                                                                                             <  2:1375756/91‑1 (MQ=255)
                                                              ccTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGa                                                                                                                                             >  1:1375756/1‑91 (MQ=255)
                                                                          gcgtTAACCAGATGCTACCGCTTCACCCCTTCATCACTAACGAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGa                                                                     <  1:510075/150‑1 (MQ=255)
                                                                           cgcTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGAt                                                                    >  1:293303/1‑150 (MQ=255)
                                                                            gcTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATa                                                                   >  2:146678/1‑150 (MQ=255)
                                                                              tAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGc                                                                 <  2:556914/150‑1 (MQ=255)
                                                                                        cTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGaa                                                                                                                          <  1:273599/84‑1 (MQ=255)
                                                                                        cTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGaa                                                                                                                          <  1:272998/84‑1 (MQ=255)
                                                                                        cTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGaa                                                                                                                          >  2:272998/1‑84 (MQ=255)
                                                                                        cTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGaa                                                                                                                          >  2:273599/1‑84 (MQ=255)
                                                                                         tACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCATCGCCCTGGGGCGaga                                                                                                                                              >  2:985181/1‑62 (MQ=37)
                                                                                         tACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCATCGCCCTGGGGCGaga                                                                                                                                              <  1:985181/64‑3 (MQ=37)
                                                                                         tACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTg                                                      <  2:1662163/150‑1 (MQ=255)
                                                                                          aCCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCgg                                                                                                            >  2:350382/1‑95 (MQ=255)
                                                                                          aCCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCgg                                                                                                            <  1:350382/95‑1 (MQ=255)
                                                                                          aCCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGcc                                                                <  1:1541006/139‑1 (MQ=255)
                                                                                          aCCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGcc                                                                >  2:1541006/1‑139 (MQ=255)
                                                                                          aCCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCAGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGt                                                     >  1:668258/1‑150 (MQ=255)
                                                                                            cGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTa                                                    <  1:1794522/149‑1 (MQ=255)
                                                                                            cGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTa                                                    >  2:1794522/1‑149 (MQ=255)
                                                                                                 cACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGc                                              <  1:1382251/150‑1 (MQ=255)
                                                                                                    cccTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCaa                                           <  1:376418/150‑1 (MQ=255)
                                                                                                    cccTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCaa                                           >  1:517768/1‑150 (MQ=255)
                                                                                                      cTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCTAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACg                                         >  1:749469/1‑150 (MQ=255)
                                                                                                        tCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGaa                                       <  1:522484/150‑1 (MQ=255)
                                                                                                           tAACTGTCCAGCGTAGCGCGTGCAAGCTCTCCCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCAAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTc                                    <  2:228184/150‑1 (MQ=255)
                                                                                                            aaCTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATTAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCt                                   >  2:702491/1‑150 (MQ=255)
                                                                                                            aaCTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCtt                                                       >  2:1521919/1‑130 (MQ=255)
                                                                                                            aaCTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCtt                                                       <  1:1521919/130‑1 (MQ=255)
                                                                                                            aaCTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTACCAACGAAGTCt                                   <  2:1154913/150‑1 (MQ=255)
                                                                                                              cTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGc                                 >  2:888223/1‑150 (MQ=255)
                                                                                                               tGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCa                                <  1:1701824/150‑1 (MQ=255)
                                                                                                               tGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCa                                <  1:1701682/150‑1 (MQ=255)
                                                                                                                     gCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCAt                               <  1:67962/145‑1 (MQ=255)
                                                                                                                     gCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCAt                               >  2:67962/1‑145 (MQ=255)
                                                                                                                              gTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGctct                 <  1:888223/150‑1 (MQ=255)
                                                                                                                               tGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTg                >  2:865372/1‑150 (MQ=255)
                                                                                                                               tGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTg                >  2:223224/1‑150 (MQ=255)
                                                                                                                                gCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGt                                                                        >  2:1500786/1‑93 (MQ=255)
                                                                                                                                gCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGt                                                                        <  1:1500786/93‑1 (MQ=255)
                                                                                                                                  aaTCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGACGTCTGCATGACCTCGAGGCACTGAtt             >  2:898964/1‑150 (MQ=255)
                                                                                                                                    tctctcGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCTAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTGGTTGa           <  2:749469/150‑1 (MQ=255)
                                                                                                                                      tctcGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTGGTTGATc         <  2:1628447/150‑1 (MQ=255)
                                                                                                                                              ccAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGt                                     >  1:1593657/1‑115 (MQ=255)
                                                                                                                                              ccAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGt                                     <  2:1593657/115‑1 (MQ=255)
                                                                                                                                              ccAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTGGTTGATCAAATCca  <  1:865372/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CCTCGTGCGGCCTGTTTCGTTTTCGCCAGCCGCGAGCACGGTCATTGTCGAAGGTCACAGTGCCTTCAGGGGTGGCGCTAACCAGGTGCTACCGCTTCACCCCTTCATAACTGTCCAGCGTATCGCGTGCAATCTCTCGCCCCAGGGCGATGAGCTCCGGCGCCTTGTAGAACTCGAAAAACCGGCACACCCGCTTGGGTACGTTGATCAGCACATCCGGTGGATAGCCGGCAATCTTGTACTGTGCCAACGAAGTCTGCATGACCTCGAAGCTCTGGTTGATCAAATCCA  >  NC_002947/2635482‑2635772

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: