Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,578,549 A→G C263C (TGT→TGC mcrB ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,578,5490AG99.2% 479.0 / NA 121C263C (TGT→TGCmcrB5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (70/50);  total (70/51)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGC  >  NC_000913/4578406‑4578695
                                                                                                                                               |                                                                                                                                                  
gTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAAt                                                                                                                                             >  2:213288/1‑151 (MQ=255)
  tAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAAtt                                                                                                                                            <  1:215844/150‑1 (MQ=255)
       gAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCGTTAGCTTGCTGGCAAAAATTGTaaaa                                                                                                                                      <  2:939724/151‑1 (MQ=255)
         aCAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAata                                                                                                                                    <  2:798246/151‑1 (MQ=255)
                 gTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCt                                                                                                                            <  2:1223834/151‑1 (MQ=255)
                  ttttCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCtt                                                                                                                           >  2:198701/1‑151 (MQ=255)
                   tttCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCtt                                                                                                                           <  2:351566/150‑1 (MQ=255)
                     tCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTa                                                                                                                         <  1:1165517/150‑1 (MQ=255)
                      cACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTAc                                                                                                                        <  2:59237/150‑1 (MQ=255)
                      cACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACg                                                                                                                       <  2:167080/151‑1 (MQ=255)
                        ccTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTc                                                                                                                     >  2:1287830/1‑151 (MQ=255)
                         cTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCg                                                                                                                    >  1:969166/1‑151 (MQ=255)
                           cGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCg                                                                                                                    <  2:144107/149‑1 (MQ=255)
                           cGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGa                                                                                                                  <  2:860103/151‑1 (MQ=255)
                            gTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGaa                                                                                                                 <  1:207533/151‑1 (MQ=255)
                            gTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGaa                                                                                                                 <  1:958090/151‑1 (MQ=255)
                            gTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGa                                                                                                                  >  1:823163/1‑150 (MQ=255)
                              ttATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGc                                                                                                               >  1:679756/1‑151 (MQ=255)
                              ttATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGc                                                                                                               >  1:368655/1‑151 (MQ=255)
                              ttATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGc                                                                                                               >  2:316572/1‑151 (MQ=255)
                               tATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGcc                                                                                                              >  2:198051/1‑151 (MQ=255)
                               tATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGc                                                                                                               <  2:591911/150‑1 (MQ=255)
                               tATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGc                                                                                                               <  2:1260510/150‑1 (MQ=255)
                                aTCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCg                                                                                                             >  1:848997/1‑151 (MQ=255)
                                aTCATGTACCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCg                                                                                                             >  1:791949/1‑151 (MQ=255)
                                  cATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGAc                                                                                                           >  2:894251/1‑151 (MQ=255)
                                   aTGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACg                                                                                                          >  1:752209/1‑151 (MQ=255)
                                      ttCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGcc                                                                                                        >  1:657498/1‑150 (MQ=255)
                                        ccATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                                     >  2:281770/1‑151 (MQ=255)
                                         cATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTc                                                                                                    <  2:716796/151‑1 (MQ=255)
                                          aTTAACATCATCACTTCGCCAAATACTTAACTGAGATTGGCACGATTGATGTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCCCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTTCGTCGGAAGCCGACGCCa                                                                                                       <  1:1149933/147‑1 (MQ=255)
                                             aaCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATCCGGAc                                                                                                <  1:137282/151‑1 (MQ=255)
                                              aCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACg                                                                                               >  1:1193183/1‑151 (MQ=255)
                                               catcatCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACg                                                                                               >  1:641196/1‑150 (MQ=255)
                                                 tcatcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGAt                                                                                             >  1:346513/1‑150 (MQ=255)
                                                 tcatcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGAt                                                                                             >  1:45375/1‑150 (MQ=255)
                                                      aCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCt                                                                                       <  1:10160/151‑1 (MQ=255)
                                                        ttCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGt                                                                                     >  2:1015766/1‑151 (MQ=255)
                                                        ttCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGt                                                                                     >  1:1079538/1‑151 (MQ=255)
                                                         tCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGt                                                                                     >  2:41703/1‑150 (MQ=255)
                                                          cGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGtat                                                                                   <  1:915383/151‑1 (MQ=255)
                                                          cGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGtat                                                                                   <  1:352512/151‑1 (MQ=255)
                                                            ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaa                                                                                 <  1:505167/151‑1 (MQ=255)
                                                             cAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaa                                                                                <  2:543264/151‑1 (MQ=255)
                                                                 tACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATcc                                                                            >  2:352089/1‑151 (MQ=255)
                                                                  aCTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCt                                                                           >  1:474506/1‑151 (MQ=255)
                                                                  aCTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCt                                                                           >  1:537038/1‑151 (MQ=255)
                                                                     ttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTc                                                                          <  1:33550/149‑1 (MQ=255)
                                                                     ttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCa                                                                         <  1:494817/150‑1 (MQ=255)
                                                                       aCTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTCTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                                        >  2:702659/1‑149 (MQ=255)
                                                                        cTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGc                                                                     >  1:216463/1‑151 (MQ=255)
                                                                          gagaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCt                                                                    >  1:1277667/1‑150 (MQ=255)
                                                                           agaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCtat                                                                  >  1:952134/1‑151 (MQ=255)
                                                                             aTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCtata                                                                 >  1:672305/1‑150 (MQ=255)
                                                                                ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCtata                                                                 >  2:907353/1‑147 (MQ=255)
                                                                                  cACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTg                                                           >  2:26629/1‑151 (MQ=255)
                                                                                     gattgattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATg                                                        <  1:486200/151‑1 (MQ=255)
                                                                                        tgattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGa                                                      <  2:365418/150‑1 (MQ=255)
                                                                                         gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                                     <  2:778987/150‑1 (MQ=255)
                                                                                         gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAAc                                                    >  2:1057497/1‑151 (MQ=255)
                                                                                            ttCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTg                                                  >  1:1242605/1‑150 (MQ=255)
                                                                                              cATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAAc                                               >  2:587425/1‑151 (MQ=255)
                                                                                                 cTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCat                                            >  1:170283/1‑151 (MQ=255)
                                                                                                   ataataAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCatat                                          >  2:217382/1‑151 (MQ=255)
                                                                                                   ataataAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCa                                             <  2:396555/148‑1 (MQ=255)
                                                                                                     aataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATAtt                                         <  2:1277708/150‑1 (MQ=255)
                                                                                                      ataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTg                                        >  1:90930/1‑150 (MQ=255)
                                                                                                       taaCAAAATACTTTCGCTCTGGCTGCTCTTTAGCTCTCTGGCAAAAATTGTAAAATATGCCGTATTTTTGTCGGAAGCAGACGCAATTCCGACGATAGCCATGTATAAAATCCTCATAGCTATAAGATTGATGGTACTGAACCATATTGAc                                      >  1:887089/1‑151 (MQ=255)
                                                                                                       taaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAc                                      >  2:771797/1‑151 (MQ=255)
                                                                                                       taaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAc                                      >  2:1017353/1‑151 (MQ=255)
                                                                                                         aaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                                    >  2:524091/1‑151 (MQ=255)
                                                                                                           aaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGt                                  >  1:806773/1‑151 (MQ=255)
                                                                                                            atatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGtt                                 <  2:504698/151‑1 (MQ=255)
                                                                                                             tataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTg                                >  1:578878/1‑151 (MQ=255)
                                                                                                              ataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGc                               >  1:532799/1‑151 (MQ=255)
                                                                                                                aCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCgg                             <  2:1147980/151‑1 (MQ=255)
                                                                                                                aCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                              <  1:213284/150‑1 (MQ=255)
                                                                                                                aCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                              <  2:536364/150‑1 (MQ=255)
                                                                                                                  tttttCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCgg                             >  2:546860/1‑149 (MQ=255)
                                                                                                                   ttttCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGa                            <  2:51779/149‑1 (MQ=255)
                                                                                                                    tttCTCTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGcc                         >  2:1197209/1‑151 (MQ=255)
                                                                                                                      tctcTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                       >  1:404862/1‑151 (MQ=255)
                                                                                                                      tctcTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                       >  2:93844/1‑151 (MQ=255)
                                                                                                                      tctcTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                       >  2:470224/1‑151 (MQ=255)
                                                                                                                      tctcTGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                       >  2:281947/1‑151 (MQ=255)
                                                                                                                       ctctGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCttt                      >  2:475060/1‑151 (MQ=255)
                                                                                                                         ctGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCttt                      <  2:243317/149‑1 (MQ=255)
                                                                                                                         ctGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCtt                       <  2:663530/148‑1 (MQ=255)
                                                                                                                          tGGCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtct                   >  2:34566/1‑151 (MQ=255)
                                                                                                                           ggCTGCTCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtc                    <  2:115313/149‑1 (MQ=255)
                                                                                                                            gctgctCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTcc                 <  2:430261/151‑1 (MQ=255)
                                                                                                                              tgctCTTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTg               >  2:1197524/1‑151 (MQ=255)
                                                                                                                                  ctTTAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTc             >  2:851801/1‑149 (MQ=255)
                                                                                                                                   tttAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGc          >  1:1271908/1‑151 (MQ=255)
                                                                                                                                    ttAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAg           <  1:309735/149‑1 (MQ=255)
                                                                                                                                     tAGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCa         >  2:871783/1‑150 (MQ=255)
                                                                                                                                      aGCTTGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCaa        <  2:1276203/150‑1 (MQ=255)
                                                                                                                                         ttGCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTa     >  1:950770/1‑150 (MQ=255)
                                                                                                                                           gCTGGCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGc  >  2:711589/1‑151 (MQ=255)
                                                                                                                                             tGACAAAAATTGTAAAATATGCCGGCTTTATGTCGGAAGCCGACGCTATTCGGACGATAGCCCTGTATAAAAGCCTCATGGCTATAAGATTGATGCAGCTGACCCATAGACACGCGTTGCTGAGCGTTTTGACCTGTCAGCCAGTAAgt  <  1:271593/149‑2 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTTAGGGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGC  >  NC_000913/4578406‑4578695

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