Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,502 C→G 5.1% intergenic (+19/‑113) 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,5020CG 5.1% 214.7 / 2.5 77intergenic (+19/‑113)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (41/32);  new base G (2/2);  total (43/34)
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.42e-01

AAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTT  >  NC_000913/3891355‑3891650
                                                                                                                                                   |                                                                                                                                                    
aagTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTAATCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAgccagcca                                                                                                                                                     <  2:40419/149‑1 (MQ=255)
   tttCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTc                                                                                                                                                  >  1:1004793/1‑149 (MQ=255)
       ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGTTTGTCATTCTCTTTCGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                              <  2:593584/149‑1 (MQ=255)
        cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                             <  1:172308/149‑1 (MQ=255)
          aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                                                                           >  2:392362/1‑149 (MQ=255)
           tCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGctt                                                                                                                                          >  1:493983/1‑149 (MQ=255)
             atttGTTCACGGTCTGGGGCTGCAATCTTATTCCGGCGATTGTCATCCTCTTTGGTAAAACCGTGATTATGTGATGGTTCGGCTACGTCTTTAACGTGATACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcgtct                                                                                                                                        <  1:765163/146‑1 (MQ=255)
                tGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                                                                     <  1:328135/149‑2 (MQ=255)
                tGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCAATCGCTGGCTTCTTGc                                                                                                                                     <  2:713908/149‑1 (MQ=255)
                 gTTCAGGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTTGTATAACCGTGATGTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                                                                    <  1:64249/149‑3 (MQ=255)
                      cGGTCTGGGGCGGCATTCTTATTCCGGCGATTGTCATTCTATTTGGTATAACCGTTATTTTGTGCTGGTTCGGCAACGTGTTTAACGTGGTACCTAAATTTGGCTAAATCCTTAAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                                                                               <  1:705426/149‑1 (MQ=255)
                      cGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                                                                               >  2:559315/1‑149 (MQ=255)
                              ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTTGTATTACCGGTATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACGTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                       <  1:770071/149‑1 (MQ=255)
                                cggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                     >  2:72949/1‑149 (MQ=255)
                                    aaTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAATTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                                                                 <  1:697180/149‑1 (MQ=255)
                                    aaTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                                                                 <  1:426959/149‑1 (MQ=255)
                                     aTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTTGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGt                                                                                                                <  2:567126/149‑1 (MQ=255)
                                          aTTCCGGCGATAGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCAAAATTTGGCTAAATCATTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                                                                           >  1:862628/1‑149 (MQ=255)
                                           ttCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCTCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                                                                          >  2:501422/1‑149 (MQ=255)
                                           ttCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGTAAATCACTTATGTCCAAAt                                                                                                          >  2:1184567/1‑149 (MQ=255)
                                           ttCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTAGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                                                                          >  2:546196/1‑149 (MQ=255)
                                             ccGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                        >  2:994081/1‑149 (MQ=255)
                                              cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTACGTGTTAACTAAATTTGTCTAAATACTTCAAGAAGCCAGCCATTCCCAGGCTTCTTGCCTCGCAGGAAATCACGTATGGACaaaatgca                                                                                                      >  2:882805/1‑144 (MQ=255)
                                              cGGCGATTGTCATTCTCTTTGGTATAACCCTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCCGCCATTCGCAGGCTTCTGTCCTCTCAGGATATCACTTATGTCAAAATGGc                                                                                                       >  2:672120/1‑149 (MQ=255)
                                              cGGCGATTGTCATTCTCTTTGGTAAAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                                                                       <  2:1094762/149‑1 (MQ=255)
                                                gCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACGTAAATTTGGCTAAATCCTTCATGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCaa                                                                                                     <  2:415029/149‑1 (MQ=255)
                                                gCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCaa                                                                                                     <  1:390261/149‑1 (MQ=255)
                                                   aTTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                                  <  1:1134486/149‑1 (MQ=255)
                                                   aTTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                                                                  <  1:49755/149‑1 (MQ=255)
                                                       tCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTTGATAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATTGCAACtatcct                                                                                              >  1:790359/1‑144 (MQ=255)
                                                         aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTATATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGATGGCTTCTTGCATCTCAGGAAATAACTTATGTCCAAATGGATACTCGCCTGa                                                                                            >  1:911492/1‑149 (MQ=255)
                                                         aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                            >  1:1059454/1‑149 (MQ=255)
                                                         aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                            >  2:1067222/1‑149 (MQ=255)
                                                         aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGACTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                            >  1:228468/1‑149 (MQ=255)
                                                         aTTCTCTTTGGTATAACCGTGATTTTGTCCTCGTTAGGCAACGTCTTTAACGTGTTACCTAAATTTGTCTAAATCCTTCAAAAAGCCATCCATTCACTTGCTTCTTGCCTATCTGGAAATCAATTATTGCCAAATTGAAAATCGCCTGa                                                                                            >  1:789250/1‑149 (MQ=255)
                                                               tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGCtcctcc                                                                                      >  1:1118219/1‑149 (MQ=255)
                                                               tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGGCTTTCACGTGTTACCTAAATTTGGCTAAATCCTTCAAGACGCCAGCCATTCGCTGGCTTCTTGCCTCTCACGAAAGCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                      <  1:1052184/149‑1 (MQ=255)
                                                                   gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                                                                  >  1:516235/1‑149 (MQ=255)
                                                                     ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCATCTCGCCTGATCCTCCTTcacc                                                                                >  2:45969/1‑149 (MQ=255)
                                                                     ataACCGTGATTTTGTGCTGCTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCATTAAAGAAGCCAGCCATTACCTGGATTCTTGCCTCTCAGGAACTCACTTATGTCCAAATGGCAACTCGCCAGATCCAACTTcacc                                                                                >  2:848015/1‑149 (MQ=255)
                                                                        aCCGTGATTTTGTGCTGGTTCGGGAACGACTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACg                                                                             <  1:85944/149‑1 (MQ=255)
                                                                          cGTGATTTTGTGCTGGTTCGGTAACGTATTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                                                                           <  1:674570/149‑1 (MQ=255)
                                                                          cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCGTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcaccacttt                                                                           >  1:89574/1‑146 (MQ=255)
                                                                           gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                                                                          >  1:390331/1‑149 (MQ=255)
                                                                                  tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                   <  1:1134324/149‑1 (MQ=255)
                                                                                    tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGt                                                                 >  1:692605/1‑149 (MQ=255)
                                                                                    tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTGTTGTGCAAATGGCAACTCGCCTGATCCGCGTTCCCCACGTATGCTTTTCGt                                                                 >  1:324472/1‑149 (MQ=255)
                                                                                     gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                                                                >  1:143001/1‑149 (MQ=255)
                                                                                      cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                                                               >  1:1142251/1‑149 (MQ=255)
                                                                                      cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGAATCCTTTGCgacc                                                               >  2:64056/1‑146 (MQ=255)
                                                                                        ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc                                                             >  1:986071/1‑149 (MQ=255)
                                                                                        ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGATGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc                                                             >  1:1037463/1‑149 (MQ=255)
                                                                                         gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCATTCAAGAAGCCATCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt                                                            <  1:71033/149‑1 (MQ=255)
                                                                                          ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCtt                                                           <  1:1137936/149‑1 (MQ=255)
                                                                                               cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCGTCTTGCCTCTCAGGAAATCAATTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTAt                                                      <  2:77772/149‑1 (MQ=255)
                                                                                                aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc                                                     <  1:588645/149‑1 (MQ=255)
                                                                                                 aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGAGTCACCTTACTATcc                                                    >  1:201924/1‑148 (MQ=255)
                                                                                                  cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg                                                   >  1:1110872/1‑149 (MQ=255)
                                                                                                      tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACg                                               >  1:729694/1‑149 (MQ=255)
                                                                                                       ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGc                                              <  2:692605/149‑1 (MQ=255)
                                                                                                         aaCGTGTTACCTAAATTAGGCTAAATCCTTCAAGAAGCCAGCCAATCGCTGGCTTCTTGACTCTCAGGAAATCACTTATGTCCAAATGGAAAAACGCCTGATCCTACTACACCACGTATGCATTGCGTCACCTTACTATAAGTACGCtt                                            >  1:258822/1‑149 (MQ=255)
                                                                                                            gtgtTACCTAAATTTGGCTAAATCCTTCCAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCTAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                         <  2:257111/149‑1 (MQ=255)
                                                                                                            gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                         <  1:1067222/149‑1 (MQ=255)
                                                                                                               ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGccc                                      <  2:261025/149‑1 (MQ=255)
                                                                                                               ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGccc                                      <  1:546196/149‑1 (MQ=255)
                                                                                                                    tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCCCGTATGTCCAAATGGCAACTCGCCTGAGGCTCCTTCACCACGGATGCTTTGCGTCACCTTACTATCAGGACGGGTTAGCCCATGTc                                 >  1:813207/1‑149 (MQ=255)
                                                                                                                    tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATACTCCTTCACCGCGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGtg                                 >  1:253016/1‑148 (MQ=255)
                                                                                                                       aTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCg                              >  1:864636/1‑149 (MQ=255)
                                                                                                                             cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGAAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTAGTATCAGTACGCTTTAGCCCATGTCCCGCttttt                        <  2:493983/149‑1 (MQ=255)
                                                                                                                             cTAAATCCTTCAAGAAGCAAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTATGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                        <  2:569963/149‑1 (MQ=255)
                                                                                                                               aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGa                      <  1:72949/149‑1 (MQ=255)
                                                                                                                                aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAACTGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAt                     <  1:338814/149‑1 (MQ=255)
                                                                                                                                aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAt                     <  2:747162/149‑1 (MQ=255)
                                                                                                                                     ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCACATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGtt                >  1:721052/1‑148 (MQ=255)
                                                                                                                                      gcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGCTCCTCCTTCACCACGTATGCTTTGCGGCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAg               >  2:234663/2‑149 (MQ=255)
                                                                                                                                       cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTAAGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGt              >  2:238310/1‑149 (MQ=255)
                                                                                                                                         agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGttt            >  1:280249/1‑149 (MQ=255)
                                                                                                                                          gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt           >  2:3089/1‑149 (MQ=255)
                                                                                                                                             gccagcgaATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGcc        >  1:161204/1‑149 (MQ=255)
                                                                                                                                              ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCAGATCCTCCTTCACCACGTATGCTTTGCGTCACCATACTATCAGGACGCTTTAGCCCATGTACCGCTTTTTGATTTGTATTTTTGccc       >  2:316602/1‑149 (MQ=255)
                                                                                                                                              ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAAACACTTATGTCCAAATGGCAACTCGCCTGAACCTCCTTCACAACGTATGCTTTGCGTCACCTTACTATCAGCACGCTTTATACCATTTCCAGCTTATTGATTTGTAGTTTTGcca       >  2:621125/1‑148 (MQ=255)
                                                                                                                                                agccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTg     >  1:209812/1‑149 (MQ=255)
                                                                                                                                                 gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACGTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTgg    >  2:760274/1‑149 (MQ=255)
                                                                                                                                                 gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTAAGCCCATGGCACGCTTTTTGATTTGTAGTTTTTCCCTgg    >  1:703691/1‑149 (MQ=255)
                                                                                                                                                  ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTCATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGt   <  2:534021/149‑1 (MQ=255)
                                                                                                                                                  ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGAAACTCGCCTGATCCTCCGTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGt   >  1:1102668/1‑149 (MQ=255)
                                                                                                                                                   caTTCGCTGGCTTGTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCTCCTCGTTTGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtt  <  1:672120/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTT  >  NC_000913/3891355‑3891650

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: