Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,504 T→A 5.3% intergenic (+21/‑111) 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,5040TA 5.3% 347.1 / 6.0 122intergenic (+21/‑111)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (62/52);  new base A (3/4);  total (66/56)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 122 total aligned reads displayed.

TTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGT  >  NC_000913/3891358‑3891649
                                                                                                                                                  |                                                                                                                                                 
     cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCAATCGCTgg                                                                                                                                            >  2:711936/1‑149 (MQ=255)
      tttcttaTGTTCTCGGTCTGGGGAGGTATTCTTATTCCGGCGATTGTCATTCTATTTGGTATAAACGTGATTTTGTGCTGGTTCGGCAAAGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAAtggctggc                                                                                                                                           <  1:1017667/143‑1 (MQ=255)
      aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                           >  1:786069/1‑149 (MQ=255)
      aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                           <  2:1278900/149‑1 (MQ=255)
      aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                           <  1:488689/149‑1 (MQ=255)
      aaTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGAGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                                                                           >  2:1097874/1‑149 (MQ=255)
       aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                                                                          >  2:711479/1‑149 (MQ=255)
        tCAGATGTGCACGGTCTGGGGCGGCTATCTTATTACGGAGATTGGCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGGCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGctt                                                                                                                                         <  1:1057728/149‑1 (MQ=255)
         cAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttc                                                                                                                                        <  2:375108/149‑1 (MQ=255)
          aGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGGGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                                                                       >  1:182216/1‑149 (MQ=255)
           gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                      <  2:406579/149‑1 (MQ=255)
           gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCAGTCGCTGGcttctt                                                                                                                                      >  1:188482/1‑149 (MQ=255)
              gTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                                                                   <  1:259665/149‑3 (MQ=255)
              gTTCACGGGCTGGGGCGGCAAACTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTAGGCTAAATCCTTCAAGAAGCCAGCGAAAGGCTGGCTTCTTgaa                                                                                                                                   <  1:499030/149‑3 (MQ=255)
               ttCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAGCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                  <  2:827704/149‑1 (MQ=255)
                  aCGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                                                               <  1:711479/149‑1 (MQ=255)
                  aCGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                                                               >  1:518947/1‑149 (MQ=255)
                  aCGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                                                               >  1:1068929/1‑149 (MQ=255)
                  aCGGTCTGGGGCGGCAATCATATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCTGCAACGTCTTTAACGTGTTACCTAAATTTGGATAAATCCTTCAAGAAGCCAGCCATTCGATAGCTTCTAGCctcta                                                                                                                               >  1:645421/1‑148 (MQ=255)
                       cTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                                                                          <  1:526601/149‑1 (MQ=255)
                           ggcggcAATGTTATTGCGGCGATTGTCATTCTCTTTGGTATTACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                      <  2:1549190/149‑1 (MQ=255)
                           ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGTTATATACGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCTCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                      <  2:563190/149‑1 (MQ=255)
                           ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGAGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                      <  2:70488/149‑1 (MQ=255)
                              ggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                    <  1:151676/148‑1 (MQ=255)
                               gcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTa                                                                                                                  >  2:1135723/1‑149 (MQ=255)
                                cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCTGGATTTTGTGCTGTTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                                                                 <  1:899731/149‑1 (MQ=255)
                                  aTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATga                                                                                                               >  2:829567/1‑148 (MQ=255)
                                   tCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTc                                                                                                              <  2:1273260/149‑1 (MQ=255)
                                   tCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGATCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTc                                                                                                              >  1:763437/1‑149 (MQ=255)
                                      tATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaa                                                                                                           <  2:219872/149‑1 (MQ=255)
                                        ttCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAATTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                                                                         <  2:1378273/149‑1 (MQ=255)
                                           cGGCGATTGTCATTCTCTTTGGTATATCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGTCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTACAAATGGc                                                                                                      <  1:605415/149‑1 (MQ=255)
                                           cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCATATGGc                                                                                                      >  2:1476329/1‑149 (MQ=255)
                                           cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                                                                      >  2:732322/1‑149 (MQ=255)
                                            ggCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                                                                     <  2:538304/149‑1 (MQ=255)
                                             gCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTGTTTTACGTGTTACCTAAATTTGTCTATATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTAAGGAAATCACTTATGTCCAAAtttcaa                                                                                                    >  2:1037033/1‑144 (MQ=255)
                                              cGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAAc                                                                                                   <  1:1329011/149‑1 (MQ=255)
                                                aTTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                                 <  1:1512387/149‑1 (MQ=255)
                                                aTTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                                 >  1:823916/1‑149 (MQ=255)
                                                       ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGGTAACTAAATTTGTCTATCTCCTTCAAGAAGAAAGCCATTAGAGTGCGTATTGCCTCTAATGATATCACTGATGTCAAAAGTGCAACTCGCGTGAt                                                                                          >  1:1021474/1‑149 (MQ=255)
                                                          tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                       <  1:341087/149‑1 (MQ=255)
                                                           ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctc                                                                                      >  1:1520820/1‑149 (MQ=255)
                                                            tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                     <  1:350594/149‑1 (MQ=255)
                                                               ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                                                                  <  2:978367/149‑1 (MQ=255)
                                                               ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                                                                  <  1:1166409/149‑1 (MQ=255)
                                                               ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCGGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                                                                  <  1:1172398/149‑1 (MQ=255)
                                                                gTATACCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTAGCTAAATTTGGCTAAATCCTTCTAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAGTGTTGTCCATATGGCAACTCGCCGGATcctactgga                                                                                 >  2:895322/1‑146 (MQ=255)
                                                                gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 >  1:426128/1‑149 (MQ=255)
                                                                gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 <  1:1241358/149‑1 (MQ=255)
                                                                gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAAATTGGCTAAATCCTTCAAGTAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 <  1:467515/149‑1 (MQ=255)
                                                                gTATAACCGTGATTTTGTGCGGGTTCGGCAACGTCGTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGGCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                 <  2:1115031/149‑1 (MQ=255)
                                                                  ataAGCGTGATTTTGTGCTGGTTCGGCATCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                                                               <  2:36494/149‑1 (MQ=255)
                                                                  ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                                                               >  2:1138191/1‑149 (MQ=255)
                                                                  ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                                                               >  1:1365395/1‑149 (MQ=255)
                                                                       ggTGATTTTGTGCTGGTTCTGCCACGTCTTTAACGTGTTACATTATTTTGTCTAAATCCTGCAAGAAGCCAGCCATTAGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCAGATCCTCCTTCACCACGTa                                                                          <  1:1057598/148‑1 (MQ=255)
                                                                             tttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCttt                                                                    <  1:1138191/149‑1 (MQ=255)
                                                                                 tgCTGGTTCGGCAACGTCTTTAAAGTGTAACCTAAATTTTGCTAAATCCTTCCTGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCCCGTATGCTTTGCGt                                                                <  2:286980/149‑1 (MQ=255)
                                                                                    tGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc                                                             >  2:449358/1‑149 (MQ=255)
                                                                                      gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCCGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGGCCAAATGGCAACTCTCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                           >  2:1389209/1‑149 (MQ=255)
                                                                                      gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                           >  1:1141184/1‑149 (MQ=255)
                                                                                      gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                           >  2:1472469/1‑149 (MQ=255)
                                                                                           gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa                                                      <  2:532020/149‑1 (MQ=255)
                                                                                           gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa                                                      >  1:724765/1‑149 (MQ=255)
                                                                                           gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa                                                      >  2:1240488/1‑149 (MQ=255)
                                                                                           gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCATACTa                                                      >  2:383308/1‑149 (MQ=255)
                                                                                           gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTGCGTGCCTCTCAGGAAAACACTAATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa                                                      >  1:1280814/1‑149 (MQ=255)
                                                                                             aaagtatttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGATGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                    <  2:935822/143‑1 (MQ=255)
                                                                                             aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCCACTCGCATGATCCACCGTCACACCGTATGCTTTGCGGCACCTTACTTTc                                                    >  2:189103/1‑149 (MQ=255)
                                                                                                gTCTGTAACGTGTTACCTAAAATTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCAGGCTTATTGCATCACAGGAAATAACTTATTACCAAATGGCAACTCGCCTGATCCTCCTTAACCACGACTGCTTGTCTTCACCTTAAAATCAgc                                                 >  2:1333867/1‑148 (MQ=255)
                                                                                                 tCTTTATCGTGGTACCTAAATTTGGCTAACTCCTTAAAGAAGCCAGCCATTCTCTGGCTTCTTGCCTCTCAGTAAATCACTTATGTCCAATTTGAAACTCGCCTGATCCCCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGa                                                <  2:645421/149‑1 (MQ=255)
                                                                                                 tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGa                                                >  1:459699/1‑149 (MQ=255)
                                                                                                 tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCGTTGCGTCACCTTACTATCAGGa                                                >  2:272748/1‑149 (MQ=255)
                                                                                                    ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGACGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGc                                             <  2:127087/149‑1 (MQ=255)
                                                                                                     tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATTTCCAAATGGAAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCt                                            >  1:316866/1‑149 (MQ=255)
                                                                                                      aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCtt                                           <  2:786069/149‑1 (MQ=255)
                                                                                                      aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGACTCTCAGGAAATCAGTTATGTCCAAAGGGCAACTCGCCTGATCCTCCTTGACCACGAATGGTTTGCGTCACCTTACTATCAGGTCGGtt                                           >  1:317880/1‑149 (MQ=255)
                                                                                                         gtgtTACCTAAATTTTGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                        <  1:196184/149‑1 (MQ=255)
                                                                                                         gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCCCCTTACTATCAGGACGCTTTAg                                        >  2:639896/1‑149 (MQ=255)
                                                                                                         gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                        >  1:1561007/1‑149 (MQ=255)
                                                                                                           gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATGCTCCTTCACCACGTATGCTTTGCGTCACCTTACTAGCAGGACGCTTTAGcc                                      >  1:537289/1‑149 (MQ=255)
                                                                                                            ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGccc                                     <  1:933566/149‑1 (MQ=255)
                                                                                                              aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCTCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATAAGGACGCTTTAGCCCAt                                   >  1:1314759/1‑149 (MQ=255)
                                                                                                               ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAATTCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCagg                                  >  1:430030/1‑147 (MQ=255)
                                                                                                               ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCTTg                                  >  1:670533/1‑149 (MQ=255)
                                                                                                               ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATg                                  <  1:711936/149‑1 (MQ=255)
                                                                                                                 tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCCCGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTc                                >  1:677819/1‑149 (MQ=255)
                                                                                                                 tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTc                                <  1:1427402/149‑1 (MQ=255)
                                                                                                                      ttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCt                           >  1:99423/1‑149 (MQ=255)
                                                                                                                         gCTAAATCCTTCAAGAAGCCAGCCATTCTCTGGCTTCTTGCCTCTCAGGAAATCCCTTATGTCCAACTGGCAACACGCCTGATCCTCCTTCACCACGTGTGCTTTGCGTCCGCGTACTGTCTGGCGGCTTGAGCACATGTCCCGCtttt                        >  2:1027818/1‑149 (MQ=255)
                                                                                                                         gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCACCTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtgtt                        >  2:915209/1‑149 (MQ=255)
                                                                                                                          cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACGATCAGGTCGGTTTAGCCCATGTCCCGCttttt                       >  2:1550819/1‑149 (MQ=255)
                                                                                                                           tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                      >  2:749831/1‑149 (MQ=255)
                                                                                                                                     atgaagCCAGCATTTCGCTGGATTCTTGGATCTCATTTAATTACTTATTTCCAATTGGCATCTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCTCCTTACTATCATGACGCTTTAGCCCATGTCCCGCTTTTTGATTTgttgtt            <  1:471657/147‑1 (MQ=255)
                                                                                                                                     aagaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTAAACACGTATGCTTTGGGTCACCTTACTATCAGGACGCTTTAGCCCATTACCCCCTTTTTGATTTGTAggt            >  2:1303537/1‑147 (MQ=255)
                                                                                                                                       gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt          >  2:1220693/1‑149 (MQ=255)
                                                                                                                                        aagACAGCGAATGGCTGGCTGCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTTATCCTCCTTCACCACGTATGCATTGCGGCACCTTACAATCAGGCCGATTTAGACCCTGACCCGCTTTTTGAATTGTAGTTTTg         >  2:1560791/1‑149 (MQ=255)
                                                                                                                                          gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGcc       <  2:220634/149‑1 (MQ=255)
                                                                                                                                             agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTg    <  2:559658/149‑1 (MQ=255)
                                                                                                                                              gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATAACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTGCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGACCCGCTTTTTGATTTGTAGTTTTGCACTgg   >  1:555379/1‑149 (MQ=255)
                                                                                                                                               ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGt  <  1:1135723/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGT  >  NC_000913/3891358‑3891649

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: