Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,612 C→T 5.6% V541V (GTC→GTT)  PP_3340 → TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,778,6120CT 5.6% 262.4 / 4.0 89V541V (GTC→GTT) PP_3340TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base C (47/37);  new base T (2/3);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GTCAAGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCATCAGC  >  NC_002947/3778463‑3778759
                                                                                                                                                     |                                                                                                                                                   
gTCAAGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTc                                                                                                                                                     >  2:193458/1‑150 (MQ=255)
 tCAAGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCa                                                                                                                                                    >  2:52080/1‑150 (MQ=255)
   aaGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAAc                                                                                                                                                  <  1:289775/150‑1 (MQ=255)
    aGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAAcc                                                                                                                                                 <  2:1442771/150‑1 (MQ=255)
     gTTCACCCCGGCCAGCTGGGTCGATGGCCCCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCt                                                                                                                                                <  1:1436045/150‑1 (MQ=255)
      ttCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTTGACCTg                                                                                                                                               >  2:1271370/1‑150 (MQ=255)
      ttCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGACCGAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTg                                                                                                                                               >  1:1534851/1‑150 (MQ=255)
           cccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTTGACCTGCaggt                                                                                                                                          <  1:1271370/150‑4 (MQ=255)
           cccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCAccc                                                                                                                                          >  1:1824483/1‑150 (MQ=255)
              ggCCAGCTGGGTAGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTTGACCTGCaggtcat                                                                                                                                       <  1:1832147/150‑7 (MQ=255)
               gCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACATCGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGac                                                                                                                                      <  2:1206670/150‑1 (MQ=255)
               gCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGac                                                                                                                                      >  2:63672/1‑150 (MQ=255)
                ccAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGacg                                                                                                                                     >  2:1723559/1‑150 (MQ=255)
                   gCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGacgaac                                                                                                                                  <  1:295054/150‑1 (MQ=255)
                   gCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGacgaac                                                                                                                                  <  2:459043/150‑1 (MQ=255)
                       ggTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACgaacgaac                                                                                                                              <  1:193458/150‑1 (MQ=255)
                           gATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGa                                                                                                                           <  1:316330/149‑1 (MQ=255)
                           gATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGa                                                                                                                           >  2:316330/1‑149 (MQ=255)
                             tGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGgt                                                                                                                        <  1:52080/150‑1 (MQ=255)
                             tGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGgt                                                                                                                        <  2:1743528/150‑1 (MQ=255)
                             tGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGgt                                                                                                                        >  2:278968/1‑150 (MQ=255)
                                  gcgTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGa                                                                                                                   >  2:1317567/1‑150 (MQ=255)
                                   cgTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGAt                                                                                                                  <  1:167345/150‑1 (MQ=255)
                                       gCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGacgaacg                                                                                                                                 <  1:1653970/131‑1 (MQ=255)
                                       gCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCAATGACCTGCAGGTCAACCTGCACCCGGacgaacg                                                                                                                                 >  2:1653970/1‑131 (MQ=255)
                                          tATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTAc                                                                                                           >  2:1266216/1‑150 (MQ=255)
                                          tATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTAc                                                                                                           >  1:87104/1‑150 (MQ=255)
                                          tATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTAc                                                                                                           >  1:970171/1‑150 (MQ=255)
                                               gCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTg                                                                                                      >  1:1249105/1‑150 (MQ=255)
                                                 gTCAGGACGCGTCTGGCGAAGGGAGCCGTCGCCTGAACGACCCCAGAGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGTTATGACCTGCATGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGc                                                                                                    <  2:361835/150‑1 (MQ=255)
                                                 gTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGCGCCGGACAACGTCGGCGAAACACGCCGCTGTGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGc                                                                                                    <  2:258378/150‑1 (MQ=255)
                                                 gTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGc                                                                                                    <  2:264061/150‑1 (MQ=255)
                                                      gACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCattgc                                                                                               >  1:790324/1‑146 (MQ=255)
                                                           ctctGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCaa                                                                                          >  1:438942/1‑150 (MQ=255)
                                                            tctGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAg                                                                                         >  1:1740703/1‑150 (MQ=255)
                                                                  gAAGTGAGCCGTCGCCTGAATGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCtt                                                                                   <  2:1610768/150‑1 (MQ=255)
                                                                    aGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGa                                                                                 <  2:894251/150‑1 (MQ=255)
                                                                         gCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAg                                                                                                                          >  2:1355729/1‑104 (MQ=255)
                                                                         gCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAg                                                                                                                          <  1:1355729/104‑1 (MQ=255)
                                                                                ccTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACAGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAgcgc                                                                     >  1:807856/1‑150 (MQ=255)
                                                                                   gAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  <  2:970171/150‑1 (MQ=255)
                                                                                   caaCGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  <  1:1848804/149‑1 (MQ=255)
                                                                                    aaCGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  >  1:869228/1‑149 (MQ=255)
                                                                                    aaCGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  <  2:869228/149‑1 (MQ=255)
                                                                                     aCGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCAAAGCATCGGGTTga                                                                >  2:1848804/1‑149 (MQ=255)
                                                                                       gACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCa                                                              <  2:1612769/150‑1 (MQ=255)
                                                                                       gACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCa                                                              <  2:438942/150‑1 (MQ=255)
                                                                                         ccccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                            <  1:1317567/150‑1 (MQ=255)
                                                                                         ccccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                            >  1:874370/1‑150 (MQ=255)
                                                                                          cccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGgc                                                           >  2:1597988/1‑150 (MQ=255)
                                                                                          cccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCAAGgc                                                           >  1:770706/1‑150 (MQ=255)
                                                                                           ccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGgca                                                          <  1:428809/150‑1 (MQ=255)
                                                                                            cAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGgcag                                                         >  1:1538408/1‑150 (MQ=255)
                                                                                            cAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGgcag                                                         >  1:1397884/1‑150 (MQ=255)
                                                                                             aGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGcc                                                                              >  1:887896/1‑128 (MQ=255)
                                                                                             aGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGcc                                                                              <  2:887896/128‑1 (MQ=255)
                                                                                              gcggcgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGAcgaacgaa                                                                                                                               <  1:1148817/78‑1 (MQ=255)
                                                                                              gcggcgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGAcgaacgaa                                                                                                                               >  2:1148817/1‑78 (MQ=255)
                                                                                                ggcgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCaaa                                                     <  2:1249105/150‑1 (MQ=255)
                                                                                                   gAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  >  1:1350834/1‑134 (MQ=255)
                                                                                                   gAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  <  2:1350834/134‑1 (MQ=255)
                                                                                                    aGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTgg                                                                                                                    <  1:806677/83‑1 (MQ=255)
                                                                                                    aGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTgg                                                                                                                    >  2:806677/1‑83 (MQ=255)
                                                                                                    aGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGcc                                                 >  2:601208/1‑150 (MQ=255)
                                                                                                      tCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAg                                               <  1:601208/150‑1 (MQ=255)
                                                                                                           cAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAAt                                          <  1:1141994/150‑1 (MQ=255)
                                                                                                                    cGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGt                                 <  1:278968/150‑1 (MQ=255)
                                                                                                                         cacGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCcac                            >  1:16398/1‑150 (MQ=255)
                                                                                                                          acGCCGCTGGGGCTATGACCTGCGGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCGGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCTAGGAAATCGATCACGTGAcaca                           >  2:949666/1‑150 (MQ=255)
                                                                                                                          acGCCGCTGGGGCTATGACCTGCGGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCGGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCGAGGAAAACGATCACGTGGcacc                           >  2:950012/1‑149 (MQ=255)
                                                                                                                          acGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCAc                                   >  2:447697/1‑142 (MQ=255)
                                                                                                                          acGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCAc                                   <  1:447697/142‑1 (MQ=255)
                                                                                                                          acGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCcaca                           >  2:1420632/1‑150 (MQ=255)
                                                                                                                          acGCCGCTGGGGCTATGACCTACAGGTCAACCAGAACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCcaca                           >  1:912172/1‑150 (MQ=255)
                                                                                                                              cgcTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCAt                       >  2:1724306/1‑150 (MQ=255)
                                                                                                                                  ggggCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCtgt                   <  2:87104/150‑1 (MQ=255)
                                                                                                                                     gCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCTAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGa                >  2:78173/1‑150 (MQ=255)
                                                                                                                                     gCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGa                >  1:1566325/1‑150 (MQ=255)
                                                                                                                                       tATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAac              >  2:201040/1‑150 (MQ=255)
                                                                                                                                          cgggtGCAGGTTGACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGa                                       <  2:510676/118‑1 (MQ=255)
                                                                                                                                          cgggtGCAGGTTGACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGa                                       >  1:510676/5‑122 (MQ=255)
                                                                                                                                             cTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGc        >  1:1371677/1‑150 (MQ=255)
                                                                                                                                             cTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGc        >  2:155841/1‑150 (MQ=255)
                                                                                                                                              tGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCa       >  1:1723519/1‑150 (MQ=255)
                                                                                                                                                cAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCATc     >  2:569501/1‑150 (MQ=255)
                                                                                                                                                 aGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGaaa                                           >  1:1558877/1‑111 (MQ=255)
                                                                                                                                                 aGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGaaa                                           <  2:1558877/111‑1 (MQ=255)
                                                                                                                                                   gTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCATCAGc  >  2:1750918/1‑150 (MQ=255)
                                                                                                                                                    tCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAgg                                                   >  2:646227/1‑100 (MQ=255)
                                                                                                                                                    tCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAgg                                                   <  1:646227/100‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTCAAGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCATCAGC  >  NC_002947/3778463‑3778759

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: ‑