Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,492,133 A→T 7.1% L110Q (CTG→CAG)  PP_3096 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,492,1330AT 7.1% 263.2 / 8.8 84L110Q (CTG→CAG) PP_3096hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base T (4/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

AAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCA  >  NC_002947/3491986‑3492280
                                                                                                                                                   |                                                                                                                                                    
aaGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagcagcagc                                                                                                                                                    >  1:1284754/1‑150 (MQ=255)
   cccGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCg                                                                                                                                                 <  1:1388082/150‑1 (MQ=255)
    ccGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCgg                                                                                                                                                >  2:1437993/1‑150 (MQ=255)
    ccGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCgg                                                                                                                                                <  2:1284754/150‑1 (MQ=255)
     cGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGgt                                                                                                                                               <  2:1509020/150‑1 (MQ=255)
      gAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCTGCTGCTCgagc                                                                                                                                              <  1:618682/150‑4 (MQ=255)
             cATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGc                                                                                                                                       >  2:1403893/1‑150 (MQ=255)
               ttGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGt                                                                                                                                     <  2:2363428/150‑1 (MQ=255)
                 gcgcTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTgg                                                                                                                                   <  1:2233601/150‑1 (MQ=255)
                         gCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGc                                                                                                                           <  2:2235132/150‑1 (MQ=255)
                          cTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAACCCGcc                                                                                                                          >  2:894639/1‑150 (MQ=255)
                             gATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGa                                                                                                                       <  2:2101306/150‑1 (MQ=255)
                             gATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGa                                                                                                                       >  2:1103916/1‑150 (MQ=255)
                              aTGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGAc                                                                                                                      <  1:756930/150‑1 (MQ=255)
                                            tttttGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGcc                                                                                                        <  1:2390508/150‑1 (MQ=255)
                                                 gCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGa                                                                                                   >  1:1939830/1‑150 (MQ=255)
                                                  tcgggCCTCTTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGaa                                                                                                  >  2:2211610/4‑150 (MQ=255)
                                                  ccAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGaa                                                                                                  <  1:219624/150‑1 (MQ=255)
                                                  ccAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGaa                                                                                                  >  1:1615980/1‑150 (MQ=255)
                                                   cAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTAGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTTTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAg                                                                                                 >  2:272326/1‑150 (MQ=255)
                                                   cAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCACCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAg                                                                                                 >  1:301412/1‑150 (MQ=255)
                                                   cAGGCCTCATAGAACAACGTGGTCATGCGGTGCGAGAAGATGTAGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCACCAGCTCGGTGTACGGGCGTCGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAg                                                                                                 <  2:407525/150‑1 (MQ=255)
                                                   cAGGCCTACTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAg                                                                                                 <  1:1309568/150‑1 (MQ=255)
                                                    aGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAgg                                                                                                >  1:1674270/1‑150 (MQ=255)
                                                    aGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAgg                                                                                                <  1:1673747/150‑1 (MQ=255)
                                                        cTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                            <  2:1363903/150‑1 (MQ=255)
                                                         tCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGgtg                                                                                                                                              <  2:811620/99‑1 (MQ=255)
                                                         tCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGgtg                                                                                                                                              >  1:811620/1‑99 (MQ=255)
                                                         tCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAcc                                                                                           >  1:971611/1‑150 (MQ=255)
                                                            tAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGc                                                                                        >  1:1919358/1‑150 (MQ=255)
                                                            tAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGc                                                                                        >  2:2008400/1‑150 (MQ=255)
                                                             aGAACATCGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCa                                                                                       <  2:2406894/150‑1 (MQ=255)
                                                               aacaacGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGccc                                                                                                                 <  2:714581/122‑1 (MQ=255)
                                                               aacaacGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGccc                                                                                                                 >  1:714581/1‑122 (MQ=255)
                                                               aacaacGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACt                                                                                     >  1:198579/1‑150 (MQ=255)
                                                                  aacGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGgtg                                                                                                                                              >  1:967246/1‑90 (MQ=255)
                                                                  aacGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGgtg                                                                                                                                              <  2:967246/90‑1 (MQ=255)
                                                                      tGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACa                                                                              >  2:1067750/1‑150 (MQ=255)
                                                                             gCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCaa                                                                       <  1:1103916/150‑1 (MQ=255)
                                                                             gCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCaa                                                                       <  2:1615980/150‑1 (MQ=255)
                                                                                 tGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAgccg                                                                   >  1:2388283/1‑150 (MQ=255)
                                                                                    gagaAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACCTCCGACAGCGGCAAgccggcc                                                                >  1:607983/1‑150 (MQ=255)
                                                                                    gagaAGATGTCGAGAAACGCATGGGCGGCATCGTCACGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAgcccgcc                                                                >  2:361803/1‑150 (MQ=255)
                                                                                             tCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACAcc                                                       <  1:1437993/150‑1 (MQ=255)
                                                                                               gagaAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCtt                                                     <  2:369242/150‑1 (MQ=255)
                                                                                               gagaAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCtt                                                     >  1:26624/1‑150 (MQ=255)
                                                                                               gagaAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCtt                                                     >  2:2220251/1‑150 (MQ=255)
                                                                                                 gaATCGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCTACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGc                                                   <  2:260971/150‑1 (MQ=255)
                                                                                                 gaAAGGCATGGGCGCCATCGTCGAGATACTGTACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGANGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGc                                                   <  1:941562/150‑1 (MQ=255)
                                                                                                  acaCGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGcc                                                 <  2:1950734/148‑1 (MQ=255)
                                                                                                  aaaCGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGcc                                                 >  2:493152/1‑150 (MQ=255)
                                                                                                      gCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTcg                                             >  1:1736657/1‑150 (MQ=255)
                                                                                                       cATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTcgc                                            <  1:2220251/150‑1 (MQ=255)
                                                                                                             cggcATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGt                                                                                               >  2:121384/1‑94 (MQ=255)
                                                                                                             cggcATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGt                                                                                               <  1:121384/94‑1 (MQ=255)
                                                                                                             cggcATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGa                                      >  1:403310/1‑150 (MQ=255)
                                                                                                              ggcATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGaa                                     >  2:1464135/1‑150 (MQ=255)
                                                                                                              ggcATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGaa                                     <  1:1403893/150‑1 (MQ=255)
                                                                                                                cATCGTCGCGATACTTCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGACGGCGACACCTTGCCTTCGCGCTCGAAGt                                   <  2:1674270/150‑1 (MQ=255)
                                                                                                                cATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGt                                   >  1:1916750/1‑150 (MQ=255)
                                                                                                                cATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGt                                   >  2:2404804/1‑150 (MQ=255)
                                                                                                                 aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGcc                                                                                                        >  2:1131741/1‑81 (MQ=255)
                                                                                                                 aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGcc                                                                                                        <  1:1131741/81‑1 (MQ=255)
                                                                                                                 aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAg                                                                      <  2:984966/115‑1 (MQ=255)
                                                                                                                 aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAg                                                                      >  1:984966/1‑115 (MQ=255)
                                                                                                                    gtcgCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCACTTCGAGCGCGAAggcaa                                <  2:971611/150‑5 (MQ=255)
                                                                                                                       gcgATACTTCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTc                            >  1:412703/1‑150 (MQ=255)
                                                                                                                        cgATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCg                           <  2:1919358/150‑1 (MQ=255)
                                                                                                                          aTACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCt                         <  1:1067750/150‑1 (MQ=255)
                                                                                                                          aTACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCt                         <  1:874513/150‑1 (MQ=255)
                                                                                                                           tACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAGTGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTg                        <  1:164397/150‑1 (MQ=255)
                                                                                                                             cTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                            <  2:1706663/81‑1 (MQ=255)
                                                                                                                             cTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                            >  1:1706663/1‑81 (MQ=255)
                                                                                                                              ggCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCGGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGAcc                     <  2:301412/149‑1 (MQ=255)
                                                                                                                                cACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                            <  1:1099574/78‑1 (MQ=255)
                                                                                                                                cACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGAc                                                                                            >  2:1099574/1‑78 (MQ=255)
                                                                                                                                   gTGGAGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCt                >  1:2153164/1‑150 (MQ=255)
                                                                                                                                   gTGGAGCCCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCt                <  2:1978106/150‑1 (MQ=255)
                                                                                                                                        gcTCGAGCAGCAGCTCGGTGTACGGGAGCGGGCGCCCGATGCCCGTGCTTCGGAGACCCATGAAGGTGACTAGCACCTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGGAGTGCAGGTCGCTGACCTCGCTGGCCg           >  1:421207/1‑150 (MQ=255)
                                                                                                                                         accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGc                                                   <  2:1975377/108‑1 (MQ=255)
                                                                                                                                         accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCgg          >  1:500040/3‑150 (MQ=255)
                                                                                                                                         accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCgg          >  2:1495556/3‑150 (MQ=255)
                                                                                                                                         accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCgg          >  2:1622769/3‑150 (MQ=255)
                                                                                                                                         accGAGCTGCTGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCGGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGc                                                   >  1:1975377/3‑110 (MQ=255)
                                                                                                                                          tCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCggg         >  1:1129874/1‑150 (MQ=255)
                                                                                                                                             agcagcagcTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAAc      <  2:2113487/150‑1 (MQ=255)
                                                                                                                                                 gcagcTCGGTGTACGAGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGACa  <  1:493152/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAACGCCA  >  NC_002947/3491986‑3492280

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