Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,613 A→G 5.7% N542D (AAC→GAC)  PP_3340 → TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,778,6130AG 5.7% 269.1 / 4.1 88N542D (AAC→GAC) PP_3340TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (46/37);  new base G (2/3);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

TCAAGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCATCAGC  >  NC_002947/3778464‑3778759
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                           >  2:316330/1‑149 (MQ=255)
                          gATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGa                                                                                                                           <  1:316330/149‑1 (MQ=255)
                            tGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGgt                                                                                                                        <  1:52080/150‑1 (MQ=255)
                            tGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGgt                                                                                                                        >  2:278968/1‑150 (MQ=255)
                            tGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGgt                                                                                                                        <  2:1743528/150‑1 (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                                                                                                           >  1:970171/1‑150 (MQ=255)
                                         tATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTAc                                                                                                           >  1:87104/1‑150 (MQ=255)
                                         tATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTAc                                                                                                           >  2:1266216/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:438942/150‑1 (MQ=255)
                                                                                      gACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCa                                                              <  2:1612769/150‑1 (MQ=255)
                                                                                        ccccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                            >  1:874370/1‑150 (MQ=255)
                                                                                        ccccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                            <  1:1317567/150‑1 (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                                                                                                                               >  2:1148817/1‑78 (MQ=255)
                                                                                             gcggcgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGAcgaacgaa                                                                                                                               <  1:1148817/78‑1 (MQ=255)
                                                                                               ggcgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCaaa                                                     <  2:1249105/150‑1 (MQ=255)
                                                                                                  gAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  <  2:1350834/134‑1 (MQ=255)
                                                                                                  gAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGtt                                                                  >  1:1350834/1‑134 (MQ=255)
                                                                                                   aGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTgg                                                                                                                    >  2:806677/1‑83 (MQ=255)
                                                                                                   aGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTgg                                                                                                                    <  1:806677/83‑1 (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                                       >  1:510676/5‑122 (MQ=255)
                                                                                                                                         cgggtGCAGGTTGACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGa                                       <  2:510676/118‑1 (MQ=255)
                                                                                                                                            cTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGc        >  2:155841/1‑150 (MQ=255)
                                                                                                                                            cTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGc        >  1:1371677/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)
                                                                                                                                                     |                                                                                                                                                  
TCAAGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGGCATCAGC  >  NC_002947/3778464‑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: ‑