Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,059 Δ1 bp 6.5% intergenic (+35/+37) PP_4203 → / ← mqsA ubiquinone oxidoreductase electron transfer flavoprotein/type II toxin‑antitoxin MqsRA system antitoxin/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,750,0590T. 6.5% 332.4 / 17.0 107intergenic (+35/+37)PP_4203/mqsAubiquinone oxidoreductase electron transfer flavoprotein/type II toxin‑antitoxin MqsRA system antitoxin/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (49/50);  new base . (3/4);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.34e-01

Only 100 of 109 total aligned reads displayed.

CGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAT  >  NC_002947/4749918‑4750208
                                                                                                                                             |                                                                                                                                                        
cgcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGa                                                                                                                                                  <  2:1564520/150‑1 (MQ=255)
 gcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAg                                                                                                                                                 <  1:681119/150‑1 (MQ=255)
    gAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCg                                                                                                                                            >  2:786342/1‑150 (MQ=255)
     aaCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                           <  2:1814163/150‑1 (MQ=255)
     aaCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                           <  2:741274/150‑1 (MQ=255)
     aaCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                           <  2:2349603/150‑1 (MQ=255)
      aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACCAGCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                          <  2:1326906/150‑8 (MQ=255)
       cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                         <  1:1254380/150‑1 (MQ=255)
        tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCAt                                                                                                                                                    <  1:2154439/140‑1 (MQ=255)
        tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCAt                                                                                                                                                    >  2:2154439/1‑140 (MQ=255)
        tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                        <  2:170800/150‑1 (MQ=255)
                 tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                              >  1:1109485/1‑150 (MQ=255)
                 ggCAAGACCTGCGACAGCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTCTctgtt                                                                                                                              >  2:833661/2‑150 (MQ=255)
                    aaGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGCCAACAAAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTttgttgcc                                                                                                                           <  2:1177560/150‑3 (MQ=37)
                      gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                        <  2:1136309/136‑1 (MQ=255)
                      gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                        >  1:1136309/1‑136 (MQ=255)
                      gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccgg                                                                                                                                           <  2:963808/133‑7 (MQ=255)
                      gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccgg                                                                                                                                           >  1:963808/1‑127 (MQ=255)
                         cTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCGTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCa                                                                                                                      <  1:698836/150‑1 (MQ=255)
                          tGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAt                                                                                                                     >  2:2014092/1‑150 (MQ=255)
                          tGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAt                                                                                                                     >  1:1801495/1‑150 (MQ=255)
                               cATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGTGAACAAAAANCCCCCGNTCCTCATGAG‑‑CCGTGGGCTTTTTGTTGGCTGCCATTATCa                                                                                                                <  1:1260787/150‑1 (MQ=255)
                                aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                               >  2:377636/1‑150 (MQ=255)
                                   aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                              >  1:2262692/1‑132 (MQ=255)
                                   aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                            >  2:855820/1‑150 (MQ=255)
                                     ggACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGt                                                                                                          >  1:1563220/1‑150 (MQ=255)
                                         ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCg                                                                                                                                            <  1:812595/113‑1 (MQ=255)
                                         ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCg                                                                                                                                            >  2:812595/1‑113 (MQ=255)
                                         ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGc                                                                                                      >  2:1694240/1‑150 (MQ=255)
                                          cGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATc                                                                                                                 <  1:1660832/138‑1 (MQ=255)
                                          cGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATc                                                                                                                 >  2:1660832/1‑138 (MQ=255)
                                              ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggg                                                                                                                                          >  1:1600049/1‑110 (MQ=255)
                                              ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggg                                                                                                                                          <  2:1600049/110‑1 (MQ=255)
                                               cAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                        <  2:1162996/111‑1 (MQ=255)
                                               cAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                        >  1:1162996/1‑111 (MQ=255)
                                                    cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGc                                                                                                                        >  1:461390/1‑121 (MQ=255)
                                                    cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGc                                                                                                                        <  2:461390/121‑1 (MQ=255)
                                                    cATCACCTGGGTCAACCCTGAAGGCGCTGGCGGGACGAACTATCCAAACATGTAAGGCCCCTGCACTTGCAGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  2:1109485/150‑1 (MQ=255)
                                                      tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                         >  2:1805486/1‑150 (MQ=255)
                                                        aCCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGt                                                                                       >  1:2211898/1‑150 (MQ=255)
                                                        aCCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGt                                                                                       >  1:1842924/1‑150 (MQ=255)
                                                           tGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCACGGTTCTCATGAG‑‑CCGCCATCGTTGGTGTTGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                    <  1:2224215/150‑1 (MQ=255)
                                                               tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           >  1:269608/1‑139 (MQ=255)
                                                               tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  2:269608/139‑1 (MQ=255)
                                                                 aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                              <  1:1605472/150‑1 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTAc                                                                                                        >  1:1744144/1‑122 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                            >  2:1275994/1‑150 (MQ=255)
                                                                    ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTa                                                                                                         <  1:710245/120‑1 (MQ=255)
                                                                    ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTa                                                                                                         >  2:710245/1‑120 (MQ=255)
                                                                     cTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTc                                                                                     >  1:399578/1‑139 (MQ=255)
                                                                     cTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTc                                                                                     <  2:399578/139‑1 (MQ=255)
                                                                     cTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGAt                                                                          <  2:1842924/150‑1 (MQ=255)
                                                                      tGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATc                                                                         <  1:1275994/150‑1 (MQ=255)
                                                                           gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                    <  1:1009310/150‑1 (MQ=255)
                                                                           gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                    <  1:701103/150‑1 (MQ=255)
                                                                             gcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCagcagc                                                                  >  1:896331/1‑150 (MQ=255)
                                                                               tGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGACAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                <  2:1738656/150‑1 (MQ=255)
                                                                                ggCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTg                                                               >  1:468007/1‑150 (MQ=255)
                                                                                   gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGc                                                            >  1:1043134/1‑150 (MQ=255)
                                                                                     gCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTACCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAg                                                          <  2:1563220/150‑1 (MQ=255)
                                                                                      ccGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                         <  1:855820/150‑1 (MQ=255)
                                                                                         aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAACCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTc                                                      >  2:695615/1‑150 (MQ=255)
                                                                                              tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgttg                                                                                                                             <  1:836399/74‑1 (MQ=255)
                                                                                              tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgttg                                                                                                                             >  2:836399/1‑74 (MQ=255)
                                                                                                 gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGc                                              >  1:1455443/1‑150 (MQ=255)
                                                                                                 gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGATTCACCAACGc                                              >  1:1766275/1‑150 (MQ=255)
                                                                                                      tGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTg                                         >  2:837635/1‑150 (MQ=255)
                                                                                                         aaGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCg                                      >  2:382067/1‑150 (MQ=255)
                                                                                                         aaGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGATCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCg                                      >  2:47474/1‑150 (MQ=255)
                                                                                                             ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggt                                  <  1:1694240/150‑1 (MQ=255)
                                                                                                                   cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCa                                            >  1:1819077/1‑134 (MQ=255)
                                                                                                                   cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCa                                            <  2:1819077/134‑1 (MQ=255)
                                                                                                                    ccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTg                           >  1:671263/1‑150 (MQ=255)
                                                                                                                      ttGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                         >  1:973240/1‑86 (MQ=255)
                                                                                                                          ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGATTCACCAACGCCACTGGCGGGGTGGTCTTGCCGttt                     <  2:1766275/150‑1 (MQ=255)
                                                                                                                          gccaacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGttt                     >  1:1571872/3‑150 (MQ=38)
                                                                                                                           gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                    >  1:40761/1‑150 (MQ=255)
                                                                                                                           gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                    >  2:576954/1‑150 (MQ=255)
                                                                                                                            caacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                   >  2:185854/1‑150 (MQ=255)
                                                                                                                            caacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                   <  1:576954/150‑1 (MQ=255)
                                                                                                                            caacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                              >  2:1951165/1‑107 (MQ=25)
                                                                                                                            caacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACCTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                              <  1:1951165/107‑1 (MQ=25)
                                                                                                                             aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                              >  2:1541535/1‑106 (MQ=255)
                                                                                                                             aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                              <  1:1541535/106‑1 (MQ=255)
                                                                                                                             aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCa                  >  1:2037668/1‑150 (MQ=255)
                                                                                                                                 aaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAAc              >  2:1827924/1‑150 (MQ=255)
                                                                                                                                  aaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           >  2:1534060/1‑72 (MQ=255)
                                                                                                                                  aaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  1:1534060/72‑1 (MQ=255)
                                                                                                                                  aaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACg             >  2:70810/1‑150 (MQ=255)
                                                                                                                                   aaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg            <  2:468007/150‑1 (MQ=255)
                                                                                                                                   aaGCCCCCGGTTCTCATGAG‑‑CCGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACggg           <  2:303215/150‑1 (MQ=255)
                                                                                                                                       ccccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGgaga        >  1:155907/1‑150 (MQ=255)
                                                                                                                                        cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCa                  >  1:317738/1‑139 (MQ=255)
                                                                                                                                        cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCa                  <  2:317738/139‑1 (MQ=255)
                                                                                                                                        cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa       <  1:2014092/150‑1 (MQ=255)
                                                                                                                                         ccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGa                                                                                 <  1:1297500/75‑1 (MQ=255)
                                                                                                                                         ccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAAc      <  2:155907/150‑1 (MQ=255)
                                                                                                                                          cGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACg     >  1:1716222/1‑150 (MQ=255)
                                                                                                                                            gTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGCACgaa   <  1:695615/150‑3 (MQ=255)
                                                                                                                                             ttCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAt  <  2:2211898/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                        
CGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAT  >  NC_002947/4749918‑4750208

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: