Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,542 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5420CA98.7% 224.3 / NA 76P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (39/36);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAG  >  NC_000913/4184402‑4184681
                                                                                                                                            |                                                                                                                                           
tCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                           >  1:83481/1‑143 (MQ=255)
   ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACACCGAGAACGTGACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                        >  2:397999/1‑143 (MQ=255)
    gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGCCGTAACCGTCAGTAATGCGTGATTTCTTCTCTTAACCCATTCGTAGATATGCCTTACGATGAAAACGGTACACCGGTAGACATCGTACTGAACAAGCTggg                                                                                                                                       >  2:700071/1‑143 (MQ=255)
       gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAAAGGTACGCCGGTAGACATCGTACAGAACCCGCGGGGCga                                                                                                                                    >  1:180636/1‑142 (MQ=255)
       gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGACGAACACGGTACGGCGGTAGACATCGTACTGAACCACCTCGGCga                                                                                                                                    >  1:515467/1‑142 (MQ=255)
       gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAAATTCTAAGATGAACCAGATCGAAGATATGCCTTACGATGAAAACGGTACGACGGTAGACATCGTACTGAACCAGCTGGGCGt                                                                                                                                    >  1:163456/1‑143 (MQ=255)
        ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGcaga                                                                                                                                   >  2:664214/1‑140 (MQ=255)
          aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAcc                                                                                                                                 >  2:788405/1‑143 (MQ=255)
            aggcactATCCAGCCTGGTTACAAGATGGCAGGTCGTCACGCGAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACTATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                               <  1:941514/137‑1 (MQ=255)
               ccGTATCCAGCCTGGTGACAAGATGTCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAACTGGGCGTACCGtctc                                                                                                                            <  2:152382/143‑1 (MQ=255)
               ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                            <  2:712259/143‑1 (MQ=255)
                 gTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                          >  2:563433/1‑143 (MQ=255)
                  tacccAGCCTGGTGCCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                         <  2:288340/140‑1 (MQ=255)
                    tCCAGCCGGGTGACAAGATGGCAGGTCGTCACGGGAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACTGTCTCGTATg                                                                                                                       <  2:826465/143‑1 (MQ=255)
                      cAGCCTGGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                     <  1:552807/143‑1 (MQ=255)
                         ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGAAATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAAcgt                                                                                                                  >  2:646637/1‑141 (MQ=255)
                          cTGGTTACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATc                                                                                                                 <  1:742938/143‑1 (MQ=255)
                           tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                                <  2:827822/143‑1 (MQ=255)
                                    gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTCCGATGAAAACGGTACGCCGGTAGACATCGTACTGAGCGAGCTGCGCGAACCGTCTCGCACGAACACCGGGCAGATcc                                                                                                       >  1:419312/1‑143 (MQ=255)
                                      ttgCAGCTCGTCACGGTAACAAGGGTGTAAGTTCTAAGATCAACACGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                     <  1:41819/141‑1 (MQ=255)
                                      tGGTAGGTCGTCAGGGGAACAAGGGTGGAATTCCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                     <  2:118386/143‑1 (MQ=255)
                                      tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                     <  2:533923/143‑1 (MQ=255)
                                       ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                                    >  1:168563/1‑143 (MQ=255)
                                          aGGTCTACACGGTATCACGGTTGTACTTTGTACGCTCAACCCGATCGAAGCTATGCCTTACGCTGAAATCCGTGCTCCTGTCCACCTCGTACTCAACCAGCGGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                 <  1:831403/143‑1 (MQ=255)
                                          aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                 >  2:112250/1‑143 (MQ=255)
                                            gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                               >  2:391689/1‑143 (MQ=255)
                                              cgtcACGGTAACGAGGGTGTAATTTCTAAGATCACCCCGCTCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                             <  1:543892/143‑1 (MQ=255)
                                              cgtcACGGTAACAAGGGTGGAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAGCGGGACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                                                             <  1:769177/143‑1 (MQ=255)
                                                tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATCAACATCGGTCAGATCCTCGAAAcccccc                                                                                           >  1:80956/1‑143 (MQ=255)
                                                taaCGGGAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGACCCTCGAAACCCAcc                                                                                           >  1:166821/3‑143 (MQ=255)
                                                        acaattttttAATTTCTAATATCAACCCGATCGAATATATGCCTTACGATGAAAACTGTACTCCGATAGACATCGTACTGAACCATCTGGTCATACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                   <  1:210495/134‑1 (MQ=255)
                                                        aCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                   <  2:202885/143‑1 (MQ=255)
                                                        aCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTCCCGTCTCGTATGATCATCCGTCACGTCCACGAAACCCACCTGGGTATg                                                                                   >  2:130544/1‑143 (MQ=255)
                                                          aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTCACGATGAAAACGGTACGCAGGTAGACATCGCACTGAACCAGCTCGGCGTAACGTCTCGTAGGGACATCGGTCGGGTGCTCGGCACCGCCCAGGGTACGGc                                                                                 >  1:922403/1‑143 (MQ=255)
                                                           aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                                <  2:726764/143‑1 (MQ=255)
                                                             ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGGCACCCACGTGGGCGTGGCTGc                                                                              >  1:31661/1‑143 (MQ=255)
                                                             ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCAGTAGACAGCGTTCTGAGCCAGCTGGGCGTACCGTCTCGTATGAACAGCGCTGAGATCCCCGAAAGCCACCAGGATATGGccca                                                                              >  1:721099/1‑140 (MQ=255)
                                                              gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCg                                                                             >  1:525695/1‑143 (MQ=255)
                                                                gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCCGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGTTCTGGCTGCGaa                                                                           >  2:878841/1‑143 (MQ=255)
                                                                ggaaTTTCGAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                           <  2:281335/141‑1 (MQ=255)
                                                                 tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                          >  2:783507/1‑143 (MQ=255)
                                                                  aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                         <  1:535926/143‑1 (MQ=255)
                                                                      tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAGGGTAt                                                                     >  1:182911/1‑143 (MQ=255)
                                                                      tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                                     >  2:829945/1‑143 (MQ=255)
                                                                        tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                   >  2:801840/1‑143 (MQ=255)
                                                                           gATCAACCCGATCGAAGATATGCCTTACGATGAAAACCGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGCAACCCACCAGGGTATGGATGCGAAAGGAATCGGCg                                                                >  1:18488/1‑143 (MQ=255)
                                                                            aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCAGGGCGTACCGTCTCGTATGAACATCCGTCAGCTCCTCGAAACCCACCCCGGTCTCGCTGCGGAAGGTATCGGCGa                                                               >  1:76268/1‑143 (MQ=255)
                                                                             tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGCATGAACATCGGTCAGATCCTCGAAACCCACCCGGGTATGGCTGCGAAAGGTATCGGCGAc                                                              >  1:120152/1‑143 (MQ=255)
                                                                               accccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACGTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACATAGATCTGGCGGCGAAAGGTAGCGGCGACaa                                                            <  1:581783/141‑1 (MQ=255)
                                                                                aCCCGATCGAGGATATGCCTTACGATGACAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                           <  1:47724/143‑1 (MQ=255)
                                                                                 cccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                          >  2:168677/1‑143 (MQ=255)
                                                                                   cGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                                                        <  2:839764/143‑1 (MQ=255)
                                                                                     aTCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCaa                                                      <  2:683908/143‑1 (MQ=255)
                                                                                      tctacgatctgtgctACTATGAAAACTGTACGCCTGTACACATCGTACTCAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCTAAAGGTATCGGCGACAAGATCAAc                                                     <  1:919243/129‑1 (MQ=255)
                                                                                      tagagcatatGCTTTGCTATCAAAACGGGAGGCGGGGATAAATCATACTGAACCAGCTGGGCGTACCGTTTCGTATGAACATCGGTCAGATCCTCGAAACCCAACTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                     <  1:61010/137‑1 (MQ=255)
                                                                                      tCGATGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCTAAAGGTATCGGCGACAAGATCAAc                                                     <  1:609401/143‑1 (MQ=255)
                                                                                         aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                                  >  2:732123/1‑143 (MQ=255)
                                                                                         aaGATATGCCTTACGATGAAAACCGTAACACGAAAAAATTGGGACTTAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                                  >  2:909262/1‑143 (MQ=255)
                                                                                          aGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCCGGGTATGGCTGCGAAAGGTATCGGCGACAAGAGCAACACCa                                                 >  2:230625/1‑143 (MQ=255)
                                                                                                  cTTACGATGAAACCGGTGCGCCGGTAGACTTCGTACTGAACCAGCTGCGCGTACCGTCTCGTATGAACATCGGTCAGATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                         <  2:420403/143‑1 (MQ=255)
                                                                                                        aTGAAAGCAGTACGCCGTTAGACATCCTACTCAACCAGCTGCTCCTACCGTCTCGTATGAACATCGGTCAGGTCCTCCAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                                   <  1:658860/143‑1 (MQ=255)
                                                                                                          gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGAATCAGCGACAAGATCAACGCCACGCTGAACcagcaacc                                 >  1:759780/1‑140 (MQ=255)
                                                                                                           aaaaCGGTACGCCGGTAGACATCGTACTGAACCACCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                                <  2:863360/143‑1 (MQ=255)
                                                                                                                 ggccggcGGTAGACATCGGACTGAACCAGCTGGGCGGACCGTCTCGGATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTc                          <  1:937413/137‑1 (MQ=255)
                                                                                                                   aCGCCGTTAGGCGTCGTACTGCACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                        <  1:688780/143‑1 (MQ=255)
                                                                                                                         gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACt                  >  1:788305/1‑143 (MQ=255)
                                                                                                                         gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGCATCGGCGACAAGATCAACGCCATGCTGAAAAAGCAGCAAGGAGTCGCGAAACt                  >  1:276909/1‑143 (MQ=255)
                                                                                                                             aCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGGCAAGATCAACGCCATGCCGAAACCGCAGCACGAAGCCACGAAACTgcgc              >  2:868897/1‑143 (MQ=255)
                                                                                                                              cATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGcgcg             <  1:130544/143‑1 (MQ=255)
                                                                                                                                 cGTACTGACCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt          <  2:296112/143‑1 (MQ=255)
                                                                                                                                  gTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAtt         <  2:525695/143‑1 (MQ=255)
                                                                                                                                        aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                   <  2:228979/79‑1 (MQ=255)
                                                                                                                                        aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                   >  1:228979/1‑79 (MQ=255)
                                                                                                                                         aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGGGAAACTGCGCGAATTCATCCAg  >  2:692196/1‑143 (MQ=255)
                                                                                                                                         aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCaacaag  >  2:211381/1‑138 (MQ=255)
                                                                                                                                         aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAg  <  2:341930/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
TCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAG  >  NC_000913/4184402‑4184681

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