Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,068:1 +A 5.1% intergenic (+44/+28) 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,0681.A 5.1% 402.8 / 14.0 98intergenic (+44/+28)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 . (43/50);  new base A (2/3);  total (45/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGA  >  NC_002947/4749919‑4750211
                                                                                                                                                       |                                                                                                                                                
gcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAg                                                                                                                                                   >  1:639212/1‑150 (MQ=255)
    aaCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                             <  2:2383086/150‑1 (MQ=255)
      cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                           >  1:619831/1‑150 (MQ=255)
      cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCGACCAAAATTCATCGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                           <  2:663486/150‑1 (MQ=255)
       tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                          <  1:2029154/150‑1 (MQ=255)
           tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCttt                                                                                                                                      >  1:982560/1‑150 (MQ=255)
             cACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCttttt                                                                                                                                    >  2:479653/1‑150 (MQ=255)
               cTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgt                                                                                                                                  <  1:479653/150‑1 (MQ=255)
                tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                                 <  1:1630548/150‑1 (MQ=255)
                tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGACCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAA‑‑CCGGGGGCTTTttgtt                                                                                                                                 <  1:29668/150‑1 (MQ=38)
                tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCGGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTttgtt                                                                                                                                 <  2:1175627/150‑1 (MQ=38)
                 gCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgttg                                                                                                                                <  1:1760857/150‑1 (MQ=255)
                  cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTgg                                                                                                                               <  2:555923/150‑1 (MQ=255)
                     gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTg                                                                                                                            <  2:1191280/150‑1 (MQ=255)
                      aCCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGc                                                                                                                           <  1:984793/150‑1 (MQ=255)
                        cTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCa                                                                                                                         >  1:1931286/1‑150 (MQ=255)
                         tGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAt                                                                                                                        <  2:1931286/150‑1 (MQ=255)
                          gCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAtt                                                                                                                       <  2:338675/150‑1 (MQ=255)
                            gagaTCAAGGACCCGGCCCAGAACATCCCCTGGGTCACTCCTGGGGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTAt                                                                                                                     <  2:982560/147‑1 (MQ=255)
                            gACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGc                                                                                                                              >  1:1165840/1‑141 (MQ=255)
                            gACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGc                                                                                                                              <  2:1165840/141‑1 (MQ=255)
                               aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                                  <  1:1079207/150‑1 (MQ=255)
                                  aaGGACCCGGCCCAGAGCATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                               <  2:1223472/150‑1 (MQ=255)
                                  aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                               >  1:605320/1‑150 (MQ=255)
                                  aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                               >  1:2195208/1‑150 (MQ=255)
                                   aGGCCCCGCCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAg                                                                                                              <  2:1641619/150‑1 (MQ=255)
                                   aGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAg                                                                                                              <  2:2076470/150‑1 (MQ=255)
                                    ggACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCa                                                                                                                   <  2:47496/144‑1 (MQ=255)
                                    ggACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCa                                                                                                                   >  1:47496/1‑144 (MQ=255)
                                      aCCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTAc                                                                                                           <  1:294717/150‑1 (MQ=255)
                                             ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTc                                                                                                    >  1:800889/1‑150 (MQ=255)
                                             ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTc                                                                                                    >  1:1185526/1‑150 (MQ=255)
                                             ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTc                                                                                                    <  2:2422898/150‑1 (MQ=255)
                                               aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                             <  1:467344/107‑1 (MQ=255)
                                               aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                             >  2:467344/1‑107 (MQ=255)
                                                gAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTc                                                                                                 <  2:639212/150‑1 (MQ=255)
                                                   cATCACCTGGGTCACCCCTGATGGAGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAGTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  2:1207410/150‑1 (MQ=255)
                                                   cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  2:2195208/150‑1 (MQ=255)
                                                   cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  1:1050860/150‑1 (MQ=255)
                                                     tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                            <  2:1185526/150‑1 (MQ=255)
                                                         ccTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                         <  2:467993/150‑1 (MQ=255)
                                                           tGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                       <  2:699316/150‑1 (MQ=255)
                                                           tGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGCTCGCTTCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTGTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                       >  1:1239488/1‑150 (MQ=255)
                                                                cACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATg                                                                                  <  1:2107915/150‑1 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAGCAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                           <  2:2415423/90‑1 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAGCAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                           >  1:2415423/1‑90 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGGTGGCg                                                                               >  1:1541585/1‑150 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                               >  1:1059106/1‑150 (MQ=255)
                                                                   cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                               >  2:546608/1‑150 (MQ=255)
                                                                    ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                              >  2:1897581/1‑150 (MQ=255)
                                                                      tGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATc                                                                            >  2:2065097/1‑150 (MQ=255)
                                                                           gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                       <  1:2065097/150‑1 (MQ=255)
                                                                              cTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCt                                                                    <  2:356373/150‑1 (MQ=255)
                                                                               tGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                   <  2:605320/150‑1 (MQ=255)
                                                                               tGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                   >  1:831818/1‑150 (MQ=255)
                                                                                  cGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAGGAGTTCTGGATGGCGATCCCGCCGCTTGAg                                                                >  2:786191/1‑150 (MQ=255)
                                                                                   gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCAc                                                                                                       <  2:249999/110‑1 (MQ=255)
                                                                                   gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCAc                                                                                                       >  1:249999/1‑110 (MQ=255)
                                                                                     gCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAg                                                             <  1:1897581/150‑1 (MQ=255)
                                                                                         aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                              <  1:2277956/129‑1 (MQ=255)
                                                                                         aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                              >  2:2277956/1‑129 (MQ=255)
                                                                                           cTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCAc                                                       <  1:1678732/150‑1 (MQ=255)
                                                                                                cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACg                                                  >  2:2133917/1‑150 (MQ=255)
                                                                                                     aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACt                                             >  2:970355/1‑150 (MQ=255)
                                                                                                      tGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcag                                                                         >  1:1758547/1‑121 (MQ=255)
                                                                                                      tGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcag                                                                         <  2:1758547/121‑1 (MQ=255)
                                                                                                       gTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                                 <  1:1262530/112‑1 (MQ=255)
                                                                                                       gTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTgg                                           >  1:742704/1‑150 (MQ=255)
                                                                                                       gTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCGTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAATAGTTCTGGATgg                                                                                 >  2:1262530/1‑112 (MQ=255)
                                                                                                               ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAg                                                                <  2:83729/121‑1 (MQ=255)
                                                                                                               ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAg                                                                >  1:83729/1‑121 (MQ=255)
                                                                                                               ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtgg                                   >  1:502855/1‑150 (MQ=255)
                                                                                                                cTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCGGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtggt                                  >  2:1493846/1‑150 (MQ=255)
                                                                                                                   cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                               >  1:371905/1‑150 (MQ=255)
                                                                                                                     cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                                 <  2:2347456/98‑1 (MQ=255)
                                                                                                                     cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                                 >  1:2347456/1‑98 (MQ=255)
                                                                                                                     cTTGCGGCAACAAAAAGCAACCGGGTGTGGTGTG‑‑CCGGGGGGTGTTTGGTGGAAACCAGTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                             >  1:1872716/1‑150 (MQ=255)
                                                                                                                          ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGttt                        >  1:439648/1‑150 (MQ=255)
                                                                                                                           gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                       <  2:433830/150‑1 (MQ=255)
                                                                                                                           ccaacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                         >  2:419706/2‑84 (MQ=18)
                                                                                                                           ccaacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                         <  1:419706/83‑1 (MQ=18)
                                                                                                                              acaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                      <  2:1422886/148‑1 (MQ=38)
                                                                                                                              acaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                      >  1:1422886/1‑148 (MQ=38)
                                                                                                                               caaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGTTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATa                   >  1:2194441/1‑150 (MQ=255)
                                                                                                                                aaaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                    >  2:184405/1‑132 (MQ=255)
                                                                                                                                aaaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                    <  1:184405/132‑1 (MQ=255)
                                                                                                                                 aaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAAc                 <  2:800889/150‑1 (MQ=255)
                                                                                                                                  aaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACg                >  2:709875/1‑150 (MQ=255)
                                                                                                                                      cccccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGgag            >  1:1907591/1‑150 (MQ=255)
                                                                                                                                        cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg               >  2:983523/1‑145 (MQ=255)
                                                                                                                                        cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg               <  1:983523/145‑1 (MQ=255)
                                                                                                                                        cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa          >  1:1822364/1‑150 (MQ=255)
                                                                                                                                         ccGNTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTCGCCGTTTTCATAACGGGAGAAc         >  2:1504054/1‑150 (MQ=255)
                                                                                                                                         ccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAAc         <  2:1059106/150‑1 (MQ=255)
                                                                                                                                         ccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAAc         <  2:831818/150‑1 (MQ=255)
                                                                                                                                            gTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCa      <  2:619831/150‑1 (MQ=255)
                                                                                                                                                tcATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGa  >  2:100681/1‑150 (MQ=255)
                                                                                                                                                  aTGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAAc         <  2:2111647/141‑1 (MQ=255)
                                                                                                                                                  aTGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAAc         >  1:2111647/1‑141 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
GCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGA  >  NC_002947/4749919‑4750211

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