Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,182,877 T→G F545V (TTC→GTC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,8770TG98.5% 230.3 / NA 66F545V (TTC→GTC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (33/32);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAA  >  NC_000913/4182735‑4183023
                                                                                                                                              |                                                                                                                                                  
gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGtt                                                                                                                                            <  1:6931/151‑1 (MQ=255)
gATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGtt                                                                                                                                            >  1:93341/1‑151 (MQ=255)
   ttCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCg                                                                                                                                          >  2:583587/1‑150 (MQ=255)
     ccGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGtt                                                                                                                                            >  2:323982/1‑146 (MQ=255)
            agTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTaca                                                                                                                                 <  1:169637/150‑1 (MQ=255)
             gTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACAcc                                                                                                                               <  1:266286/151‑1 (MQ=255)
                 aaGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGAc                                                                                                                           >  2:780127/1‑151 (MQ=255)
                   gagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTc                                                                                                                         >  1:116841/1‑151 (MQ=255)
                    agTTCTTCGGTTCCAGCCAGTTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTc                                                                                                                         >  1:566472/1‑150 (MQ=255)
                     gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGCCCAGAACACCCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCa                                                                                                                        <  1:323558/150‑1 (MQ=255)
                      ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACt                                                                                                                      >  2:611628/1‑151 (MQ=255)
                         ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACg                                                                                                                   >  1:761458/1‑151 (MQ=255)
                          tcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTa                                                                                                                     <  2:810127/148‑1 (MQ=255)
                            ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTc                                                                                                                >  1:356204/1‑151 (MQ=255)
                             gTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcg                                                                                                               <  2:162604/151‑1 (MQ=255)
                                  agccggcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTGCGAGACGTACCCCCGACTCACTACGGTCGCGTAt                                                                                                          <  1:48906/151‑1 (MQ=255)
                                   gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAt                                                                                                          >  2:253249/1‑150 (MQ=255)
                                    ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGt                                                                                                        >  2:471119/1‑151 (MQ=255)
                                        cTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCa                                                                                                     >  1:593267/1‑150 (MQ=255)
                                         tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa                                                                                                    >  2:358959/1‑150 (MQ=255)
                                         tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt                                                                                                   >  2:740588/1‑151 (MQ=255)
                                              cAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg                                                                                                 <  2:279790/148‑1 (MQ=255)
                                               aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGa                                                                                                <  2:700459/148‑1 (MQ=255)
                                               aGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAc                                                                                             >  1:437316/1‑151 (MQ=255)
                                                 tttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAccc                                                                                           >  2:604614/1‑151 (MQ=255)
                                                   tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                         >  2:179610/1‑151 (MQ=255)
                                                    aTGGACCAGAACAACCCGCTGTCGGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                         <  1:552638/150‑1 (MQ=255)
                                                      ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                                                                                      <  1:404288/151‑1 (MQ=255)
                                                       gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                                                                       >  2:8351/1‑149 (MQ=255)
                                                         ccAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt                                                                                   >  2:804630/1‑151 (MQ=255)
                                                               caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                                                                             <  1:365148/151‑1 (MQ=255)
                                                                aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACa                                                                            <  1:455053/151‑1 (MQ=255)
                                                                   ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCg                                                                         >  1:537283/1‑151 (MQ=255)
                                                                       tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                                                                     >  1:70535/1‑151 (MQ=255)
                                                                       tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                                                                     >  1:208174/1‑151 (MQ=255)
                                                                        gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                                                                     <  1:81033/150‑1 (MQ=255)
                                                                                   cGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                         <  1:893300/151‑1 (MQ=255)
                                                                                    gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc                                                          >  2:550483/1‑149 (MQ=255)
                                                                                     cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGt                                                       >  2:469744/1‑151 (MQ=255)
                                                                                       caAACGTCGTATCTCCGCACTCTTCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTGGAAGATCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAGATCTGTCTGATCAACTCTCTGTc                                                      >  1:2550/1‑150 (MQ=255)
                                                                                       caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                      <  1:314072/150‑1 (MQ=255)
                                                                                           cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                 <  1:211708/151‑1 (MQ=255)
                                                                                            gtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                 >  1:215876/1‑150 (MQ=255)
                                                                                             tcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtg                                                  >  1:563215/1‑148 (MQ=255)
                                                                                              cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                              <  2:133546/151‑1 (MQ=255)
                                                                                                 aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                            <  1:728649/150‑1 (MQ=255)
                                                                                                   ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                           >  1:317609/1‑149 (MQ=255)
                                                                                                    tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                           >  2:461981/1‑148 (MQ=255)
                                                                                                     ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                           <  2:781205/147‑1 (MQ=255)
                                                                                                      cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                      <  1:546172/151‑1 (MQ=255)
                                                                                                            cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                  <  1:312365/149‑1 (MQ=255)
                                                                                                             ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                               >  2:449004/1‑151 (MQ=255)
                                                                                                             ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                               <  2:752913/151‑1 (MQ=255)
                                                                                                               cccAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                             >  2:47622/1‑151 (MQ=255)
                                                                                                                ccAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                            >  1:116281/1‑151 (MQ=255)
                                                                                                                 cAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                           >  1:745380/1‑151 (MQ=255)
                                                                                                                       gTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                     <  1:703075/151‑1 (MQ=255)
                                                                                                                          tGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                    >  1:837230/1‑149 (MQ=255)
                                                                                                                           gACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                    <  2:844040/148‑1 (MQ=255)
                                                                                                                             cccGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc               <  2:166936/151‑1 (MQ=255)
                                                                                                                                  gAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGt          <  2:130560/151‑1 (MQ=255)
                                                                                                                                   aaCGTGCAGGCGTCGAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa         <  1:631733/151‑1 (MQ=255)
                                                                                                                                    aCGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTAt        <  2:748693/151‑1 (MQ=255)
                                                                                                                                    aCGTGCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTAt        <  1:628966/151‑1 (MQ=255)
                                                                                                                                        gCAGGCGTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa    <  2:240577/151‑1 (MQ=255)
                                                                                                                                          aGGCGTCGAAGTTCGAGACGTACACCCGAATCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa  <  2:561833/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAA  >  NC_000913/4182735‑4183023

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