Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,504 T→A 7.8% intergenic (+21/‑111) tnaB → / → mdtL tryptophan:H(+) symporter TnaB/efflux pump MdtL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5040TA 7.8% 180.6 / 7.1 64intergenic (+21/‑111)tnaB/mdtLtryptophan:H(+) symporter TnaB/efflux pump MdtL
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base A (2/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.86e-01

CAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCT  >  NC_000913/3891367‑3891646
                                                                                                                                         |                                                                                                                                               
cAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGCAGCCAGCCATTCACt                                                                                                                                            >  1:363950/1‑143 (MQ=255)
cAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCt                                                                                                                                            <  2:540089/143‑1 (MQ=255)
       tcctggTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTGTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                                                                     <  2:256275/139‑1 (MQ=255)
       tCACGGGCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                                                                     <  2:94280/143‑1 (MQ=255)
          cGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttcttcc                                                                                                                                  >  1:542163/1‑141 (MQ=255)
            gTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                                                                <  1:486365/143‑4 (MQ=255)
            gTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                <  1:530516/143‑1 (MQ=255)
            gTCTGGGGCGGCAATCTTAGTCCGTTGATTGTTATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCTGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                <  2:237163/143‑1 (MQ=255)
                ggggCTGCAATCTTATTTCGTCGGTTGTCATTGTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCGTTAACGTGTTACCTATATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                                                                            <  2:100031/143‑1 (MQ=255)
                   gcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                                                                         >  1:407830/1‑143 (MQ=255)
                   gcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                                                                         <  2:614118/143‑1 (MQ=255)
                     ggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                                                                       <  1:34904/143‑1 (MQ=255)
                     ggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                                                                       >  1:515793/1‑143 (MQ=255)
                          tCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                  <  1:332262/143‑1 (MQ=255)
                           ctcaTTCCGGCGATTGTCATTCCCTTTGGTATACCCGTGATTTCGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTCCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                                                                 <  1:648266/140‑1 (MQ=255)
                             tatggcGGCGATTGTCATTCTCTCTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                               <  2:524886/138‑1 (MQ=255)
                             tATTCCGGCGATTGTCATTCTCTTTGGTATACCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                               <  2:335780/143‑1 (MQ=255)
                             gaTTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                               <  1:376712/142‑1 (MQ=255)
                                  cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCGAAATGCTTCGAGAAGCCAGCCATCCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                                                                          >  1:126503/1‑143 (MQ=255)
                                    gCGATTGTCATTCTCTTTGGTATAACAGTGCTTTTGGGCTGGCTCGGCAACGTCTTTAACGTGTTACCTGAATTAGGCTGAACGCTTCGAGAAGCCAGCCATTCGCTTGCGCCTTACCCGTCAGGAAATCACGGATGTCCaaa                                                                                                        >  2:12512/1‑143 (MQ=255)
                                     gcaTTGTCGGTCTCTTTGTGATAGCCTTGCTTTTGTGCTTGTTCGGCGACGTGTTTACCGTGTTACCTACATTTGGCTCAATCCTTCAAGAAGCCACCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCatat                                                                                                       <  2:576509/141‑1 (MQ=255)
                                       tttGTCTTTTTGTTTTGTTTTACCGTTGTTTTGTTCTGGTTCGGCCTCGTCTTTAACGTGTTTCCTGATTTTGGCTAAATCCTTCAAGTAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                     <  2:417203/142‑1 (MQ=255)
                                          gTCATTCTCTGTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCaa                                                                                                  <  1:242580/143‑1 (MQ=255)
                                            cATTGTCTTTGGTATAGCCGTGATTTTGTGCTGGTTCGGCAACGGCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCATCt                                                                                                <  2:632429/143‑1 (MQ=255)
                                             aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                               <  1:379126/143‑1 (MQ=255)
                                                ctctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                                                                            <  1:378903/143‑1 (MQ=255)
                                                 tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCACATGGCCACTCGCCt                                                                                           >  1:227596/1‑143 (MQ=255)
                                                   tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTgc                                                                                         >  1:375484/1‑142 (MQ=255)
                                                   tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                         >  1:259492/1‑143 (MQ=255)
                                                    ttGGTATAACCGTGATTTTGTGCTGCTTCGGCTACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                        >  2:83621/1‑143 (MQ=255)
                                                     tGGTATAACTGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGAAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                                                                       <  2:262623/143‑1 (MQ=255)
                                                        tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctc                                                                                    <  2:430762/143‑1 (MQ=255)
                                                        gataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctc                                                                                    <  2:93474/142‑1 (MQ=255)
                                                                     ttctgCTTGTTCGGCAACGTCTTTACCGTGTTCCCTAAACTTGGCTAAATCCTTCATGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCACGAAATCACTTATGTCCAAATTGCAACTCGCCTGATCCTCCTTCACCACGTAt                                                                       <  2:281903/140‑1 (MQ=255)
                                                                            ggTTCGGCAACGTCTTTAACGTGTTACCTCAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGc                                                               <  1:10742/143‑1 (MQ=255)
                                                                             gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAACCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTCCg                                                              >  2:494513/1‑143 (MQ=255)
                                                                              ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTGAGGGCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGt                                                             <  1:143323/143‑1 (MQ=255)
                                                                               tCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTAGCCTGATCCTCCTTCACCACGTATG‑CTTTGAGTc                                                            >  2:491540/1‑143 (MQ=255)
                                                                                     aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTGAAGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTa                                                      >  2:324949/1‑143 (MQ=255)
                                                                                       gTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACt                                                    >  1:52116/1‑143 (MQ=255)
                                                                                        tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTAAGGAAATCACGAATGCCCAAAGGGCAACCCGCCTGATCCTCCTTCACCACGTATG‑CTTCGCGTCACCTCAcca                                                   >  2:12309/1‑141 (MQ=255)
                                                                                          tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATc                                                 <  1:379261/143‑1 (MQ=255)
                                                                                            tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATAGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CATTCCGTCACTTTACTGTCAg                                               >  2:2715/1‑143 (MQ=255)
                                                                                             aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAGCTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAgg                                              <  1:678074/143‑1 (MQ=255)
                                                                                               cGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGAc                                            <  1:571658/143‑1 (MQ=255)
                                                                                                 tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTTCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGCTCCTCCTTCACCACGTATG‑CTTTGCGACACCTCACTATCAGGAAGc                                          >  2:559078/1‑143 (MQ=255)
                                                                                                      ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATGCTACCAC‑CGAACTACGTCTGTGCGTCACCTTAATATCACGACGCtttgg                                     >  2:139933/1‑141 (MQ=255)
                                                                                                       cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGCATT‑CTTTCCGTCACCTTACTATCAGGACGTCTTAGc                                    >  1:244979/1‑143 (MQ=255)
                                                                                                       cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTGATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGc                                    >  2:501643/1‑143 (MQ=255)
                                                                                                       cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTACCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCGTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGc                                    >  2:364004/1‑143 (MQ=255)
                                                                                                        tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTCCGTCACCTTACTCTCAGCACGCtttcccc                                   >  1:474859/1‑139 (MQ=255)
                                                                                                         gactttGGCTAGATCCTTCACGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAACTGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGccc                                  <  2:20012/140‑1 (MQ=255)
                                                                                                           aTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCAt                                <  1:277010/143‑1 (MQ=255)
                                                                                                                  tAAATCCTTCAAGAACCCAGGCATTCTCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCATCACCTTACTATCAGGACGCTTTAGCCCCTGTCCCGc                         >  2:272383/1‑143 (MQ=255)
                                                                                                                        cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                   >  2:279555/1‑143 (MQ=255)
                                                                                                                         ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCGTTAGCCCACGTCCCGCTTTTCGa                  >  2:372225/1‑143 (MQ=255)
                                                                                                                          tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAt                 <  1:635073/143‑1 (MQ=255)
                                                                                                                          gcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCGAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAt                 >  1:352565/2‑143 (MQ=255)
                                                                                                                           cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGTTCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAtt                <  1:607181/143‑1 (MQ=255)
                                                                                                                              gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGCCCAAATGGCAACTCGCCTGATCCTCCTTCACAACATATG‑ATTTGCGTCACGTTACTATCAGGACGCTTTAGGCCAGGTCACGCTTTTTGATTTgg             >  1:267155/1‑142 (MQ=255)
                                                                                                                                   cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCGACTCGCCTGATCCTCCTGCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGttt        <  2:395218/143‑1 (MQ=255)
                                                                                                                                    agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCTCGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt       <  1:243510/143‑1 (MQ=255)
                                                                                                                                    agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTGTCAGGGCGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt       >  1:281541/1‑143 (MQ=255)
                                                                                                                                        aTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCCCCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTGCCGCTTCTTGATTTGTAGCTTTGccg   >  2:44501/1‑142 (MQ=255)
                                                                                                                                         ttCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTGGATGTGTAGTTTTGCCCt  >  1:395157/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                               
CAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCT  >  NC_000913/3891367‑3891646

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