Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 871,189 G→C 5.2% intergenic (+34/+41) PP_0753 → / ← PP_0754 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947871,1890GC 5.2% 311.4 / 2.8 97intergenic (+34/+41)PP_0753/PP_0754lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/46);  new base C (2/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTTCGGCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTAGACCAGCCACTG  >  NC_002947/871041‑871337
                                                                                                                                                    |                                                                                                                                                    
tttCGGCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGGTGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGt                                                                                                                                                     >  1:306471/1‑150 (MQ=255)
tttCGGCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGt                                                                                                                                                     <  2:203373/150‑1 (MQ=255)
  tCGGCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGa                                                                                                                                                   <  2:1854004/150‑1 (MQ=255)
    ggCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACa                                                                                                                                                 >  1:1903520/1‑150 (MQ=255)
    ggCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACa                                                                                                                                                 <  1:775976/150‑1 (MQ=255)
      cATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAgg                                                                                                                                               <  2:986020/150‑1 (MQ=255)
            gcCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCgg                                                                                                                                         >  1:298727/1‑150 (MQ=255)
                  ccAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGttttt                                                                                                                                   >  2:430401/1‑150 (MQ=255)
                    aGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGtt                                                                                                                                      <  2:306013/145‑1 (MQ=255)
                    aGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGtt                                                                                                                                      >  1:306013/1‑145 (MQ=255)
                     ggACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAAAAAACCCGCTGCCTGTCACTCAGGCAGCGGGttttcttg                                                                                                                                <  1:581274/150‑5 (MQ=38)
                      gACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAAAAAACCCGCTGCCTGTCACTCAGGCAGCGGGttttcttgc                                                                                                                               <  1:1883680/150‑6 (MQ=37)
                        cacCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTtgttgt                                                                                                                             <  1:667818/150‑1 (MQ=255)
                         acCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTc                                                                                                                            >  2:1292613/1‑150 (MQ=255)
                           cAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTg                                                                                                                          >  1:996659/1‑150 (MQ=255)
                              cacCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCCGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAg                                                                                                                                                      <  2:1009098/119‑1 (MQ=255)
                              cacCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCCGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAg                                                                                                                                                      >  1:1009098/1‑119 (MQ=255)
                                       ggtgggtgATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGAt                                                                                                              <  1:998618/150‑1 (MQ=255)
                                        gtgggtgATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGTTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTg                                                                                                                          <  2:1440579/137‑1 (MQ=255)
                                        gtgggtgATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTg                                                                                                                          >  1:1440579/1‑137 (MQ=255)
                                              gATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAcgcc                                                                                                       >  2:1059922/1‑150 (MQ=255)
                                               aTGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGa                                                                                                                                                   >  1:1757104/1‑105 (MQ=255)
                                               aTGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGa                                                                                                                                                   <  2:1757104/105‑1 (MQ=255)
                                                tgtgCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCgg                                                                                                                                         >  2:252674/1‑114 (MQ=255)
                                                tgtgCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCgg                                                                                                                                         <  1:252674/114‑1 (MQ=255)
                                                  tgCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCtt                                                                                                   <  2:1597404/150‑1 (MQ=255)
                                                     tGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTa                                                                                                >  2:1430272/1‑150 (MQ=255)
                                                      gAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTAt                                                                                               >  2:867011/1‑150 (MQ=255)
                                                        gCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTAtgt                                                                                             <  1:397028/150‑1 (MQ=255)
                                                         cGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTAtgtg                                                                                            <  1:430401/150‑1 (MQ=255)
                                                           aTGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGa                                                                                          >  2:379168/1‑150 (MQ=255)
                                                            tGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGaa                                                                                         >  2:883969/1‑150 (MQ=255)
                                                              cccTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCg                                                                                                                                          >  2:1533219/1‑99 (MQ=255)
                                                              cccTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCg                                                                                                                                          <  1:1533219/99‑1 (MQ=255)
                                                               ccTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTAACAGGCAGCGGGTTTTttgtt                                                                                                                               >  1:1940131/1‑109 (MQ=255)
                                                               ccTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTAACAGGCAGCGGGTTTTttgtt                                                                                                                               <  2:1940131/109‑1 (MQ=255)
                                                                cTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATcc                                                                                     <  2:349825/150‑1 (MQ=255)
                                                                 tGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    <  1:1923325/150‑1 (MQ=255)
                                                                 tGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    <  1:1484563/150‑1 (MQ=255)
                                                                 tGACCCGCCTGCTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    <  1:1212270/150‑1 (MQ=255)
                                                                    cccGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGt                                                                                                                                       >  1:45000/1‑96 (MQ=255)
                                                                    cccGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGt                                                                                                                                       <  2:45000/96‑1 (MQ=255)
                                                                     ccGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCaa                                                                                >  2:1221044/1‑150 (MQ=255)
                                                                     ccGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCaa                                                                                <  1:1336318/150‑1 (MQ=255)
                                                                               aaGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGATTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCaa                                                                                >  1:1623421/1‑140 (MQ=255)
                                                                               aaGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGATTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCaa                                                                                <  2:1623421/140‑1 (MQ=255)
                                                                               aaGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCgg                                                                                                                                         >  1:1397847/1‑83 (MQ=255)
                                                                               aaGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCgg                                                                                                                                         <  2:1397847/83‑1 (MQ=255)
                                                                                    cccGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGATTATATGAAAAAATCCAGAAATTTGAACCTAgatcc                                                                 >  2:1867699/1‑146 (MQ=255)
                                                                                     ccGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    >  2:591523/1‑130 (MQ=255)
                                                                                     ccGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    <  1:591523/130‑1 (MQ=255)
                                                                                     ccGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCg                                                                >  1:651998/1‑150 (MQ=255)
                                                                                       gAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGc                                                              >  2:1198997/1‑150 (MQ=255)
                                                                                           aTCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAcgc                                                                                                        <  1:905513/104‑1 (MQ=255)
                                                                                           aTCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAcgc                                                                                                        >  2:905513/1‑104 (MQ=255)
                                                                                            tCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCa                                                         >  2:1319732/1‑150 (MQ=255)
                                                                                                 cAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCg                                                    <  1:379168/150‑1 (MQ=255)
                                                                                                  aGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCgg                                                   <  1:1867699/150‑1 (MQ=255)
                                                                                                  aGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCgg                                                   <  1:1430272/150‑1 (MQ=255)
                                                                                                     tGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAc                                                                                                          >  1:547204/1‑92 (MQ=255)
                                                                                                     tGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAc                                                                                                          <  2:547204/92‑1 (MQ=255)
                                                                                                     tGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCa                                                >  1:370054/1‑150 (MQ=255)
                                                                                                      gCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCaa                                               <  2:298727/150‑1 (MQ=255)
                                                                                                          ggTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCa                                                                      >  1:2593/1‑123 (MQ=255)
                                                                                                          ggTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCa                                                                      <  2:2593/123‑1 (MQ=255)
                                                                                                             cAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGACATCCACTTCGGGCAAGTCCTCa                                        <  2:717527/150‑1 (MQ=255)
                                                                                                              aGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCAt                                       <  2:1361499/150‑1 (MQ=255)
                                                                                                              aGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGCCATCGCATCCACTTCGGGCAAGTCCTCAt                                       <  1:1718194/150‑1 (MQ=255)
                                                                                                                tAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCg                                                    <  2:1618063/135‑1 (MQ=255)
                                                                                                                tAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCg                                                    >  1:1618063/1‑135 (MQ=255)
                                                                                                                    cTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACAt                                 >  2:1620867/1‑150 (MQ=255)
                                                                                                                        aGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTcc                             <  2:1903520/150‑1 (MQ=255)
                                                                                                                             gtgCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAc                                                                                                          <  1:1168300/68‑1 (MQ=255)
                                                                                                                             gtgCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAc                                                                                                          >  2:1168300/1‑68 (MQ=255)
                                                                                                                             acaacaaaaaaCCCGCTGACTGTCACTCAGGCAGCGGGTTTTTCGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCCGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCTCATGTCCGCCaa                        >  2:718583/8‑150 (MQ=37)
                                                                                                                              tgCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAg                       >  2:770196/1‑150 (MQ=255)
                                                                                                                               gCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTc                                         <  1:1312875/131‑1 (MQ=255)
                                                                                                                               gCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTc                                         >  2:1312875/1‑131 (MQ=255)
                                                                                                                               gCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAgc                      <  1:1198997/150‑1 (MQ=255)
                                                                                                                               gCAAGAAAACCCGCTGCCTGACTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAgc                      <  2:996659/150‑1 (MQ=255)
                                                                                                                                cAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCg                                                    >  2:365860/1‑119 (MQ=255)
                                                                                                                                cAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCg                                                    <  1:365860/119‑1 (MQ=255)
                                                                                                                                cAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAgcg                     <  1:1059922/150‑1 (MQ=255)
                                                                                                                                cAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAgcg                     >  2:1041805/1‑150 (MQ=255)
                                                                                                                                 aaGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCa                                        <  1:876000/130‑1 (MQ=255)
                                                                                                                                 aaGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCa                                        >  2:876000/1‑130 (MQ=255)
                                                                                                                                 aaGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAgcgc                    <  1:1319732/150‑1 (MQ=255)
                                                                                                                                 aaGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGGCCTCATGCACATGTCCGCCAAgcgc                    >  1:1093701/1‑150 (MQ=255)
                                                                                                                                  aaaaaaCCCGCTGCCTGTCACTCAGGCAGCGAGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCa                                        >  2:972741/3‑129 (MQ=38)
                                                                                                                                  aGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    >  2:1031853/1‑85 (MQ=255)
                                                                                                                                  aGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCa                                                                                    <  1:1031853/85‑1 (MQ=255)
                                                                                                                                  aGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGca                   >  1:1359076/1‑150 (MQ=255)
                                                                                                                                       aCCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTa              >  1:1764301/1‑150 (MQ=255)
                                                                                                                                              gCCTGAGTGACAGGCAGCGGGTTTTTTGTGGGCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTAGACCAGc       >  1:1022168/1‑150 (MQ=255)
                                                                                                                                               ccTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTAGa            <  2:83845/144‑1 (MQ=255)
                                                                                                                                               ccTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTAGa            >  1:83845/1‑144 (MQ=255)
                                                                                                                                                  gAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAcg                                                                                                         >  1:992049/1‑48 (MQ=255)
                                                                                                                                                  gAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTAcg                                                                                                         <  2:992049/48‑1 (MQ=255)
                                                                                                                                                   aGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTAGACCAGCCACTg  <  1:867011/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCGGCATGGCGCCGAACCAGGACACCAACACCCAGCTGGTGGGTGATGTGCTGAGCGATGCCCTGACCCGCCTGTTCAAGGACCCGAGCATCGGGCAGGTGCTGGGTCAGTAAGCTTCAGGTGGTGCAAGAAAACCCGCTGCCTGAGTGACAGGCAGCGGGTTTTTTGTTGTCTGGCCGGTACAGATCTACGCCGCTTCTATGTGGAAATCCAGCAAGCTGAACCCAGTGCCGGCATCCACTTCGGGCAAGTCCTCATGCACATGTCCGCCAAGCGCACAGTAGACCAGCCACTG  >  NC_002947/871041‑871337

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: