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,5420CA97.4% 219.8 / ‑4.8 77P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (33/41);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184404‑4184680
                                                                                                                                          |                                                                                                                                          
tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGCTACTCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                        <  2:847062/143‑1 (MQ=255)
       tAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                                 <  1:73835/143‑1 (MQ=255)
         aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCg                                                                                                                               >  2:685340/1‑143 (MQ=255)
          aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                              >  1:872392/1‑143 (MQ=255)
          aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGAAATTTCTAAGATCAACCCGATCGAAGATACGCCTTACGATGAAACCGGTACCCCGGCCGACCTCGTCCTGCACCAGCTGCGCGTCCCGt                                                                                                                              >  1:731061/1‑143 (MQ=255)
              cataaccAGCCTGGTGAGAAGAAGGCAGGGCGTCACGGTAACAATGGTGTACTTTCTAAGATCAACCCGATCGAAGAGATGCCTCACGATGAAAACGGTACCCCGGTATACATCGTACTGAACCAGCTGGGCGTACCGTCTCg                                                                                                                          <  1:97962/138‑1 (MQ=255)
                 aTCCAGCCTGGTGACAAGAGGGCCGGTCGTCGCGGAACCAATGTTTTTATTTCTATGATCAACCCTATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                       <  1:150556/143‑1 (MQ=255)
                     aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTCAACCAGCTGGGCGTACCGTCTCGTATGAAc                                                                                                                   >  2:40881/1‑143 (MQ=255)
                      gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACa                                                                                                                  >  1:808987/1‑143 (MQ=255)
                      gCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGAGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCACGTATGAAcg                                                                                                                  >  1:94013/1‑142 (MQ=255)
                        cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACATc                                                                                                                >  1:500103/1‑143 (MQ=255)
                           gTGGCAAGATGGCAGGCCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCCTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                             <  2:744675/143‑1 (MQ=255)
                           gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                             >  2:453446/1‑143 (MQ=255)
                            tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                            <  2:204138/143‑1 (MQ=255)
                            agccAACATGGCAGATAACCACGTTAACAAGCGTGTAATTTCTAAGATCAACCCGATCTAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                            <  1:748026/140‑1 (MQ=255)
                              aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAGTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                          <  1:28441/143‑1 (MQ=255)
                              aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACGTCGGTCAg                                                                                                          >  1:1090406/1‑143 (MQ=255)
                                aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTATCATGAAAACGGTACGCCGGTAGACATCGGAATGAATCAGATAGGCGTAGCGTCCGGTATGGACGTCGGTcggca                                                                                                        >  2:314051/1‑139 (MQ=255)
                                  gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGAAGACATCGTACTGAACCAGCTGGGCGTACCGACCCGTATGAACATccaacagaacc                                                                                                      >  2:8389/1‑133 (MQ=255)
                                   gtaggaGGTGGGCACGGGAACAAGGGGGGAATGTCTAAGATCGACCCGATCGAAGATATGCCTTACGATGAAGGCGGTACGCCGGTAGAGATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCt                                                                                                     <  2:494729/138‑1 (MQ=255)
                                    tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGAGCCTc                                                                                                    >  2:375355/1‑143 (MQ=255)
                                        aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  1:948406/143‑1 (MQ=255)
                                        aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  1:394882/143‑1 (MQ=255)
                                        aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                <  1:1007348/143‑1 (MQ=255)
                                         ggTCGTTACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                               <  1:1008319/143‑1 (MQ=255)
                                          gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGAACGCCGGTAGCGCTGGAGCAGAGGCCGCAGGTCGTACCGCCTCGTATGGACATCGGCCAGATCCtacaaccc                                                                                              >  1:509954/1‑136 (MQ=255)
                                              tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCGTTAACATGAAAACGAAACGAAGGTAGACACCGAACCGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGTTTCTCGAAACCCAct                                                                                          >  1:446941/1‑142 (MQ=255)
                                              tcACGGTAACAAGGGCGTAATTGCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                                                          <  1:244846/143‑1 (MQ=255)
                                                aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                        >  1:637650/1‑143 (MQ=255)
                                                aCGGTAACAAGGGTGTAATTTCCAAGCTCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                        <  1:201663/143‑1 (MQ=255)
                                                 cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGTTGAAAACGGTACGCCAGTAGACATCGTACTGAACCAGCTGGGCGAACAGTCTCGTATGGACATCGCTAGGATCGTCGAAACCCAcacag                                                                                       >  2:1069195/1‑139 (MQ=255)
                                                     aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGACGATATGCCTTACGATGGAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTTTAGGATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTAt                                                                                   <  2:338297/143‑1 (MQ=255)
                                                      acagtcttgtatcttCTAAGATGGGCCCGATCGAGTATATGCTTTCCGTTGTGGACGGTACGCCGGTAGACATCGTACTGAGCCAGCTGGGTGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                  <  2:1048981/130‑1 (MQ=255)
                                                        aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                                >  1:1074830/1‑143 (MQ=255)
                                                        aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                                >  2:357174/1‑143 (MQ=255)
                                                        aaGGGTGTAATTCCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGGAGACATCGCACTGACCCAGCTGGGCGTACCGTCTCGTATGAACGTCGGTCAGATCCTCGAAACCAACCGCGGTAGGGc                                                                                >  2:1039613/1‑143 (MQ=255)
                                                         cgggTGTACTTTCTACCATCCACCCGATCGAAGTTATGCCTTACGATGAAACCGCTACCCCCGTAGACTTCGTACTGAACCAGCTTGGCTTCCCACCCCGTACGATCACCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                               <  2:464503/142‑1 (MQ=255)
                                                          ggctggaaTTTCTAAGATCAACCCGATCGAAGATACGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTg                                                                              <  1:843947/137‑1 (MQ=255)
                                                             tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                           >  1:952148/1‑143 (MQ=255)
                                                               tcatgtCTAAGATCAACGCGATCGTAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCTACCTGGGTATGGCTGCGaaa                                                                         <  2:189917/138‑1 (MQ=255)
                                                               tactttCTAAGATCCACCCGATCGAAGATATGCCTTACGATGAAGACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                         <  2:324048/140‑1 (MQ=255)
                                                               tAAGTTCCAAGATCAACCCGAGCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                         <  2:599748/143‑1 (MQ=255)
                                                               tAAGTTCCAACATCACCCCGATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                         <  1:679251/143‑1 (MQ=255)
                                                                aatcacTAAGACCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                        <  1:725212/138‑1 (MQ=255)
                                                                aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGGTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                        <  2:568427/143‑1 (MQ=255)
                                                                aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                        <  1:685340/143‑1 (MQ=255)
                                                                aaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                                                        >  1:622918/1‑143 (MQ=255)
                                                                  tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAACCTCGAAACCCACAGGGTGATGGCTGCggacgat                                                                      >  2:638008/1‑137 (MQ=255)
                                                                    ccTACGATCACCGCGATCGAGCATACGCCTTACGCGGACAAGGCTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGCACCGTCACGTATGAACATCGGCCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                                    <  2:713551/142‑1 (MQ=255)
                                                                      tAAGATCTATCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                  <  1:423548/143‑1 (MQ=255)
                                                                      tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                  >  1:939828/1‑143 (MQ=255)
                                                                       aggATCTACCCGAGCGAAGATATGCCCTACGACGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                                 <  2:108542/141‑1 (MQ=255)
                                                                        aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGc                                                                >  2:255347/1‑143 (MQ=255)
                                                                          aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                              >  1:662893/1‑143 (MQ=255)
                                                                              aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAg                                                          >  1:355900/1‑143 (MQ=255)
                                                                               cttgATCGAAGATATGCCGTACGATCAAAACGTTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                         <  1:921032/140‑1 (MQ=255)
                                                                                    tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCGCGCACGCAATCCGGGTATGGCTCCGCAAGGTATGGGCGagaagacccac                                                    >  2:229352/1‑138 (MQ=255)
                                                                                    tCGAAGATAGGCCGTAGGATGCAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                    <  1:83328/143‑1 (MQ=255)
                                                                                      gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAgcgc                                                  >  1:1022483/1‑143 (MQ=255)
                                                                                       aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGCCGACAAGACCAACGcc                                                 >  2:925346/1‑143 (MQ=255)
                                                                                                cTTACGTGGAAAACGGTACGCCGGTAGACACCGTACTGAGCCGGCTGGGCGTCCCGTCTCGTATGAACATCGGTCAGCTCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                        <  2:181448/143‑1 (MQ=255)
                                                                                                cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                        <  1:609689/143‑1 (MQ=255)
                                                                                                   aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                     <  2:628990/143‑1 (MQ=255)
                                                                                                         aacaccgTACGCCGGGAGACATCGGACTGACCCAGCTGGGCGTACCGTCTCGTATGGACATCGGTCAGATTCTGGACTCTCATCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                               <  2:536574/137‑1 (MQ=255)
                                                                                                            atggTACGGCGGTAGACATCTTACTTAACCAGCTGGCCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                            <  2:94013/141‑1 (MQ=255)
                                                                                                            aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGCAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagca                            >  2:275930/1‑141 (MQ=255)
                                                                                                              ggTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAACGGTATCGGCGACAAGATCAACGCCATGCTGAAAAAGCAGCaagaaat                          >  1:1030496/1‑141 (MQ=255)
                                                                                                                tACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAACCAGCAGCAAGAAGTcg                        >  1:817544/1‑143 (MQ=255)
                                                                                                                 aCGCCGGTAGACATCGTACTGAACCAGCTGGTCGTACCGGCGCGTATGCACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgc                       <  2:713224/143‑1 (MQ=255)
                                                                                                                   gCCGGTAGACGGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                     <  2:119213/143‑1 (MQ=255)
                                                                                                                         ggACATCGTAGAGCACCAGGTTGGCGTACCGGCTCATATGAACATCGGGCAGAGGCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGTTCAACGCCATGCTGAATCAGCAGCAAGAAGTCGCGAAACTgc               <  2:414431/142‑1 (MQ=255)
                                                                                                                             aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCACGAAGTCGCGAAACTGCGCGa           >  1:817701/1‑143 (MQ=255)
                                                                                                                              tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa          <  2:622918/143‑1 (MQ=255)
                                                                                                                               cGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt         <  1:453446/143‑1 (MQ=255)
                                                                                                                                 tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAACGACATGCAAAAACCGAAGCGAGGAGTCGGGAAACCGCGCGAAATc       >  2:992453/1‑143 (MQ=255)
                                                                                                                                    tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc    <  1:174659/143‑1 (MQ=255)
                                                                                                                                      aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCTATGTTGAAACAGAAGCAAGAAGTAGCGAAACTGCGCGAATTGATACa  >  1:194778/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184404‑4184680

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: