Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,973,159 A→G 5.6% intergenic (+153/+32) PP_3502 → / ← PP_3503 transposase/sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,973,1590AG 5.6% 266.1 / 2.4 88intergenic (+153/+32)PP_3502/PP_3503transposase/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base G (3/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTGCAACGCAGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGCACATCGCACAGTGCCTGG  >  NC_002947/3973015‑3973308
                                                                                                                                                |                                                                                                                                                     
cTGCAACGCAGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGc                                                                                                                                                  <  2:1519434/150‑1 (MQ=32)
 tacAACGCAGCGCGCCCAAACTGACGGGCAGTATCCAAAGCACGGACATCACTGGCGAGCACGCTGACTCGTACCTCGAGACTGCAGAGAGCCCCCGATGGCGTTTTGAATGAAAAGTTGGCTGCGGCTATAGAATCTCCTACAGGGGcg                                                                                                                                                 >  2:339667/3‑149 (MQ=2)
    aaCGCAGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTtgt                                                                                                                                              <  1:1218384/150‑1 (MQ=35)
    aaCGCAGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTtgt                                                                                                                                              <  1:2392916/150‑1 (MQ=35)
     aCGCAGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTtgtg                                                                                                                                             <  1:141557/150‑1 (MQ=35)
         aGCCCGCCCAAACTGACGGGCAGTATCCAATTCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGttt                                                                                                                                         <  1:1578774/150‑1 (MQ=255)
         aGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGttt                                                                                                                                         <  1:2398912/150‑1 (MQ=39)
             cgcccAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCt                                                                                                                                     <  1:1199397/150‑1 (MQ=255)
              gcccAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCtt                                                                                                                                    <  1:1875997/150‑1 (MQ=255)
              gcccAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCttt                                                                                                                                   <  2:1321003/150‑1 (MQ=255)
               cccAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCttt                                                                                                                                   <  2:2067730/150‑1 (MQ=255)
                 cAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCtt                                                                                                                                    <  2:977727/147‑1 (MQ=255)
                 cAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCtt                                                                                                                                    >  1:977727/1‑147 (MQ=255)
                  aaaCTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTgag                                                                                                                                >  2:2253782/1‑150 (MQ=255)
                        aCGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCNTACCTCGAGGTTGCGGAGAGAGGCCGATGGCGTTTTTGCTGAAAATTTACCTCCGTCCAGAGGAGCTCCTACAGGGGCTCGAGGTTTAGCTTTGagagagat                                                                                                                          >  1:976748/1‑146 (MQ=11)
                        aCGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGcgcc                                                                                                                          <  2:2233031/150‑1 (MQ=255)
                                  tCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTa                                                                                                                <  1:2262172/150‑1 (MQ=255)
                                   ccAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTAt                                                                                                               >  1:546848/1‑150 (MQ=255)
                                     aaTGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCt                                                                                                             >  2:696880/1‑150 (MQ=255)
                                            gggAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCtt                                                                                                      <  2:2141745/150‑1 (MQ=255)
                                              gagaAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGCTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTtgt                                                                                                    >  1:1737116/1‑150 (MQ=255)
                                               agaAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTtgtg                                                                                                   <  2:1739483/150‑1 (MQ=255)
                                                  aaTTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGcaa                                                                                                <  2:387999/150‑1 (MQ=255)
                                                    ttCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGcaaca                                                                                              >  1:1913272/1‑150 (MQ=255)
                                                      ctctGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTCACTTCTGCGCAGGGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGcaacatc                                                                                            <  2:1771523/150‑3 (MQ=21)
                                                        ctGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCg                                                                                          >  2:164250/1‑150 (MQ=255)
                                                           gCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATa                                                                                       >  1:1199073/1‑150 (MQ=255)
                                                             aCGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTga                                                                                                                                 >  2:1762506/1‑106 (MQ=255)
                                                             aCGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTga                                                                                                                                 <  1:1762506/106‑1 (MQ=255)
                                                                 tGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGt                                                                                 >  1:585116/1‑150 (MQ=255)
                                                                        tACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGt                                                                          >  2:1070774/1‑150 (MQ=255)
                                                                           cTCGAGGTTGCTTAGAGAGGGCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGACTCTCCTACAGGGACTTCTGCCTTCCCTTTGAGAGCGCCGCTCAGCGTGTCTGGTTCTGGTGCAACAACCGATAGCTCGTCGGCCGTGaa                                                                       >  2:66733/1‑150 (MQ=21)
                                                                           cTCGAGGTTGCTGAGAGAGGACGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACGGGGGCatccaggctggtcaatctcttctgtaacggcaaagcggcgacttggcaccccaccccttcagggttatgccgagc                                                                       >  2:1894140‑M2/1‑75 (MQ=255)
                                                                            tCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGaaa                                                                      <  1:164250/150‑1 (MQ=255)
                                                                            tCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGaaa                                                                      >  1:1583283/1‑150 (MQ=255)
                                                                             cGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAAt                                                                     <  2:1539895/150‑1 (MQ=255)
                                                                                ggTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGcc                                                                  <  1:696880/150‑1 (MQ=255)
                                                                                  ttGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCa                                                                <  2:1583283/150‑1 (MQ=255)
                                                                                      tGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACGGGGGCatccaggctggtcaatctgttctgtaacgacaaagcggcgacttggcaccccaccccttcagggttatgccgagc                                                                       <  1:885507‑M2/139‑76 (MQ=255)
                                                                                      tGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACGGGGGCatccaggctggtcaatctgttctgtaacgacaaagcggcgacttggcaccccaccccttcagggttatgccgagc                                                                       >  2:885507‑M2/1‑64 (MQ=255)
                                                                                      tGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGc                                                                                                                   <  1:614027/95‑1 (MQ=255)
                                                                                      tGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGc                                                                                                                   >  2:614027/1‑95 (MQ=255)
                                                                                            agGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGTCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGc                                                      >  2:1216049/1‑150 (MQ=255)
                                                                                            agGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGc                                                                                                                   <  1:2066513/89‑1 (MQ=255)
                                                                                            agGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGc                                                                                                                   >  2:2066513/1‑89 (MQ=255)
                                                                                               ccGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGc                                                   <  2:597078/150‑1 (MQ=255)
                                                                                                       gTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCa                                           >  2:1460658/1‑150 (MQ=255)
                                                                                                        tttttACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTCGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAg                                          <  1:2346508/150‑1 (MQ=255)
                                                                                                              cTGAAAATTTGCCTCCGGCCATAGAATCTCCTACGGGGGCatccaggctggtcaatctgttctgtaacgacaaagcggcgacttggcaccccaccccttcagggttatgccgagcc                                                                      <  2:111913‑M2/116‑77 (MQ=255)
                                                                                                              cTGAAAATTTGCCTCCGGCCATAGAATCTCCTACGGGGGCatccaggctggtcaatctgttctgtaacgacaaagcggcgacttggcaccccaccccttcagggttatgccgagcc                                                                      >  1:111913‑M2/1‑40 (MQ=255)
                                                                                                              cTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATGTGGGGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTc                                    >  2:1771631/1‑150 (MQ=255)
                                                                                                              cTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCGGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTc                                    <  1:2299117/150‑1 (MQ=255)
                                                                                                              cTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTc                                    <  2:1689173/150‑1 (MQ=255)
                                                                                                               tGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCg                                   >  1:1343186/1‑150 (MQ=255)
                                                                                                               tGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCg                                   >  2:679279/1‑150 (MQ=255)
                                                                                                                gAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCg                                                                               <  2:904418/105‑1 (MQ=255)
                                                                                                                gAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCg                                                                               >  1:904418/1‑105 (MQ=255)
                                                                                                                gAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCGCCTTGGCTGCCACGCTCAGGTTGTCGc                                  <  2:546848/150‑1 (MQ=255)
                                                                                                                 aaatgttGCCCCCGGCCATAGAAGCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCt                                 <  1:339667/145‑1 (MQ=255)
                                                                                                                    aTTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCt                                 <  1:1253411/147‑1 (MQ=255)
                                                                                                                    aTTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCt                                 >  2:1253411/1‑147 (MQ=255)
                                                                                                                     tttGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGc                             >  1:1982689/1‑150 (MQ=255)
                                                                                                                         ccTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAg                                                                                      <  1:423043/89‑1 (MQ=255)
                                                                                                                         ccTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAg                                                                                      >  2:423043/1‑89 (MQ=255)
                                                                                                                          cTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGt                                                                          >  1:1944222/1‑100 (MQ=255)
                                                                                                                          cTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGt                                                                          <  2:1944222/100‑1 (MQ=255)
                                                                                                                          cTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGtt                        <  1:730491/150‑1 (MQ=255)
                                                                                                                          cTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGtt                        <  2:1913272/150‑1 (MQ=255)
                                                                                                                          cTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGtt                        >  1:1882003/1‑150 (MQ=255)
                                                                                                                          cTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGtt                        >  2:153171/1‑150 (MQ=255)
                                                                                                                           tCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGcaacaa                                                                                             <  2:1077918/80‑1 (MQ=255)
                                                                                                                           tCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGcaacaa                                                                                             >  1:1077918/1‑80 (MQ=255)
                                                                                                                           tCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCaacaac                                                                                            <  2:1907119/81‑1 (MQ=255)
                                                                                                                           tCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCaacaac                                                                                            >  1:1907119/1‑81 (MQ=255)
                                                                                                                           tCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTg                       >  2:2217375/1‑150 (MQ=255)
                                                                                                                           tCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTg                       >  2:1806255/1‑150 (MQ=255)
                                                                                                                              ggCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGATAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGc                    >  1:1395732/1‑150 (MQ=255)
                                                                                                                              ggCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGc                    >  2:1724244/1‑150 (MQ=255)
                                                                                                                                 cATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACTCTCAGGTTGTCGCTGTGCCGGTTGAGcaca                 <  1:1070774/150‑1 (MQ=255)
                                                                                                                                  aTAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGCACAt                >  1:1474834/1‑150 (MQ=255)
                                                                                                                                  aTAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGCACAt                <  2:1982689/150‑1 (MQ=255)
                                                                                                                                  aTAGAATCTCCTACAGGAGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCa                                           >  2:1650269/1‑123 (MQ=255)
                                                                                                                                  aTAGAATCTCCTACAGGAGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCa                                           <  1:1650269/123‑1 (MQ=255)
                                                                                                                                          tcCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTg                                                        >  2:2363776/1‑102 (MQ=255)
                                                                                                                                          tcCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTg                                                        <  1:2363776/102‑1 (MQ=255)
                                                                                                                                              aCAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGtt                                       <  1:724773/115‑1 (MQ=255)
                                                                                                                                              aCAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGtt                                       >  2:724773/1‑115 (MQ=255)
                                                                                                                                               cAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGCACCTCGCACAGTGCCTg   <  2:1199073/150‑1 (MQ=255)
                                                                                                                                                aGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGCACATCGCACAGTGCCTgg  <  1:1806255/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CTGCAACGCAGCCCGCCCAAACTGACGGGCAGTATCCAATGCTGGGGAGAAATTCTCTGGCACGCTGACTCGTACCTCGAGGTTGCTGAGAGAGGCCGATGGCGTTTTTACTGAAAATTTGCCTCCGGCCATAGAATCTCCTACAGGGGCTTGTGGTTTCCCTTTGAGAGCGCCGCTCAGCGTATCTGGTGCTTGTGCAACAACCGATAGAACGTCGGCCGTGAAATGCCCAGCACCTTGGCTGCCACGCTCAGGTTGTCGCTGTGCCGGTTGAGCACATCGCACAGTGCCTGG  >  NC_002947/3973015‑3973308

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