Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,492,130 A→T 7.1% L111Q (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,1300AT 7.1% 262.9 / 9.5 84L111Q (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 = 5.52e-01

CATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAAC  >  NC_002947/3491981‑3492276
                                                                                                                                                     |                                                                                                                                                   
cATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagca                                                                                                                                                     >  1:7968/1‑150 (MQ=255)
  tCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATATCGAGAAACGCATGGGCGGCATCGTCACGATACTGCACGTGGCGCTCGagcagc                                                                                                                                                   >  2:874513/1‑150 (MQ=255)
     aaGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGagcagcagc                                                                                                                                                >  1:1284754/1‑150 (MQ=255)
        cccGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCg                                                                                                                                             <  1:1388082/150‑1 (MQ=255)
         ccGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCgg                                                                                                                                            <  2:1284754/150‑1 (MQ=255)
         ccGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCgg                                                                                                                                            >  2:1437993/1‑150 (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:1615980/1‑150 (MQ=255)
                                                       ccAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGaa                                                                                              <  1:219624/150‑1 (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:1673747/150‑1 (MQ=255)
                                                         aGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAgg                                                                                            >  1:1674270/1‑150 (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                                                                                    >  2:2008400/1‑150 (MQ=255)
                                                                 tAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGc                                                                                    >  1:1919358/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                                                                   <  2:1615980/150‑1 (MQ=255)
                                                                                  gCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCaa                                                                   <  1:1103916/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:2220251/1‑150 (MQ=255)
                                                                                                    gagaAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCtt                                                 <  2:369242/150‑1 (MQ=255)
                                                                                                    gagaAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCtt                                                 >  1:26624/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                               >  2:2404804/1‑150 (MQ=255)
                                                                                                                     cATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGt                               >  1:1916750/1‑150 (MQ=255)
                                                                                                                      aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGcc                                                                                                    >  2:1131741/1‑81 (MQ=255)
                                                                                                                      aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGcc                                                                                                    <  1:1131741/81‑1 (MQ=255)
                                                                                                                      aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAg                                                                  >  1:984966/1‑115 (MQ=255)
                                                                                                                      aTCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAg                                                                  <  2:984966/115‑1 (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:874513/150‑1 (MQ=255)
                                                                                                                               aTACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCt                     <  1:1067750/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)
                                                                                                                                                     |                                                                                                                                                   
CATCGAAGCCCGAGGTGCCATTGCGCTCGTGCTCGATGTGGAATTTGTATTTTTGCCAGGCCTCGTAGAACAACGTGGTCATGCGGTGCGAGAAGATGTCGAGAAACGCATGGGCGGCATCGTCGCGATACTGCACGTGGCGCTCGAGCAGCAGCTCGGTGTACGGGCGTGGCAGCACGCCCGACGGCCCGACCAGGCCCATGAAGGTGACCTGCACTTCCGACAGCGGCAAGCCGGCCGGCGACACCTTGCCTTCGCGCTCGAAGTGCAGGTCGCTGACCTCGCTGGCCGGGAAC  >  NC_002947/3491981‑3492276

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: