Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,183,227 A→G 6.1% V661V (GTA→GTG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,183,2270AG 6.1% 194.5 / 2.8 66V661V (GTA→GTG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/2);  total (35/31)
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.49e-01

CAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTT  >  NC_000913/4183080‑4183374
                                                                                                                                                   |                                                                                                                                                   
cAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCTCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTAt                                                                                                                                                    <  6:233120/149‑1 (MQ=255)
 aaCTACGTTATCGCCCAGGCGAACTCCAACTTGGATCAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATc                                                                                                                                                   >  3:98097/1‑149 (MQ=255)
      cGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTc                                                                                                                                               <  4:29959/148‑1 (MQ=255)
      cGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTc                                                                                                                                               <  5:138877/148‑1 (MQ=255)
       gTTATCGCCCAGGCGAACTCCAACTTGGATGCAGCAGGCCACTTCGTAGAAGACCTGGTACCTTGCCGTAGCAAAGGCGCATCCCGCGTGTTCAGCCGCGACCAGGTTGCCTACATGGACGTATCCACCCCGCAGGTGGTATCCGTCgg                                                                                                                                             >  7:105150/1‑149 (MQ=255)
       gTTATCGCCCAGGCGAACTCCAAATTGGATGAAGAAGGGCACTTCGTAGAAGACCGGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGAAGCGACCAGGATGACTACATGGACGTATCCACCCAGCAGGTGGAAACCGTCgg                                                                                                                                             >  5:101179/1‑149 (MQ=255)
        ttATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGAGCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGTCTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGt                                                                                                                                            <  3:86676/149‑1 (MQ=255)
             gCCCAGGCGAACTCCTACTTGTATGAAGAAGGCCACTTCGTAGAAGACCTGGTATCTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCATGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTc                                                                                                                                       <  4:63681/149‑1 (MQ=255)
             gCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGt                                                                                                                                        <  2:198055/148‑1 (MQ=255)
             gCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTc                                                                                                                                       >  5:107135/1‑149 (MQ=255)
                 aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTGTCCGTCGGTGCGTCCCt                                                                                                                                    >  5:215589/1‑148 (MQ=255)
                 aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTGTCCGTCGGTGCGTCCCt                                                                                                                                    >  5:219216/1‑148 (MQ=255)
                 aGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCt                                                                                                                                    >  3:38680/1‑148 (MQ=255)
                    cGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTTGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATc                                                                                                                                <  1:35017/149‑1 (MQ=255)
                         tCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGtt                                                                                                                           >  5:138702/1‑149 (MQ=255)
                           cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTcc                                                                                                                         >  3:147900/1‑149 (MQ=255)
                           cAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTcc                                                                                                                         >  3:49793/1‑149 (MQ=255)
                            aaCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCt                                                                                                                        <  4:49793/149‑1 (MQ=255)
                             aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACTTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTg                                                                                                                       <  6:107135/149‑1 (MQ=255)
                             aCTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTg                                                                                                                       <  7:217405/149‑1 (MQ=255)
                                    tGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGCCCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAacac                                                                                                                 >  1:35211/1‑148 (MQ=255)
                                     gaagaaGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGt                                                                                                                                            >  7:177201/1‑120 (MQ=255)
                                     gaagaaGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGt                                                                                                                                            <  8:177201/120‑1 (MQ=255)
                                         aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTGTCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGAc                                                                                                           <  6:215589/149‑1 (MQ=255)
                                         aagGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTGTCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGAc                                                                                                           <  6:219216/149‑1 (MQ=255)
                                             ccACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGcc                                                                                                        >  8:80160/1‑148 (MQ=255)
                                             ccACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGcc                                                                                                        >  4:223325/1‑148 (MQ=255)
                                              cACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGTCGAATCCAGCTTGTTCAGCCGCGACCAGGTAGACTACATGGACGTATCCACCCAGCAGGTGGTCTCAGTCGGTGCTTCCCTGATCCCGTTCCTGGAACACGATGACGCCa                                                                                                       >  4:39555/1‑148 (MQ=255)
                                                   cGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGt                                                                                                  <  1:132088/148‑1 (MQ=255)
                                                         agaCCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGa                                                                                           <  3:78622/149‑1 (MQ=255)
                                                            ccTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATgg                                                                                        <  2:46209/149‑1 (MQ=255)
                                                            ccTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGGTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATgg                                                                                        <  8:176619/149‑1 (MQ=255)
                                                              tGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGt                                                                                      >  7:66642/1‑149 (MQ=255)
                                                                    cTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGTTGGGTGCGAAc                                                                                <  4:38680/149‑1 (MQ=255)
                                                                     ttGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGAGCAGGTTGACTACATGGACGTATCCACCCCGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGACCa                                                                               >  6:177704/1‑149 (MQ=255)
                                                                        ccGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACACCGAAGGTGGTATCAGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCTTGCATTGATGGGTGCGAACATg                                                                             >  7:173336/1‑148 (MQ=255)
                                                                                aaGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAg                                                                    >  8:179078/1‑149 (MQ=255)
                                                                                      aaTCCAGCTAGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCCGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCACGCCGtt                                                              >  4:72305/1‑149 (MQ=255)
                                                                                          cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                                          <  6:138702/149‑1 (MQ=255)
                                                                                          cAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCATGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGa                                                          >  3:54597/1‑149 (MQ=255)
                                                                                           aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                              <  1:86706/112‑1 (MQ=255)
                                                                                           aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCAt                                                                                              >  2:86706/1‑112 (MQ=255)
                                                                                           aGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGAc                                                         >  6:34747/1‑149 (MQ=255)
                                                                                                     gCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCATCGTCAGGCCGTTCCGACTCTGCGCGCt                                               <  8:66642/149‑1 (MQ=255)
                                                                                                       cgcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGa                                             >  1:17451/1‑149 (MQ=255)
                                                                                                       cgcgACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGa                                             <  7:80160/149‑1 (MQ=255)
                                                                                                                     cTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttg                               <  2:35211/149‑1 (MQ=255)
                                                                                                                      tACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg                              >  4:196134/1‑149 (MQ=255)
                                                                                                                       aCATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCtggttgg                              >  8:131240/1‑148 (MQ=255)
                                                                                                                        cATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTa                            >  2:101476/1‑149 (MQ=255)
                                                                                                                         aTGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTAc                           >  7:137806/1‑149 (MQ=255)
                                                                                                                          tGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACt                          >  3:80322/1‑149 (MQ=255)
                                                                                                                            gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg                         <  8:137806/148‑1 (MQ=255)
                                                                                                                            gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg                         <  2:17451/148‑1 (MQ=255)
                                                                                                                            gACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTg                         <  1:101476/148‑1 (MQ=255)
                                                                                                                            gACGTATCCACCCAGCAGGGGTGATCAGTAGGGGCGTCACTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCGGATAAGCCGCTGGTTGGTACTgg                        <  8:173336/149‑1 (MQ=255)
                                                                                                                                 aTCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATg                    >  4:17022/1‑148 (MQ=255)
                                                                                                                                    ctcccTGCAGGTGGTATCCGGCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAAGATGGAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTTGTACTGGTATGGAAc                <  4:54597/147‑1 (MQ=255)
                                                                                                                                    cACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAAc                >  3:143873/1‑149 (MQ=255)
                                                                                                                                     aCCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACg               >  4:105689/1‑149 (MQ=255)
                                                                                                                                           cagGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGt         >  7:108623/1‑149 (MQ=255)
                                                                                                                                              gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                                          >  6:224004/1‑113 (MQ=255)
                                                                                                                                              gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATaa                                          <  5:224004/113‑1 (MQ=255)
                                                                                                                                              gtggtATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGc      >  5:212032/1‑149 (MQ=255)
                                                                                                                                                  tATCCGTGGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCATCGTCAGGCCGTTCCGACTCTGGGCGCTGATAAGCCGCTTGTTGGTACTTGTATGGAACGTGCTGTTGCCGtt  <  4:98097/149‑1 (MQ=255)
                                                                                                                                                  tATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGtt  <  3:17022/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAACTACGTTATCGCCCAGGCGAACTCCAACTTGGATGAAGAAGGCCACTTCGTAGAAGACCTGGTAACTTGCCGTAGCAAAGGCGAATCCAGCTTGTTCAGCCGCGACCAGGTTGACTACATGGACGTATCCACCCAGCAGGTGGTATCCGTCGGTGCGTCCCTGATCCCGTTCCTGGAACACGATGACGCCAACCGTGCATTGATGGGTGCGAACATGCAACGTCAGGCCGTTCCGACTCTGCGCGCTGATAAGCCGCTGGTTGGTACTGGTATGGAACGTGCTGTTGCCGTT  >  NC_000913/4183080‑4183374

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑