Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,520,390 G→A 5.3% A68A (GCG→GCA cspD → DNA replication inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,520,3900GA 5.3% 220.2 / 2.2 75A68A (GCG→GCAcspDDNA replication inhibitor
Reads supporting (aligned to +/- strand):  ref base G (42/29);  new base A (2/2);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

GTAAATCCGATGACCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGAT  >  NC_002947/4520251‑4520539
                                                                                                                                           |                                                                                                                                                     
gTAAATCCGATGACCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCgg                                                                                                                                             <  1:1135895/150‑1 (MQ=255)
gTAAATCCGATGACCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCgg                                                                                                                                             >  1:500710/1‑150 (MQ=255)
    aTCCGATGACCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCaga                                                                                                                                         <  2:903796/150‑1 (MQ=255)
            aCCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCa                                                                                                                                 <  1:1297988/150‑1 (MQ=255)
            aCCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCAGCTGTGGCAGAGGCTGCCa                                                                                                                                 <  2:80665/150‑1 (MQ=255)
              cTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCAGCGGCCGCGGCatctt                                                                                                                                       >  2:185121/1‑138 (MQ=255)
              cTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCAGCGGCCGCGGCatctt                                                                                                                                       <  1:185121/142‑5 (MQ=255)
                   cGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCtt                                                                                                                          <  1:1546774/150‑1 (MQ=255)
                    gCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGACACCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTg                                                                                                                         >  2:1910239/1‑150 (MQ=255)
                     cccATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGTGGCAGAGGCTGCCAGAGCCTTgg                                                                                                                        >  1:913574/1‑150 (MQ=255)
                     cccATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCaga                                                                                                                                         >  2:370218/1‑133 (MQ=255)
                     cccATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCaga                                                                                                                                         <  1:370218/133‑1 (MQ=255)
                      ccATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTggg                                                                                                                       >  1:1369048/1‑150 (MQ=255)
                      ccATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTggg                                                                                                                       >  1:990203/1‑150 (MQ=255)
                             tCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCg                                                                                                                <  2:1105937/150‑1 (MQ=255)
                                gCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCagag                                                                                                                                        <  2:162017/123‑1 (MQ=255)
                                gCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCagag                                                                                                                                        >  1:162017/1‑123 (MQ=255)
                                 ccATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCaaa                                                                                                            >  2:1730004/1‑150 (MQ=255)
                                   aTCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCAGCGGCCGCGGCatctt                                                                                                                                       <  2:1409939/121‑5 (MQ=255)
                                   aTCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCAGCGGCCGCGGCatctt                                                                                                                                       >  1:1409939/1‑117 (MQ=255)
                                       aaagtgACGGACACAAGCCCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGaaa                                                                                                      <  1:140355/145‑1 (MQ=255)
                                       aaaTGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGaaa                                                                                                      >  2:126845/1‑150 (MQ=255)
                                       aaaTGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGaaa                                                                                                      >  1:1514822/1‑150 (MQ=255)
                                            gACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCa                                                                                                 >  1:143269/1‑150 (MQ=255)
                                            gACGCATAAAAAACCCTGACGCCGGCCCAGTCCGTGGTATTAAACATCGGCCAGGGCCCCAAGGGCTCGCACACCGTTGAAATCAAAGCTGCCGCGGCCGCTGCGGCAGAGCCTGCCAGAGCTTTAGGCGAGCCACAAAGCGAAACTGCa                                                                                                 >  2:1540237/1‑150 (MQ=255)
                                             aCGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCaa                                                                                                >  2:1333641/1‑150 (MQ=255)
                                               ggATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCaaaa                                                                                              >  2:1883803/1‑150 (MQ=255)
                                               ggATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATACAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAGATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCaaaa                                                                                              >  1:292178/1‑150 (MQ=255)
                                                 aTACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGc                                                                                            >  2:1734578/1‑150 (MQ=255)
                                                 aTACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGc                                                                                            >  2:129438/1‑150 (MQ=255)
                                                   aCAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTTGGCGAGCCGCAAAGCGAAACTGCAAAAGCCt                                                                                          <  1:1910239/150‑1 (MQ=255)
                                                     aaGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGTGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGa                                                                                        <  2:913574/150‑1 (MQ=255)
                                                          ccTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAgg                                                                                                                                       >  1:1445770/1‑98 (MQ=255)
                                                          ccTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAgg                                                                                                                                       <  2:1445770/98‑1 (MQ=255)
                                                          ccTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACt                                                                                                    <  1:1751356/133‑1 (MQ=255)
                                                          ccTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACt                                                                                                    >  2:1751356/1‑133 (MQ=255)
                                                          ccTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGACCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGActgcc                                                                                   >  2:559930/1‑150 (MQ=255)
                                                             gAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGActgcgtgc                                                                                >  1:944446/1‑150 (MQ=255)
                                                                 gCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCaaaa                                                                                              <  1:639530/132‑1 (MQ=255)
                                                                 gCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCaaaa                                                                                              >  2:639530/1‑132 (MQ=255)
                                                                    ggCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGc                                                                         <  2:1663079/150‑1 (MQ=255)
                                                                       cAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGa                                                                      >  1:1526867/1‑150 (MQ=255)
                                                                        aGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGAt                                                                     <  2:500710/150‑1 (MQ=255)
                                                                          gCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAg                                                                   >  1:1238895/1‑150 (MQ=255)
                                                                              tGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACaa                                                               >  2:954401/1‑150 (MQ=255)
                                                                                gTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGTCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACaaaa                                                             >  1:660951/1‑150 (MQ=255)
                                                                                   ttCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACaaaaaaa                                                          <  1:1333641/150‑1 (MQ=255)
                                                                                   ttCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGAAAGACaaaaaaa                                                          >  1:1267763/1‑150 (MQ=255)
                                                                                      aaCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGGTGCGGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACCAAAAAATCg                                                       >  1:538493/1‑150 (MQ=255)
                                                                                      aaCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCg                                                       <  1:768616/150‑1 (MQ=255)
                                                                                      aaCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCg                                                       >  2:986850/1‑150 (MQ=255)
                                                                                      aaCATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCg                                                       >  2:768616/1‑150 (MQ=255)
                                                                                        cATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGc                                                     >  2:691100/1‑150 (MQ=255)
                                                                                                   gCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTaa                                          <  1:126845/150‑1 (MQ=255)
                                                                                                     cccAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTaaaa                                        >  2:558093/1‑150 (MQ=255)
                                                                                                       cAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCtgt                                                                           >  1:1848291/1‑113 (MQ=255)
                                                                                                       cAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCtgt                                                                           <  2:1848291/113‑1 (MQ=255)
                                                                                                        aaGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaa                                     <  1:762827/150‑1 (MQ=255)
                                                                                                        aaGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaa                                     <  1:1260028/150‑1 (MQ=255)
                                                                                                         aGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCt                                                                                          >  2:698716/1‑96 (MQ=255)
                                                                                                         aGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCt                                                                                          <  1:698716/96‑1 (MQ=255)
                                                                                                         aGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAaagaag                                    >  1:1255673/1‑150 (MQ=255)
                                                                                                              ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCg                                                   <  1:1695033/130‑1 (MQ=255)
                                                                                                              ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCg                                                   >  2:1695033/1‑130 (MQ=255)
                                                                                                              ttGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGAcc                               >  2:1284119/1‑150 (MQ=255)
                                                                                                                      cGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCt                       <  1:1284119/150‑1 (MQ=255)
                                                                                                                           aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGc                  >  2:296364/1‑150 (MQ=255)
                                                                                                                              tCAAAGATGCCGCGGCCGCGGCGGCAGAGGCGGCTAGGGCCGTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCtt               <  1:1883803/150‑1 (MQ=255)
                                                                                                                               cAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTa              >  1:1584329/1‑150 (MQ=255)
                                                                                                                                aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGtttt                         <  1:1387880/138‑1 (MQ=255)
                                                                                                                                aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGtttt                         >  2:1387880/1‑138 (MQ=255)
                                                                                                                                      gCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATgtgt       >  1:242059/1‑150 (MQ=255)
                                                                                                                                       ccGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTaaaa                                        <  1:1336493/116‑1 (MQ=255)
                                                                                                                                       ccGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTaaaa                                        >  2:1336493/1‑116 (MQ=255)
                                                                                                                                       ccGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt      <  2:292178/150‑1 (MQ=255)
                                                                                                                                          cgGCCGCTGCGGCAGAGGCTGCCAGATCTTTTGGCGAGCCGCAAAGCTAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtttaa   <  1:559930/150‑3 (MQ=255)
                                                                                                                                          cgGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGa   <  1:1734578/150‑1 (MQ=255)
                                                                                                                                           ggCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGAt  <  2:242059/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
GTAAATCCGATGACCTGTTCGCCCATTATTCCGCCATCCAAATGGACGGATACAAGACCCTGAAGGCTGGCCAGGCCGTGGTATTCAACATCGTCCAGGGCCCCAAGGGCTTGCACGCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGAT  >  NC_002947/4520251‑4520539

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