Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,866,051:1 +A 48.4% coding (871/1020 nt) PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,866,0511.A48.4% 100.5 / 144.6 62coding (871/1020 nt)PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (15/17);  new base A (13/17);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

AGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGC  >  NC_002947/3865904‑3866193
                                                                                                                                                    |                                                                                                                                              
aGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCatat                                                                                                                                              >  1:1234799/1‑150 (MQ=255)
      tctcGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGcc                                                                                                                                         >  1:1583039/1‑150 (MQ=255)
          gcgcACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCAAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCACCCCAACGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAgg                                                                                                                                     >  2:859639/1‑150 (MQ=255)
            gcACCATCAGCTCAAAGCCCAGTTCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                                                 <  1:1834976/150‑1 (MQ=255)
              aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCggg                                                                                                                                >  2:1420548/1‑150 (MQ=255)
                  tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCTCCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGAt                                                                                                                            >  2:1110890/1‑150 (MQ=255)
                  tCAGCTCAAAGCCCAGACCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATc                                                                                                                           >  1:483625/1‑150 (MQ=255)
                      cTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                                                        <  2:1234799/150‑1 (MQ=255)
                        cAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGGGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGa                                                                                                                      <  2:344547/150‑1 (MQ=255)
                         aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGa                                                                                                                      >  1:77563/1‑150 (MQ=255)
                         aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGaa                                                                                                                     >  1:1418297/1‑150 (MQ=255)
                             cccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAAcc                                                                                                                  <  2:1583039/150‑1 (MQ=255)
                              ccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCg                                                                                                                 >  1:660780/1‑150 (MQ=255)
                              ccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAg                                                                                                                                      >  2:851319/1‑128 (MQ=255)
                              ccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAg                                                                                                                                      <  1:851319/128‑1 (MQ=255)
                                aGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGAcc                                                                                                              <  2:1315382/150‑1 (MQ=255)
                                     ctgctgCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGc                                                                                                         <  2:1295459/150‑1 (MQ=255)
                                       gctgcGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCa                                                                                                       <  1:760909/150‑1 (MQ=255)
                                               ttCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGt                                                                                               <  2:1823778/150‑1 (MQ=255)
                                                  gTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGc                                                                                            >  2:1139811/1‑150 (MQ=255)
                                                      cccGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCt                                                                                        >  2:1554237/1‑150 (MQ=255)
                                                            tCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAccc                                                                                  <  2:1200467/150‑1 (MQ=255)
                                                               cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCtt                                                                               <  2:1418297/150‑1 (MQ=255)
                                                                   cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTcc                                                                                      <  1:344177/140‑1 (MQ=255)
                                                                   cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTcc                                                                                      >  2:344177/1‑140 (MQ=255)
                                                                     gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGc                                                                          >  1:1430792/1‑150 (MQ=255)
                                                                     gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGc                                                                          <  2:1691050/150‑1 (MQ=255)
                                                                      cAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCaa                                                                        >  2:1879175/1‑150 (MQ=255)
                                                                        aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                                                       <  2:77563/150‑1 (MQ=255)
                                                                        aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                                                       <  1:1429136/150‑1 (MQ=255)
                                                                        aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAgg                                                                      >  1:1838117/1‑150 (MQ=255)
                                                                          gCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                                                        >  1:1228697/1‑98 (MQ=255)
                                                                          gCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                                                        <  2:1228697/98‑1 (MQ=255)
                                                                              gcagcGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCAGCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgca                                                                <  2:177534/150‑1 (MQ=255)
                                                                               cagcGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGAGAGCATCAAGGTCAGCCGCGCCACCAAAACACGCGGTTTGCGTGTGATCGTGTAAGCCTACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgca                                                                <  2:218599/150‑1 (MQ=255)
                                                                                   gcCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCAc                                                            <  1:1033717/150‑1 (MQ=255)
                                                                                   gcCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCAc                                                            <  2:660780/150‑1 (MQ=255)
                                                                                     cTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACcg                                                          >  2:1009660/1‑150 (MQ=255)
                                                                                        tgCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGcc                                                      <  1:1412657/150‑1 (MQ=255)
                                                                                               gcccgccccACCGCAGCCAGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGAGCGTTGAAACCGACCATCGATACCTGGTGTGGCACCTACACCCCTTGATGAAAGGCGTGTCGCACCTCGCCCTGAgac                                               >  2:827418/2‑148 (MQ=255)
                                                                                                cccgccccACCGCAGACCTCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCa                                              >  1:66670/1‑150 (MQ=255)
                                                                                                  cgccccTCCGCAGCCCGCGGTCGGCGAATGGAGGTCAGCCGTGCCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCGCCGCCCCGTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAg                                             <  2:625408/150‑1 (MQ=255)
                                                                                                  cgccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAg                                             <  1:698852/150‑1 (MQ=255)
                                                                                                  caccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCGCGTGGCTTAAGGCCTGCAGCACCGCGCCCTGAGCCAg                                             <  1:859639/148‑1 (MQ=255)
                                                                                                    ccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGcc                                                                                                                                         <  2:11932/56‑1 (MQ=255)
                                                                                                    ccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGcc                                                                                                                                         >  1:11932/1‑56 (MQ=255)
                                                                                                        aCCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGtcgtcg                                      >  2:986604/1‑150 (MQ=255)
                                                                                                         ccGCAGCCCGCTGTGGGCGAAGGGGTGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAGCCCGCCCTCGCTACCTGGCGCGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTcgtcgt                                     <  2:66670/150‑1 (MQ=255)
                                                                                                           gCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTGGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTa                                   <  1:1554237/150‑1 (MQ=255)
                                                                                                            cAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTa                                   <  1:1617337/150‑1 (MQ=255)
                                                                                                                    ggTGTGCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACagaagaa                           <  2:77258/150‑1 (MQ=255)
                                                                                                                      tgtgCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGt                                                                                               >  1:1374355/1‑80 (MQ=255)
                                                                                                                      tgtgCGAATGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGt                                                                                               <  2:1374355/80‑1 (MQ=255)
                                                                                                                        tgCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATAc                      >  2:384204/1‑150 (MQ=255)
                                                                                                                            aaTGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGt                   >  1:1514572/1‑150 (MQ=255)
                                                                                                                            aaTGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTc                  <  1:1139811/150‑1 (MQ=255)
                                                                                                                             aTGGAGGTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTc                  >  2:1426221/1‑150 (MQ=255)
                                                                                                                                  ggTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCt                                                    <  2:1417376/111‑1 (MQ=255)
                                                                                                                                  ggTCAGCCGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCt                                                    >  1:1417376/1‑111 (MQ=255)
                                                                                                                                     cAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCgg         >  1:55858/1‑150 (MQ=255)
                                                                                                                                        ccGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTg                                                                                              >  1:189868/1‑63 (MQ=255)
                                                                                                                                        ccGTGGCACCATAATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTg                                                                                              <  2:189868/63‑1 (MQ=255)
                                                                                                                                            ggCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCTTCGACATCCGGCGCctgc  <  2:1838117/150‑1 (MQ=255)
                                                                                                                                                  atatGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgc                                                                 >  2:762235/1‑81 (MQ=255)
                                                                                                                                                  atatGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgc                                                                 <  1:762235/81‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAGAAGAAGATACCGTCGACATCCGGCGCCTGC  >  NC_002947/3865904‑3866193

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