Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,506 C→G 6.2% intergenic (+23/‑109) tnaB → / → mdtL tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5060CG 6.2% 156.9 / 2.2 60intergenic (+23/‑109)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (2/2);  total (29/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 = 8.10e-01

TTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTT  >  NC_000913/3891359‑3891651
                                                                                                                                                   |                                                                                                                                                 
ttCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTTGCTAAATCCTTCAAGAAGCCAGCCATTCg                                                                                                                                                  >  2:989338/1‑149 (MQ=255)
  cccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCt                                                                                                                                                <  2:188550/149‑1 (MQ=255)
     aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                             <  1:353579/149‑1 (MQ=255)
     aaTCAGATGTTCACGGTCTGGGGCGGCAAGCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                             <  1:908949/149‑1 (MQ=255)
     aaTCAGATGTTAACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                             >  1:725103/1‑149 (MQ=255)
      aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                                                                            >  2:236214/1‑149 (MQ=255)
       tCAGATGTTAACTGTCTGGGGCGGCAATCTTATTCCTGCGATTGTAATTCTCTTTGGTATACCCGTGATTTTGAGCTGGTACGGCAACGTCTTTAACGTGTTACATAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTAGctt                                                                                                                                           >  2:36438/1‑149 (MQ=255)
         agagGTTCACGGTCTGGGGCGTCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATCACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTAAAGAAGCCAGCGAATGGCTGGcttct                                                                                                                                         <  2:136967/149‑1 (MQ=255)
          gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                        <  1:893230/149‑1 (MQ=255)
            tGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                                                                      <  2:53610/149‑2 (MQ=255)
              ttCTCGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTGATTCTCTTTTGTTTAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTATCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                    <  2:301289/149‑1 (MQ=255)
              ttCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                    >  1:275784/1‑149 (MQ=255)
                 aCGGTCTGGGGCGGCAATCTTAGTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                                                                 >  1:533486/1‑149 (MQ=255)
                  cGGTCTGGGGCGGCAAGCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                                                                                <  2:767816/149‑1 (MQ=255)
                        ggggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATGTTGTGCTGGGTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCACGAACCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAAt                                                                                                                          <  2:678010/149‑1 (MQ=255)
                          ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                        <  1:57617/149‑1 (MQ=255)
                          ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                        <  1:236214/149‑1 (MQ=255)
                               cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                                   <  1:1009343/149‑1 (MQ=255)
                                         ccGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                         >  1:247790/1‑149 (MQ=255)
                                         ccGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                         <  1:643452/149‑1 (MQ=255)
                                         ccGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGGCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                         <  1:163525/149‑1 (MQ=255)
                                         ccGGCGATTGTCATTCTCTTTGGTATAACCGCGTTTTTGTGCTGGTTCGGCCACGGCTTTAACGTGATACCTAAATGTGGCTAAATCCTTCAATAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                         <  2:668518/149‑1 (MQ=255)
                                          cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                                                                        <  2:172683/149‑1 (MQ=255)
                                           ggCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCTGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTGATGTCCAAATGGCa                                                                                                       >  2:589472/1‑149 (MQ=255)
                                           ggCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATACTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                                                                       >  2:569810/1‑149 (MQ=255)
                                                      ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTAGCTGGCTTCTTTCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                            <  1:703561/149‑1 (MQ=255)
                                                        ctctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCAGGAtcc                                                                                          <  2:777544/149‑1 (MQ=255)
                                                          ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctc                                                                                        <  2:725103/149‑1 (MQ=255)
                                                           tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                       <  2:292291/149‑1 (MQ=255)
                                                            ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATcctcct                                                                                      >  1:321831/1‑149 (MQ=255)
                                                                     ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGt                                                                             >  1:621646/1‑149 (MQ=255)
                                                                              tgatcTGGGTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCCCGTATGCTTTGc                                                                    <  2:449403/145‑1 (MQ=255)
                                                                                tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGt                                                                  >  1:524798/1‑149 (MQ=255)
                                                                                 tcTGGTTCGGCAACTTCTTTAAGGTGTTAACTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGATTCTTGCCTCTCAGGAAATCACTAATGTCCAAATGTCAACTCGCCTGATCCTCCTTCACCCCGTATGCTTTGCGTc                                                                 <  1:564535/148‑1 (MQ=255)
                                                                                  cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCTAATGGCAACTCGCCTGATGCTCCTGCACCACGTATGCTTTGCggcc                                                                >  2:996911/1‑146 (MQ=255)
                                                                                   tGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCATTGAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc                                                               <  1:36438/149‑1 (MQ=255)
                                                                                          gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCAGATCATCCTTCACCACGTATGCTTTGCGTCACCTTACTa                                                        >  2:281978/1‑149 (MQ=255)
                                                                                            aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                      >  1:337183/1‑149 (MQ=255)
                                                                                              cGTCTTTAACGTGTTAACTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg                                                    <  2:321831/149‑1 (MQ=255)
                                                                                              cGTCATTAACGTGTTACCTAAATTAGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg                                                    <  2:272825/149‑1 (MQ=255)
                                                                                               gTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAgg                                                   >  1:97834/1‑149 (MQ=255)
                                                                                                 cTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAc                                                 >  1:377952/1‑149 (MQ=255)
                                                                                                     aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTGATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGcgt                                             >  2:902748/1‑147 (MQ=255)
                                                                                                         tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGc                                         <  2:533486/149‑1 (MQ=255)
                                                                                                             aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCGTTAGCCCAt                                     >  2:100964/1‑149 (MQ=255)
                                                                                                                       ggCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttt                           >  2:779134/1‑149 (MQ=255)
                                                                                                                        gCTAAATCCTTCAATAAGGCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCCCGTATGCTTTGCTTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                          <  1:91101/149‑1 (MQ=255)
                                                                                                                        gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTAACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                          >  2:965805/1‑149 (MQ=255)
                                                                                                                         cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCCCCACGTATGCTTTGCGTCACCTTACTATCAGGCCGCTTTAGCGCCTGTCCCGCttttt                         >  2:238200/1‑149 (MQ=255)
                                                                                                                         cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                         >  1:4222/1‑149 (MQ=255)
                                                                                                                         cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTGCTTGCCGCTCAGGAACTCCCTTATGTCCAAATGGCAACTCGCCTGCTCCTCCTGCACCGCGTCTGCGTTGCGTCGCCTTGCTATCAGGGCGCGGTGGCCCATGTCCCGCtttgg                         >  2:668164/1‑147 (MQ=255)
                                                                                                                          tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCCCTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                        >  2:183863/1‑149 (MQ=255)
                                                                                                                              tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                    <  1:281792/149‑1 (MQ=255)
                                                                                                                              tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                    <  1:26527/149‑1 (MQ=255)
                                                                                                                                attCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACATCGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTACCGCTTTTTGATTTGt                  <  1:930030/148‑1 (MQ=255)
                                                                                                                                  gcAAGAAGCCAGCGAATGGCTGGCTTCTAGCCTCTCAGGAAATCCCTTAAGTCCAAATCGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCATTTGCCCATGTCCCGCTTTTTGATTTTTat                >  2:383551/2‑148 (MQ=255)
                                                                                                                                      gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt            >  1:44599/1‑149 (MQ=255)
                                                                                                                                       aagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTg           >  1:167910/1‑149 (MQ=255)
                                                                                                                                        agccagccATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGc          <  1:848018/149‑1 (MQ=255)
                                                                                                                                                aTTCGCTGGCATCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGttt  >  2:547231/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTT  >  NC_000913/3891359‑3891651

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