Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,184,497 A→C 23.5% M1085L (ATG→CTG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,4970AC23.5% 231.8 / 64.3 122M1085L (ATG→CTG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (42/51);  new base C (15/14);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAG  >  NC_000913/4184354‑4184646
                                                                                                                                               |                                                                                                                                                     
cccAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCtt                                                                                                                                                 >  1:1890203/1‑150 (MQ=255)
  cAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGCTCAACCCGATCGAAGATATGCCTTAc                                                                                                                                               >  1:428903/1‑150 (MQ=255)
  cAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTAc                                                                                                                                               >  1:1882632/1‑150 (MQ=255)
  cAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTAc                                                                                                                                               <  2:1197435/150‑1 (MQ=255)
  cAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTAc                                                                                                                                               >  2:853524/1‑150 (MQ=255)
   aGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAATCGCCGTCTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGAAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACg                                                                                                                                              >  1:370716/1‑150 (MQ=255)
    gggCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGa                                                                                                                                             >  1:537300/1‑150 (MQ=255)
     ggCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACACGATCGAAGATATGCCTTACGAt                                                                                                                                            >  2:2114837/1‑150 (MQ=255)
     ggCGACGATCTGGCACCGGGCGAGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGTCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCTACCCGATCGAAGATATGCCTTACGAt                                                                                                                                            <  1:132035/150‑1 (MQ=255)
       cgacgaTCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGa                                                                                                                                          <  2:1160441/150‑1 (MQ=255)
        gacgaTCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGaa                                                                                                                                         <  1:1241083/150‑1 (MQ=255)
           gaTCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAAc                                                                                                                                      <  1:1914641/150‑1 (MQ=255)
           gaTCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAAc                                                                                                                                      <  2:1476299/150‑1 (MQ=255)
               tGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTa                                                                                                                                  <  1:2094293/150‑1 (MQ=255)
                 gCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACg                                                                                                                                >  1:466714/1‑150 (MQ=255)
                        gCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAg                                                                                                                         <  2:639456/150‑1 (MQ=255)
                             cTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATc                                                                                                                    >  1:853097/1‑150 (MQ=255)
                             cTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATc                                                                                                                    >  2:1638607/1‑150 (MQ=255)
                               gtaGATTGTTAAGGTATATCTGGAGGTTAAACGCCGTATCCAGCCTTTTGACAAGATGGCATGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACACGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGt                                                                                                                  <  1:332097/148‑1 (MQ=255)
                               gAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGt                                                                                                                  >  1:846708/1‑150 (MQ=255)
                                 aGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAc                                                                                                                >  1:358113/1‑150 (MQ=255)
                                 aGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAc                                                                                                                >  1:1627925/1‑150 (MQ=255)
                                 aGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAc                                                                                                                >  1:1485335/1‑150 (MQ=255)
                                 aGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTAc                                                                                                                >  2:76312/1‑150 (MQ=255)
                                  gATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACt                                                                                                               <  2:699797/150‑1 (MQ=255)
                                      gTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAc                                                                                                           >  1:640474/1‑150 (MQ=255)
                                      gTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAc                                                                                                           >  1:298855/1‑150 (MQ=255)
                                       ttAAGGTATATCTGGAGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTTTATTTTCTATGTTCAACGCGATCGAAGATATGCCTTACGATGATTACGGTACGCCGGTAGTCATCGTACTGAAcc                                                                                                          <  2:1275704/150‑1 (MQ=255)
                                          aGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGc                                                                                                       <  1:561234/150‑1 (MQ=255)
                                           ggTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCt                                                                                                      >  1:161759/1‑150 (MQ=255)
                                                atCTGGCTGTTAAACGCCGATTCCAGCCTGGTGACAAGATGGCAGGTCGTCTCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAATACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCg                                                                                                 <  1:1759711/150‑1 (MQ=255)
                                                     gCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCg                                                                                            <  2:1641583/150‑1 (MQ=255)
                                                          tAAACGCCGTAACCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATCCCTAACGATGAAAACGGTACGCCGGTAGACATCGAACTGAACCCGCTGGGCGTACCGTCTCg                                                                                       >  1:2180650/1‑150 (MQ=255)
                                                            aaCGCCGTATCTAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTa                                                                                     >  2:1457879/1‑150 (MQ=255)
                                                            aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTa                                                                                     <  2:527853/150‑1 (MQ=255)
                                                            aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTa                                                                                     <  2:261507/150‑1 (MQ=255)
                                                            aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTa                                                                                     >  2:1169752/1‑150 (MQ=255)
                                                            aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGGGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTa                                                                                     <  1:720803/150‑1 (MQ=255)
                                                              cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATg                                                                                   <  1:1003181/150‑1 (MQ=255)
                                                               gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGa                                                                                  <  1:819207/150‑1 (MQ=255)
                                                                ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGACAACGGTACGCCGGTAGACGTCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGaa                                                                                 >  2:104567/1‑150 (MQ=255)
                                                                 cGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAAc                                                                                <  1:1535766/150‑1 (MQ=255)
                                                                     tCCAGCCTGGTGACAACATGGCAGGGCGTCACGGAAACAAGGGGGTAATTGCTAAGACCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCg                                                                            <  2:318874/150‑1 (MQ=255)
                                                                     tCCAGCCTGGGGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGAACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCg                                                                            <  2:1309097/150‑1 (MQ=255)
                                                                        aGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                         <  2:1485335/150‑1 (MQ=255)
                                                                          ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATAGAATATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                       <  1:1083996/150‑1 (MQ=255)
                                                                            tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGAAATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                     <  1:1075354/150‑1 (MQ=255)
                                                                             ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                    >  1:496045/1‑150 (MQ=255)
                                                                             ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                    >  1:1788976/1‑150 (MQ=255)
                                                                             ggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATc                                                                    <  2:1428499/150‑1 (MQ=255)
                                                                              gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTAGTATGAACATCGGTCAGATcc                                                                   >  2:1958240/1‑150 (MQ=255)
                                                                                gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                 <  1:605612/150‑1 (MQ=255)
                                                                                gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAACAACTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCACTGGGCGTACAGTCTCGTATTAACATCTGTCAGATCCTc                                                                 >  2:1395308/1‑150 (MQ=255)
                                                                                gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGAGCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                 <  1:423613/150‑1 (MQ=255)
                                                                                 aCAAGATGGCAGGTCGTCACGGAAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                <  1:1333736/150‑1 (MQ=255)
                                                                                 aCAAGAGGGCAGGTCGTCACGCTAACAAGGCTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGGACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                <  2:1308026/150‑1 (MQ=255)
                                                                                  cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCGGGGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGa                                                               >  2:1365421/1‑150 (MQ=255)
                                                                                  cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                               >  2:2071716/1‑150 (MQ=255)
                                                                                     gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                            <  1:766024/150‑1 (MQ=255)
                                                                                     gATGGCAGGTCGGCACGGTAACTAGGGTGTAATTTCTAAGATCAACCCTATCGATGTTCTGCCTTACGATGAATACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCTAAAc                                                            <  2:1842751/150‑1 (MQ=255)
                                                                                      aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                           <  1:1689410/150‑1 (MQ=255)
                                                                                      aTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                           >  1:1262914/1‑150 (MQ=255)
                                                                                       tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAccc                                                          >  1:1247249/1‑150 (MQ=255)
                                                                                        ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCa                                                         <  2:1249264/150‑1 (MQ=255)
                                                                                          cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcc                                                       >  1:478839/1‑150 (MQ=255)
                                                                                             gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                    >  1:2107652/1‑150 (MQ=255)
                                                                                              tcgtcACTGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                   <  1:1819835/150‑1 (MQ=255)
                                                                                              gcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                   <  1:406687/149‑1 (MQ=255)
                                                                                               cgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGt                                                  <  2:1787368/150‑1 (MQ=255)
                                                                                                  cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                               >  2:1146808/1‑150 (MQ=255)
                                                                                                    cGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                             <  2:298855/150‑1 (MQ=255)
                                                                                                        aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCg                                         <  2:578326/150‑1 (MQ=255)
                                                                                                        aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCg                                         <  1:497397/150‑1 (MQ=255)
                                                                                                        aaCAAGGGTGTAATTTCTAAGATAAACCCGATCGAAGATCTGCGTTTCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGTGTATGGCTGCg                                         <  2:548706/150‑1 (MQ=255)
                                                                                                          cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                       >  2:625123/1‑150 (MQ=255)
                                                                                                            aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                     <  2:900671/150‑1 (MQ=255)
                                                                                                            aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAg                                     >  1:1384458/1‑150 (MQ=255)
                                                                                                             gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAgg                                    <  2:2035340/150‑1 (MQ=255)
                                                                                                             gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCGGCGAAAgg                                    <  2:1218296/150‑1 (MQ=255)
                                                                                                               gtgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTa                                  >  1:2179463/1‑150 (MQ=255)
                                                                                                                tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                 >  1:1032860/1‑150 (MQ=255)
                                                                                                                 gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                >  2:1613510/1‑150 (MQ=255)
                                                                                                                    aTTTCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                             >  1:709295/1‑150 (MQ=255)
                                                                                                                       tCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCTAAAGGTATCGGCGAc                          <  2:161759/150‑1 (MQ=255)
                                                                                                                       tCTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                          <  2:1882632/150‑1 (MQ=255)
                                                                                                                       tCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGAc                          <  2:1021601/150‑1 (MQ=255)
                                                                                                                        cTAAGATCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                         >  2:1079860/1‑150 (MQ=255)
                                                                                                                        cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                         >  2:1173753/1‑150 (MQ=255)
                                                                                                                        cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACa                         >  2:2178042/1‑150 (MQ=255)
                                                                                                                            taTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                     <  1:1324840/149‑1 (MQ=255)
                                                                                                                              tCAACCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa                   <  2:894158/150‑1 (MQ=255)
                                                                                                                              tCAACCCGATCGAAGATCTGCCTTACGATGAAAAAGGTACGCCGGTTGACTTCTTACTTAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACGCACCTTGGTATGGCTGCGAAATGTATCGGCGACAAGATCa                   <  2:894161/150‑1 (MQ=255)
                                                                                                                              tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa                   >  1:673061/1‑150 (MQ=255)
                                                                                                                                aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                 >  1:581826/1‑150 (MQ=255)
                                                                                                                                 aCCCGATCGAAGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACg                >  2:1801171/1‑150 (MQ=255)
                                                                                                                                       tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc          <  1:1250575/150‑1 (MQ=255)
                                                                                                                                         gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGGATGCACATCGGTCAGATCCTCGAATCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGACATGCTg        >  1:1972290/1‑150 (MQ=255)
                                                                                                                                          aaGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGa       <  2:537300/150‑1 (MQ=255)
                                                                                                                                           aGATCTGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaa      >  1:904958/1‑150 (MQ=255)
                                                                                                                                              tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAca   >  1:603098/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CCCAGGGCGACGATCTGGCACCGGGCGTGCTGAAGATTGTTAAGGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAG  >  NC_000913/4184354‑4184646

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: