Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,182,928 A→T 12.2% E562V (GAA→GTA)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,182,9280AT12.2% 246.0 / 23.7 97E562V (GAA→GTA) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (50/34);  new base T (5/7);  total (56/41)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

Only 100 of 104 total aligned reads displayed.

GTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAA  >  NC_000913/4182782‑4183076
                                                                                                                                                  |                                                                                                                                                    
gTTTATGGCCCAGAACACCCCGCTGTCTGGGATTACGCACAAACGTCGTACCTCCGCACTCGGCCCAGGCGGTCTGACCCGTCAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAATACGGTCGCGTATGTCCAATCGaaa                                                                                                                                                    <  1:728816/149‑1 (MQ=255)
gTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGATTCGAAGGTCGAGACGTACACCCGACTCACTACGGTCCCGTATGTCCAATCGaaa                                                                                                                                                    >  1:1075403/1‑149 (MQ=255)
gTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGGGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGCCGCGTATGTCCAATAGaaa                                                                                                                                                    >  2:721001/1‑149 (MQ=255)
       gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                                                                                                                             >  1:796154/1‑149 (MQ=255)
       gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                                                                                                                             >  2:66111/1‑149 (MQ=255)
        aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                                                                                                                                            <  1:1389241/149‑1 (MQ=255)
          cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGCAACCCCTGAAgg                                                                                                                                          >  2:60325/1‑149 (MQ=255)
            gAACAACCCGCTGTCGGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTc                                                                                                                                        <  1:1062697/149‑1 (MQ=255)
             aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                                                                                                                                       <  1:1307834/149‑1 (MQ=255)
                aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACGCCTGAAGGTCCgtg                                                                                                                                    >  2:1084842/1‑147 (MQ=255)
                aacCCGCTCTCTGAGATTACGCACAAACGTCGTATCTCCGCGCACGGCCCAGGCGGTCTGTCCCGTGAACATGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGaa                                                                                                                                    <  2:70/149‑1 (MQ=255)
                   ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGCCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACAt                                                                                                                                 >  1:181971/1‑149 (MQ=255)
                      cTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCACTGAAGGGCCGAACATCgg                                                                                                                              >  2:65937/1‑149 (MQ=255)
                           tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAAGTCCAATCGAAACCCCTGAAGGTCCGAGCATCGGTCTGa                                                                                                                         >  1:1032655/1‑149 (MQ=255)
                           tGAGATTACGCACAAACGTCGTATCTACGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTTTGa                                                                                                                         >  2:394196/1‑149 (MQ=255)
                            gggaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAt                                                                                                                        <  2:1066214/147‑1 (MQ=255)
                            gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGCTCGCGTATGTCCAATCGTAACCCCTGCAGGTCCGCACCTCGGTCCGac                                                                                                                        >  2:888154/1‑148 (MQ=255)
                             agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATc                                                                                                                       <  1:1220971/149‑1 (MQ=255)
                                ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCGTCGAAGTTCGAGACGTACACCCGCCTACAGACGATCGCGCAAGTCCACCCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAAc                                                                                                                    >  2:524253/1‑149 (MQ=255)
                                  aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTAGCGTATGTCCAATCGAAACCCCGGAAGGTCCGAACATCGGTCTGATCAActa                                                                                                                  >  2:416869/1‑148 (MQ=255)
                                   cGCACAAACGTCGTATCCCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActct                                                                                                                 <  1:763590/149‑1 (MQ=255)
                                    gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCGTGAAGGTCCGAACATCGGTCTGATCAActctc                                                                                                                >  2:817144/1‑149 (MQ=255)
                                    gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc                                                                                                                >  2:609533/1‑149 (MQ=255)
                                    gCACAAACGTCCTATCTCCGCACTAGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGACACCCCTGAAGGTCCGAACATCGGTCTGATCCActcac                                                                                                                >  2:1105539/1‑147 (MQ=255)
                                     cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                                                                               <  2:1303056/149‑1 (MQ=255)
                                      acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGACACGTACACCCGACTCACTACGGTCGCGTATGTCCACTCGAAACCCCTGAAGGTCCGAACATCGGTCGGATCAACCCTCTg                                                                                                              >  2:1199107/1‑149 (MQ=255)
                                       caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGt                                                                                                             <  1:1312923/149‑1 (MQ=255)
                                        aaccgccgtATCCCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGGGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                                                                            <  1:1249682/143‑1 (MQ=255)
                                         aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc                                                                                                           <  2:1357188/149‑1 (MQ=255)
                                         aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCAAAACAGCGGTCTGATCAACTCCCTGTcc                                                                                                           >  2:235091/1‑149 (MQ=255)
                                           cgtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGAGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt                                                                                                         <  1:61450/149‑1 (MQ=255)
                                            gtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtg                                                                                                        >  1:983678/1‑149 (MQ=255)
                                              cgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTa                                                                                                      <  2:492661/149‑1 (MQ=255)
                                                   ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGACCGCGCTGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                                 <  1:1123863/149‑1 (MQ=255)
                                                   ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                                 >  1:462423/1‑149 (MQ=255)
                                                   ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                                                                 >  1:826692/1‑149 (MQ=255)
                                                    tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcgca                                                                                                >  1:185117/1‑149 (MQ=255)
                                                        cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCCCAGACt                                                                                            >  2:1227246/1‑149 (MQ=255)
                                                        cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                                                                            >  1:518554/1‑149 (MQ=255)
                                                        cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTAAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                                                                            <  1:999692/149‑1 (MQ=255)
                                                           tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                         >  1:999347/1‑149 (MQ=255)
                                                           tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                         <  1:948435/149‑1 (MQ=255)
                                                            cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                                                                        >  2:1280012/1‑149 (MQ=255)
                                                                 cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTCTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGTCCGCGTACGCACAGACTAACGAATAc                                                                                   >  2:867659/1‑149 (MQ=255)
                                                                 cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTCACGAATAc                                                                                   >  1:908507/1‑149 (MQ=255)
                                                                 cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCCTGTACCCACAGACTAACGGATAc                                                                                   >  2:1229382/1‑149 (MQ=255)
                                                                      ggTCTGACCCGTGAACGTGCAGGCTGCGAAGTTCGCGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGACCGCGTACGCCAAGACTAACGAAAACGGCtt                                                                              >  2:380908/1‑149 (MQ=255)
                                                                       gTCTGACCCGGGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                                                                             <  1:435110/149‑1 (MQ=255)
                                                                        tCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                                                                            <  1:1408313/149‑1 (MQ=255)
                                                                          tGACCCGTGACCGTGCAGGCTTCGATGTTCGACACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                          <  2:1421334/149‑1 (MQ=255)
                                                                          tGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCCCTACGGTCGCGTATGTCCACTCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCCCTGTCCGTGGACGCACGGACTAACGACTACGGCTTCCtt                                                                          >  1:491753/1‑149 (MQ=255)
                                                                           gACCCGTGACCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                                                                         <  2:772709/149‑1 (MQ=255)
                                                                            aaaacaaaaCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                                                                        >  2:968157/8‑149 (MQ=255)
                                                                            aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAGCATCGGTCTGATCAACTCTCTGTCCGTGGTCGCACAGACTAACGAATACGGCTTCCTTga                                                                        <  2:1035135/149‑1 (MQ=255)
                                                                            aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                                                                        <  1:1095178/149‑1 (MQ=255)
                                                                              ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                                                                      >  2:1406617/1‑149 (MQ=255)
                                                                               cGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCGCAGACTAACGAATACGGCTTCCTTGAGAc                                                                     >  1:1430557/1‑149 (MQ=255)
                                                                                   aaCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                                                                 >  1:526761/1‑149 (MQ=255)
                                                                                       tGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGATCACCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATc                                                             <  2:859115/149‑1 (MQ=255)
                                                                                        gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg                                                            >  2:1131662/1‑149 (MQ=255)
                                                                                        gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTACCg                                                            >  1:324923/1‑149 (MQ=255)
                                                                                          aGGCTTCGAAGTTCGGGACGTACACCCGACTCACTACGGTCGCGTATGTCCATTCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa                                                          <  2:1054279/149‑1 (MQ=255)
                                                                                            gCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                                        >  2:1191431/1‑149 (MQ=255)
                                                                                             cTTCGAAGTTCGAGACGTACAGCAGGCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCGCTGTCCGTGTACGCACAGCCTAACGAATACGGCTTCCTTGAGACGCCGTATCGTGAAg                                                       >  2:530146/1‑149 (MQ=255)
                                                                                                cGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTCCGCTTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa                                                    >  1:1426316/1‑149 (MQ=255)
                                                                                                    gTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCACCGACCGGCCGGACAGCGGGCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTCTCGTAAAGTGACAGa                                                >  2:528778/1‑149 (MQ=255)
                                                                                                     ttCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGAc                                               >  2:212575/1‑149 (MQ=255)
                                                                                                      tCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACg                                              <  1:1235616/149‑1 (MQ=255)
                                                                                                       cGAGCCGTACACCGCCCTCACTTCGGTGGCGTATGTCCAATCGAACCCGCTTAAGGTCCGAACATCTGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATCCGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACgg                                             <  2:287056/149‑1 (MQ=255)
                                                                                                         agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGgtg                                           >  2:1348388/1‑149 (MQ=255)
                                                                                                           ccgcacacCCGACTCACTCCGGTCGCGTATGTCCAATCGAACCCCCTGAAGGTTCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACTCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTTAAGTGACCGACGGtgtt                                         <  2:1075403/145‑1 (MQ=255)
                                                                                                              tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACCGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGGCGGTGGTGTa                                      >  2:458034/1‑149 (MQ=255)
                                                                                                               acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAGCCCCTGAAGGACCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTaa                                     <  1:112781/149‑1 (MQ=255)
                                                                                                               acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTaa                                     <  2:826692/149‑1 (MQ=255)
                                                                                                                  cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCCTGAGACCCCGAATCGTAAAGTGACCGACGGTGTTGTAACTg                                  >  1:401135/1‑149 (MQ=255)
                                                                                                                   ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATGAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATAGTCACGTGACAGACGGTGTTGTAACGGa                                 >  1:555370/1‑149 (MQ=255)
                                                                                                                   ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGa                                 >  1:1265336/1‑149 (MQ=255)
                                                                                                                        tCACTACGGTCGCGTATGTCCAATCGTAAACCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGaaa                            >  1:890101/1‑149 (MQ=255)
                                                                                                                          aatACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTACCGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                          <  1:1313494/147‑1 (MQ=255)
                                                                                                                          aCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAAtt                          <  1:1162971/149‑1 (MQ=255)
                                                                                                                          aCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGGCACTGACggcctt                          >  1:1220979/1‑144 (MQ=255)
                                                                                                                            tACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACGGGCGGTGTTGTAACTGACGAAATTCa                        >  1:240241/1‑149 (MQ=255)
                                                                                                                            tACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCa                        <  1:524253/149‑1 (MQ=255)
                                                                                                                             tcggctgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCAc                       <  2:1015263/143‑1 (MQ=255)
                                                                                                                             gcgggcgggTACGTCTAATCGTGACCCCTTAAGGTCCGCACATCGGTGTGATCATCTCTCTTTCCGTGTACGCACAGACTAACGAATACGGCGTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCAc                       <  2:1136798/148‑1 (MQ=255)
                                                                                                                               ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTGACTGACGAAATCCACtc                     >  2:100110/1‑148 (MQ=255)
                                                                                                                               ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTa                     <  1:609533/149‑1 (MQ=255)
                                                                                                                                gTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAGCTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAc                    <  2:1415618/149‑1 (MQ=255)
                                                                                                                                 tCGCGTACGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAcc                   >  2:262190/1‑149 (MQ=255)
                                                                                                                                  cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCt                  <  2:135120/149‑1 (MQ=255)
                                                                                                                                     gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCGTCCTTGAGACTCCGAATCGTAAAGTGACCGACGGTGTTGAAACTGACGAAATCCACCACCTgcc               >  1:967764/1‑147 (MQ=255)
                                                                                                                                      tATGTCCAATCGTAACCCCTGAAGGTCCGAACACCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCt              <  1:421819/149‑1 (MQ=255)
                                                                                                                                       aggTCCAATCGAAACCCCTGAAGGTCCGAACATCCGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTg             <  2:1141402/147‑1 (MQ=255)
                                                                                                                                            cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATc        >  2:452567/1‑149 (MQ=255)
                                                                                                                                              aTCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCga      >  1:1257729/1‑149 (MQ=255)
                                                                                                                                               tCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCgaa     >  1:1127413/1‑149 (MQ=255)
                                                                                                                                               tCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGACACTCCGTATCGTAAAGTGACCGACGGTGTAGTAACTGACGAAAATCAACCCCTGTCTGCTATCgaa     >  2:609478/1‑149 (MQ=255)
                                                                                                                                                  aaaCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCgaagaa  <  1:721001/149‑1 (MQ=255)
                                                                                                                                                  aaaCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCgaagaa  <  1:1040276/149‑1 (MQ=255)
                                                                                                                                                  aaaCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCgaagaa  >  1:1016851/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTATCGAAGAA  >  NC_000913/4182782‑4183076

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