Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 446,787 G→A 6.2% intergenic (‑121/+28) yajR ← / ← cyoE putative transporter/protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913446,7870GA 6.2% 875.5 / 29.2 259intergenic (‑121/+28)yajR/cyoEputative transporter/protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (102/140);  new base A (7/9);  total (109/150)
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

Only 100 of 262 total aligned reads displayed.

CACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGT  >  NC_000913/446639‑446936
                                                                                                                                                    |                                                                                                                                                     
   cTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGc                                                                                                                                                   <  2:656913/150‑1 (MQ=255)
   cTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGc                                                                                                                                                   <  2:2168218/150‑1 (MQ=255)
        gTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtt                                                                                                                                              >  1:3726099/1‑150 (MQ=255)
        gTCATTTAATAATCGTTCATTACCACCTCTGAATTAAATACGCGACTAGTGTAAAGCGGAAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCAAAGAAAAAAGGTGCCGTAGCACAtt                                                                                                                                              >  1:652349/1‑150 (MQ=255)
         tCATTTTATACTCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAACTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCttt                                                                                                                                             <  1:845067/150‑1 (MQ=255)
         tCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCttt                                                                                                                                             <  2:1946602/150‑1 (MQ=255)
         tCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            <  2:5371213/150‑1 (MQ=255)
          cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCtttt                                                                                                                                            <  2:2786847/150‑1 (MQ=255)
          cATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCtttt                                                                                                                                            <  1:2597510/150‑1 (MQ=255)
            ttttCTACTCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGAAAAGCGTCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCTACGGCACCtttttt                                                                                                                                          <  2:4216742/150‑2 (MQ=255)
            ttttATAATCGTTCATTACCACCTCTGTTTTAAATTCTCGACTAGTGTAAAGCGGCAAACGAATGTGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTa                                                                                                                                          <  1:3141484/150‑1 (MQ=255)
              ttATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAAt                                                                                                                                        <  1:5611125/150‑1 (MQ=255)
              ttATAATCGTTCATTACCACCTCTGTTTTAAAGTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTACATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAAt                                                                                                                                        <  2:4761180/150‑1 (MQ=255)
               tataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATa                                                                                                                                       <  1:3315273/150‑1 (MQ=255)
                ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATag                                                                                                                                      <  1:2801173/150‑1 (MQ=255)
                ataATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATag                                                                                                                                      >  1:3964292/1‑150 (MQ=255)
                 taATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATaga                                                                                                                                     <  1:5428643/150‑1 (MQ=255)
                 taATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATaga                                                                                                                                     <  2:5156679/150‑1 (MQ=255)
                   aTCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagaga                                                                                                                                   <  1:2332090/150‑1 (MQ=255)
                    tCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagagag                                                                                                                                  >  1:1956685/1‑150 (MQ=255)
                    tCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATagagag                                                                                                                                  <  1:3609339/150‑1 (MQ=255)
                       ttCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGtt                                                                                                                               >  2:2566390/1‑150 (MQ=255)
                            tACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGtt                                                                                                                          <  1:2373388/150‑1 (MQ=255)
                              ccaccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAc                                                                                                                        <  1:871439/150‑1 (MQ=255)
                               caccTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTAcc                                                                                                                       <  2:4975378/150‑1 (MQ=255)
                                   tctGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcacac                                                                                                                   >  1:4912015/1‑150 (MQ=255)
                                    ctGTTTTAAATTCGCGACTAGTGTACAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACcacaca                                                                                                                  <  1:5658387/150‑1 (MQ=255)
                                      gTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagc                                                                                                                <  2:2625131/150‑1 (MQ=255)
                                         ttAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagc                                                                                                             >  2:1014571/1‑150 (MQ=255)
                                              ttCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagca                                                                                                        <  2:1698170/150‑1 (MQ=255)
                                              ttCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagca                                                                                                        >  2:2388798/1‑150 (MQ=255)
                                                cgcgACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagcagc                                                                                                      <  2:3224803/150‑1 (MQ=255)
                                                cgcgACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACAcagcagccagcagc                                                                                                      <  2:770454/150‑1 (MQ=255)
                                                  cgACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGt                                                                                                    >  2:528211/1‑150 (MQ=255)
                                                     cTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATg                                                                                                 <  1:5185397/150‑1 (MQ=255)
                                                      tAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGc                                                                                                <  1:2666569/150‑1 (MQ=255)
                                                         tgtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGt                                                                                            >  2:1926629/1‑150 (MQ=255)
                                                          gtAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAATAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGt                                                                                            >  1:5379011/1‑150 (MQ=255)
                                                           tAAAGCGGCAAACGACTGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTc                                                                                           <  1:4758479/150‑1 (MQ=255)
                                                           tAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTc                                                                                           <  2:3448246/150‑1 (MQ=255)
                                                           tAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTc                                                                                           >  1:4896381/1‑150 (MQ=255)
                                                             aaGCGGCAAACGAATGGGGTAAACATGTGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCg                                                                                         <  1:719521/150‑1 (MQ=255)
                                                             aaGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTTGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCg                                                                                         <  1:3150803/150‑1 (MQ=255)
                                                             aaGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCg                                                                                         >  2:3504096/1‑150 (MQ=255)
                                                             aaGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCg                                                                                         <  1:224151/150‑1 (MQ=255)
                                                               gCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGt                                                                                       <  2:4139024/150‑1 (MQ=255)
                                                                cGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGtc                                                                                      >  1:1235594/1‑149 (MQ=255)
                                                                cGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTa                                                                                      <  2:5130992/150‑1 (MQ=255)
                                                                 ggCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTAc                                                                                     >  2:2419130/1‑150 (MQ=255)
                                                                   cAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCa                                                                                   <  2:1571201/150‑1 (MQ=255)
                                                                   cAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCa                                                                                   >  1:2157103/1‑150 (MQ=255)
                                                                   cAAACGAATGGGGTAAAAATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGACGTAGCACCTTTTTAATAGATAGGTTTTGTTACCACACAGAAGCCAGCAGCGTATGCGAGTCCGGTACCa                                                                                   >  2:141963/1‑150 (MQ=255)
                                                                    aaaCGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCAt                                                                                  >  1:1601196/1‑150 (MQ=255)
                                                                       cGAATGGTGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATaaa                                                                               <  2:5521327/150‑1 (MQ=255)
                                                                        gAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATaaaa                                                                              <  2:4912015/150‑1 (MQ=255)
                                                                           tGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCa                                                                           >  2:4623224/1‑150 (MQ=255)
                                                                            ggggTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCaa                                                                          <  1:966001/150‑1 (MQ=255)
                                                                             gggTAAACAGGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAAc                                                                         <  2:589863/150‑1 (MQ=255)
                                                                                tAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGa                                                                      <  2:311691/150‑1 (MQ=255)
                                                                                 aaaCATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGAc                                                                     <  2:934393/150‑1 (MQ=255)
                                                                                   aCATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGAcat                                                                   <  1:3937614/150‑1 (MQ=255)
                                                                                    cATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTATTGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGAcatc                                                                  <  1:1521825/150‑1 (MQ=255)
                                                                                        cGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACAtcatca                                                              >  1:5815572/1‑150 (MQ=255)
                                                                                              ttttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCgaga                                                        >  1:4365069/1‑150 (MQ=255)
                                                                                               tttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACAACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGagag                                                       >  1:2271262/1‑150 (MQ=255)
                                                                                                ttAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAgg                                                      >  1:1575194/1‑150 (MQ=255)
                                                                                                   gTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCa                                                   <  2:3265472/150‑1 (MQ=255)
                                                                                                   gTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCa                                                   <  2:3897651/150‑1 (MQ=255)
                                                                                                    tAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGc                                                                                                <  1:2001346/104‑1 (MQ=255)
                                                                                                    tAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAg                                                  >  1:2151287/1‑150 (MQ=255)
                                                                                                     aaaTATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGt                                                 >  1:5097005/1‑150 (MQ=255)
                                                                                                        tatTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGAt                                              >  1:3577964/1‑150 (MQ=255)
                                                                                                          tttCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATgg                                            <  1:2407417/150‑1 (MQ=255)
                                                                                                              gagaGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCgat                                        >  1:5700581/1‑150 (MQ=255)
                                                                                                              gagaGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCgat                                        <  2:4297919/150‑1 (MQ=255)
                                                                                                                   ggATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATag                                   >  1:3814435/1‑150 (MQ=255)
                                                                                                                    gATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATaga                                  <  2:3655018/150‑1 (MQ=255)
                                                                                                                       atGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGaa                               <  1:4563735/150‑1 (MQ=255)
                                                                                                                       atGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGaa                               <  2:785053/150‑1 (MQ=255)
                                                                                                                       atGTTTATAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGaa                               <  1:5261083/150‑1 (MQ=255)
                                                                                                                            tCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGGGAGGGCAGTGATGGCGATGATAGAGAAGCCGa                          >  1:3471414/1‑150 (MQ=255)
                                                                                                                             cTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGaa                         <  1:4645229/150‑1 (MQ=255)
                                                                                                                              tAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAAc                        >  1:4026983/1‑150 (MQ=255)
                                                                                                                                  gCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCt                    >  1:1152308/1‑150 (MQ=255)
                                                                                                                                  gCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCt                    <  2:106355/150‑1 (MQ=255)
                                                                                                                                     aaGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgc                 >  1:4125119/1‑150 (MQ=255)
                                                                                                                                       gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc               >  1:4246009/1‑150 (MQ=255)
                                                                                                                                       gAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgc               >  1:537174/1‑150 (MQ=255)
                                                                                                                                        taaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg              >  1:2926200/2‑150 (MQ=255)
                                                                                                                                        aaaaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcg              >  2:5799696/1‑150 (MQ=255)
                                                                                                                                         aaaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTgcgcgc             >  2:303956/1‑150 (MQ=255)
                                                                                                                                          aaaaGGTGCTACGGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAATCCGAACAGCTTGCGCGcc            >  1:2029799/1‑150 (MQ=255)
                                                                                                                                          aaaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGcc            >  1:1011592/1‑150 (MQ=255)
                                                                                                                                           aaaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGccc           <  1:3378367/150‑1 (MQ=255)
                                                                                                                                            aaGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCa          <  1:4478840/150‑1 (MQ=255)
                                                                                                                                             aGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGACGCCGAACAGCTTGCGCGCCCAg         >  1:3034278/1‑150 (MQ=255)
                                                                                                                                              ggTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGa        <  1:1014571/150‑1 (MQ=255)
                                                                                                                                                  ccGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCt    <  2:4405150/150‑1 (MQ=255)
                                                                                                                                                   cGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTg   <  1:5508390/150‑1 (MQ=255)
                                                                                                                                                   cGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTg   <  1:286087/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACCTGGCGTCATTTTATAATCGTTCATTACCACCTCTGTTTTAAATTCGCGACTAGTGTAAAGCGGCAAACGAATGGGGTAAACATGCGGATTTTTTAGTAAATATTTCGAGAGGGATATGTTTCTAATGCTAAGAAAAAAGGTGCCGTAGCACCTTTTTAATAGAGAGGTTTTGTTACCACACAGCAGCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGT  >  NC_000913/446639‑446936

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