Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,502 C→G 8.4% intergenic (+19/‑113) 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,5020CG 8.4% 165.3 / 9.1 60intergenic (+19/‑113)tnaB/mdtLtryptophan:H(+) symporter TnaB/efflux pump MdtL
Reads supporting (aligned to +/- strand):  ref base C (25/30);  new base G (2/3);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-02

TCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTT  >  NC_000913/3891360‑3891641
                                                                                                                                              |                                                                                                                                            
tCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAagccaccc                                                                                                                                              >  2:316802/1‑143 (MQ=255)
 cccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAgccagcca                                                                                                                                             >  1:430762/1‑143 (MQ=255)
   cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAGgccagccagt                                                                                                                                           >  1:379598/1‑141 (MQ=255)
       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)
                                                                                                                                              |                                                                                                                                            
TCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATG‑CTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTT  >  NC_000913/3891360‑3891641

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: ‑