Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 163,417 A→T 10.7% K438I (AAA→ATA) ‡ hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913163,4170AT10.7% 175.5 / 8.6 71K438I (AAA→ATA) ‡hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base T (4/4);  total (35/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 = 9.82e-01

GTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGC  >  NC_000913/163276‑163563
                                                                                                                                             |                                                                                                                                                  
gTCTGTTACAAATGCTGGGGGAACTGGATGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgcc                                                                                                                                             <  2:583933/149‑1 (MQ=255)
 tCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgccg                                                                                                                                            >  2:557189/1‑149 (MQ=255)
  cTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTgccgc                                                                                                                                           <  1:426343/149‑1 (MQ=255)
   tGTCACAAATGCCGGGGGCACTGGAGGGGGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCACCGCGGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcg                                                                                                                                          <  2:580386/149‑2 (MQ=255)
    gTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGAGCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcgg                                                                                                                                         <  1:504024/149‑3 (MQ=255)
    gTTACAAAGGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCa                                                                                                                                         <  1:51926/149‑1 (MQ=255)
      aacacATGCTGGGGGCACGGGAGGGGGACCGGCGTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGCTCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcggta                                                                                                                                       <  1:791762/148‑5 (MQ=255)
         aaaTGCTGGGGGGACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGGTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAATTTTTgccgcggtagca                                                                                                                                    <  2:552996/149‑8 (MQ=255)
         aaaTGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATtctc                                                                                                                                    >  2:852152/1‑149 (MQ=255)
           attcggtgggCAATGGAGTGTGAACGGCTTATTGCGCAAGGGCAAAAACTGGCAGCGCTGCTTAACGATCCCAGTTTAGCTGCAATGCTGGTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGaa                                                                                                                                 <  1:214230/142‑1 (MQ=255)
             gCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAag                                                                                                                                >  2:511443/1‑149 (MQ=255)
             gCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGATGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAag                                                                                                                                >  1:939949/1‑149 (MQ=255)
                    gCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAgccgcc                                                                                                                         >  1:788928/1‑149 (MQ=255)
                          gAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGAt                                                                                                                   >  2:435091/1‑149 (MQ=255)
                           aGGGGGAACGGCTTAGTGCGCAAGGGCAAAAACTGGCCGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATg                                                                                                                  <  2:739790/149‑1 (MQ=255)
                            gggTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCAGCCACGGATgg                                                                                                                 >  1:511182/1‑149 (MQ=255)
                                gAACGGCTTCGTGCGCAAGGGCAAAAAATGCCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGATACCGCGAAGCACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCaa                                                                                                             <  2:549224/149‑1 (MQ=255)
                                 aaCGGCGTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTCACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAAt                                                                                                            <  1:735730/149‑1 (MQ=255)
                                           tGCGAAAGGGCAAAAAATCGCAGCCCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCAGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTgg                                                                                                  <  1:706537/149‑1 (MQ=255)
                                               cAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGCCGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGCTGGGCAATAGTGACCTGggcgt                                                                                              >  1:828298/1‑149 (MQ=255)
                                               cAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGggcgt                                                                                              >  1:607195/1‑149 (MQ=255)
                                                 aGGGCAAAAAATGGCAGCTCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTATCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAATAGCCGCCACGGATGGGCAATAGTGACCTGggcgtgg                                                                                            <  2:648478/149‑1 (MQ=255)
                                                  gcgCAAAAAATGGCAGCGCTGGGCACCGAGCCGCGTTTAGCGGCAATGCTGGTTAGCGCGACGAACGACGACGCAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGAGGGGCAATAGTGACCTGggcgtggc                                                                                           <  2:952432/147‑1 (MQ=255)
                                                      aaaaaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGttt                                                                                       >  1:987508/1‑149 (MQ=255)
                                                      aaaaaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGGTGCGttt                                                                                       >  2:210556/1‑149 (MQ=255)
                                                        aataTGGCAGCGCTGCGGAACGATCCGCGATTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCGGCTACCGCGGCAAAAATTGCCGCGATTCTCGACGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTc                                                                                     <  1:506960/146‑1 (MQ=255)
                                                        aaaaTGTCAGCGCGGTGTAACGAGCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTc                                                                                     <  2:929112/149‑1 (MQ=255)
                                                        aaaaTGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAATGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGACATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTc                                                                                     >  2:901311/1‑149 (MQ=255)
                                                          aaTGGCAGCGCTGGGTAACGATCCGCGGTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTcgc                                                                                   <  1:284230/149‑1 (MQ=255)
                                                               cAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAAt                                                                              >  2:878312/1‑149 (MQ=255)
                                                                       ggTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGATGACGAAGCTGCTACCGCGGCAAAAAATGACGCCAATATCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGTGCGTGGCGTTTTCGCGCAATCAACCAGc                                                                      >  2:618139/1‑149 (MQ=255)
                                                                         tAACGATCCGCGTTTAGCGGCAATGCTGGGTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCt                                                                    <  2:889992/149‑1 (MQ=255)
                                                                             gATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGgca                                                                >  2:941716/1‑149 (MQ=255)
                                                                             gATCCGCGCTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGGGGCGTTTTCGCGCAATCAACCAGCCTGgca                                                                <  1:599254/149‑1 (MQ=255)
                                                                               tCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGAGTTTTCGCGCCATCAACCAGCCTGgcagc                                                              >  1:368135/1‑149 (MQ=255)
                                                                                 cgcgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAGTCAACCAGCCTGGCAGCaa                                                            >  2:468749/1‑149 (MQ=255)
                                                                                 cgcgTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCaa                                                            <  1:435091/149‑1 (MQ=255)
                                                                                     tttAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTa                                                        >  2:786151/1‑149 (MQ=255)
                                                                                          cGGCAATGCTGGTGAGCGCCAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAg                                                   <  1:89799/149‑1 (MQ=255)
                                                                                           ggCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCGGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGCCCTGGGCGTGGCGTGTTCGCGCAATCAACCAGCCTGGCAGCAACGTCGTCAGc                                                  >  1:62370/1‑149 (MQ=255)
                                                                                               aTGCTGGGTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGTATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACt                                              <  2:888342/149‑1 (MQ=255)
                                                                                                tGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTg                                             >  2:739548/1‑149 (MQ=255)
                                                                                                 gCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGt                                            <  2:656518/149‑1 (MQ=255)
                                                                                                 gCTGGTTAGCGCGAAGAACGACGACCAAGCTGCTACCGCGGCAATAATTGCCGCCATTATCGAAGAGCCGCCACGGATGGGCAATACTGCCCTCTGAGTCGCGTTTACGCGCCAGGCTCCAACCGGGCAGCACCCTCGACAGCAACGGt                                            >  1:243290/1‑149 (MQ=255)
                                                                                                       tAGCGCGAAGAACGACGACGACGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAAc                                      <  2:99934/149‑1 (MQ=255)
                                                                                                             gaagaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTACGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGATAAAACGCGTaa                                >  1:50803/1‑149 (MQ=255)
                                                                                                             gaagaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCACGAAGAGCCGCCACGCATGGGCAATAGTGACATCGGCGTGGCGGTTTCGCGCAAACAACCAGACAGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTaaa                                >  1:672179/1‑149 (MQ=255)
                                                                                                              aagaaCGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTaaa                               <  1:374482/149‑1 (MQ=255)
                                                                                                                 aaagacgacgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAacgt                            <  1:436616/146‑1 (MQ=255)
                                                                                                                        ccgaAGCTGCTACCGCGTCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGc                     <  1:929202/148‑1 (MQ=255)
                                                                                                                         cgaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTTGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCg                    >  1:264224/1‑149 (MQ=255)
                                                                                                                          gaAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCgg                   >  1:551428/1‑149 (MQ=255)
                                                                                                                           aaGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAATAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGt                  <  2:379230/149‑1 (MQ=255)
                                                                                                                            aGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCACAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGGGGTg                 >  2:555787/1‑149 (MQ=255)
                                                                                                                             gctgctACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGa                <  1:453325/149‑1 (MQ=255)
                                                                                                                               tgctACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAgg              <  2:281827/149‑1 (MQ=255)
                                                                                                                                gctACCGCGGCAAAAATAGCCGCCATTCTCGAAGAGCCGCCACGGACGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGc             <  1:76619/149‑1 (MQ=255)
                                                                                                                                 ctACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCa            <  2:52480/149‑1 (MQ=255)
                                                                                                                                  tACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTACAACGCTTAAACGTACGTGGCGGTGAGGCAg           >  2:977825/1‑149 (MQ=255)
                                                                                                                                   aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGa          <  1:852152/149‑1 (MQ=255)
                                                                                                                                   aCCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGa          <  2:607195/149‑1 (MQ=255)
                                                                                                                                    tggcgGCAATTTTTGCCTCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACATGGGCGTGGCGTTTTCGCGCAATCAACCATCCTGGCAGAAACGTAGGCAGCAACTGTTAAAACGCTGAAACGTACGTGGCGGTCAGGCAGAc         >  2:705175/3‑149 (MQ=255)
                                                                                                                                    tggcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCCACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGAc         >  1:917696/3‑149 (MQ=255)
                                                                                                                                    ccGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGAc         <  1:511443/149‑1 (MQ=255)
                                                                                                                                    ccGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGGGCCCTGGGCGTGGCGTTTTCGCGCAATCAACCAGGCTGGGCGGAAGGTAGTCAGCAACTGTTAAAGCGCTTAAACGTACGTGGCGGTGAGGCAgcc         >  1:116720/1‑147 (MQ=255)
                                                                                                                                     cgcgGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTGTTCGCGCAATCACCCAGCCTGGCAGCAACGGAGTCAGCAACTGGTAAGCCGCTTACACGTCCGGGGCGGTGAGGCAGACa        >  1:567128/1‑149 (MQ=255)
                                                                                                                                      gcgGCAATTTTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACATGGGCGTGGCGTTTTAGCGCAATCAACCCGCCTGGCAGCAACGTAGTCACCAACAGATAAAACGCTTAAACGTAAGTAGCGGTCAGGCAGACAg       >  2:787412/1‑149 (MQ=255)
                                                                                                                                         gCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTc    >  1:877748/1‑149 (MQ=255)
                                                                                                                                         gCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGACATAGTGCCCTGCTCGTTGCGGGTTCGCGCAATCACGCATGCTGGCAGGAACGTAGGCCGCAACTGGTATAAAGCTTAAAACTACGTGGCGGTGAGTACGACAGtta    >  2:208519/1‑148 (MQ=255)
                                                                                                                                          cAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATCGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCg   <  1:123755/149‑1 (MQ=255)
                                                                                                                                          cAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCg   >  2:649565/1‑149 (MQ=255)
                                                                                                                                           aaaaaTTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGc  <  2:404483/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCAAGGGCAAAAAATGGCAGCGCTGGGTAACGATCCGCGTTTAGCGGCAATGCTGGTTAGCGCGAAGAACGACGACGAAGCTGCTACCGCGGCAAAAATTGCCGCCATTCTCGAAGAGCCGCCACGGATGGGCAATAGTGACCTGGGCGTGGCGTTTTCGCGCAATCAACCAGCCTGGCAGCAACGTAGTCAGCAACTGTTAAAACGCTTAAACGTACGTGGCGGTGAGGCAGACAGTTCGC  >  NC_000913/163276‑163563

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