Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,578,548 G→A Q264* (CAG→TAG)  mcrB ← 5‑methylcytosine‑specific restriction enzyme subunit McrB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,578,5480GA99.1% 785.9 / ‑3.1 230Q264* (CAG→TAG) mcrB5‑methylcytosine‑specific restriction enzyme subunit McrB
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (121/107);  total (122/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 231 total aligned reads displayed.

GGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTC  >  NC_000913/4578411‑4578684
                                                                                                                                         |                                                                                                                                        
gggAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAAc                                                                                                                                        >  2:2195821/1‑140 (MQ=255)
    aCAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaaa                                                                                                                                    <  2:320427/140‑1 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATCCTTTTTCTCTGCCTACTCTTTCGCTTGCTAACaaaaa                                                                                                                                   >  1:914525/1‑140 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaaaa                                                                                                                                   >  2:129621/1‑140 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaaaa                                                                                                                                   <  2:2489451/140‑1 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaaaa                                                                                                                                   >  1:2302921/1‑140 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaaaa                                                                                                                                   >  1:447153/1‑140 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaaa                                                                                                                                    >  2:2475063/1‑139 (MQ=255)
     cAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACaaa                                                                                                                                     >  1:93132/1‑138 (MQ=255)
        aCCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATCAAAATATACTTTCTCTCTGGCTGCTCTTTACCTTGCTAACAAAAATTg                                                                                                                                >  2:679268/1‑140 (MQ=255)
           aGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTaa                                                                                                                             <  2:35944/140‑1 (MQ=255)
                  aCCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGGTGCTCTTCAGCGTGCTAACAAAAATTGTAAAATATGc                                                                                                                      >  1:808113/1‑140 (MQ=255)
                  aCCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGc                                                                                                                      >  1:1379146/1‑140 (MQ=255)
                       gTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTTTTTAGCTTGCTAACAAAAATTTTAAAATATGCCGTCt                                                                                                                 >  1:845532/1‑140 (MQ=255)
                       gTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCt                                                                                                                 >  1:263459/1‑140 (MQ=255)
                       gTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCt                                                                                                                 >  2:2139243/1‑140 (MQ=255)
                                 ttCCATTAACATCGTCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGa                                                                                                       >  1:1519574/1‑140 (MQ=255)
                                 ttCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTATAATATGCCGTCTTTACGTCGGa                                                                                                       >  1:1011985/1‑140 (MQ=255)
                                 ttCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGcc                                                                                                                     <  1:2482818/126‑1 (MQ=255)
                                 ttCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGcc                                                                                                                     >  2:2482818/1‑126 (MQ=255)
                                 ttCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCg                                                                                                         >  1:1986972/1‑138 (MQ=255)
                                  tCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGACTGCTCTTTAGCTTGCGAACAAAATTTGGAAAATATGCTGTCTTTACGTctgtc                                                                                                      >  1:763081/1‑136 (MQ=255)
                                      ttAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGt                                                                                                           <  2:2256135/131‑1 (MQ=255)
                                      ttAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGt                                                                                                           >  1:2256135/1‑131 (MQ=255)
                                      ttAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCg                                                                                                  >  1:1871339/1‑140 (MQ=255)
                                      ttAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCg                                                                                                  >  2:622333/1‑140 (MQ=255)
                                       tAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCTCGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGa                                                                                                 <  2:828427/140‑1 (MQ=255)
                                         aCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTCCGTCGGAAGCCGACg                                                                                               <  1:522169/140‑1 (MQ=255)
                                         aCATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACg                                                                                               >  1:1437194/1‑140 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          >  2:1485051/1‑140 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          >  2:386051/1‑140 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          >  2:187035/1‑140 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          <  1:1542488/140‑1 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          >  1:133942/1‑140 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCAGAAGCCGACGCCAtt                                                                                          <  2:447153/140‑1 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGAAAAATATCCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          >  1:1194140/1‑140 (MQ=255)
                                              atcaCTTCGCCAAATACTTTACCGAGATTGGTAGGTTTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          <  2:1132457/140‑1 (MQ=255)
                                                caCTTCGCCAAATACTTTACTGTGATTGGCACGATGGATTTCATCTATTATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCg                                                                                        <  2:904437/140‑1 (MQ=255)
                                                   ttCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGAc                                                                                     >  1:374439/1‑140 (MQ=255)
                                                     cGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCAtt                                                                                          <  2:78002/133‑1 (MQ=255)
                                                     cGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGa                                                                                   >  2:295244/1‑140 (MQ=255)
                                                       ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGACGTCTTTGCGTTGGAAGCCGAGGCCATTCGGACGATa                                                                                 >  1:675406/1‑140 (MQ=255)
                                                        cAAATACTCTACTGAGTGTGGCACGATTGATCTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAg                                                                                <  2:1194140/140‑1 (MQ=255)
                                                         aaaTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGc                                                                               <  1:622333/140‑1 (MQ=255)
                                                          aaTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                              <  2:411067/140‑1 (MQ=255)
                                                          aaTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                              <  2:1579000/140‑1 (MQ=255)
                                                          aaTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                              <  2:237640/140‑1 (MQ=255)
                                                          aaTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGATGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                              <  1:2373589/140‑1 (MQ=255)
                                                          aaTACTTTACTGAGATTGGCACCATTGATTTAATCCATACTAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                              <  2:738886/140‑1 (MQ=255)
                                                           aTACTTTACTGAGATTGGCACGCTTGATTTCTTCTATAATAAAAATATCCTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGccc                                                                             <  2:926628/140‑1 (MQ=255)
                                                           aTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGccc                                                                             <  1:129621/140‑1 (MQ=255)
                                                           aTACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTAGGAAGCCGACGcc                                                                                             <  2:564104/124‑1 (MQ=255)
                                                               tttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCTACCCCATTCAGACGATAGCCCTGta                                                                         >  1:1054524/1‑140 (MQ=255)
                                                               tttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGta                                                                         >  1:1860853/1‑140 (MQ=255)
                                                                 tACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGtata                                                                       >  1:1330083/1‑140 (MQ=255)
                                                                   cTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaa                                                                     <  1:1688277/140‑1 (MQ=255)
                                                                   cTGAGATTGGCACGACTGATTTCATCTATAATAAAAACATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaa                                                                     <  1:1121754/140‑1 (MQ=255)
                                                                    tGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaaa                                                                    >  2:1519356/1‑140 (MQ=255)
                                                                     gagaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAt                                                                   >  1:1996910/1‑140 (MQ=255)
                                                                      agaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGccc                                                                             <  1:2454652/129‑1 (MQ=255)
                                                                      agaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGc                                                                               >  2:1728343/1‑127 (MQ=255)
                                                                      agaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATc                                                                  >  2:1693352/1‑140 (MQ=255)
                                                                        aTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCt                                                                            <  2:1674698/128‑1 (MQ=255)
                                                                        aTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCt                                                                >  2:2078625/1‑140 (MQ=255)
                                                                        aTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCt                                                                <  2:507522/140‑1 (MQ=255)
                                                                        aTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCt                                                                >  1:1355138/1‑140 (MQ=255)
                                                                           ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             <  2:1379146/140‑1 (MQ=255)
                                                                           ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             <  1:2078610/140‑1 (MQ=255)
                                                                           ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             <  1:393318/140‑1 (MQ=255)
                                                                           ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             <  2:1330083/140‑1 (MQ=255)
                                                                           ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             <  2:2039533/140‑1 (MQ=255)
                                                                           ggCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCAACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             >  1:878263/1‑140 (MQ=255)
                                                                             cACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGccc                                                                             <  2:2346452/122‑1 (MQ=255)
                                                                             cACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCAt                                                             <  2:1993422/138‑1 (MQ=255)
                                                                             cACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGAAGCCATTCGGACCATTGCCCTGGATAAAATCCTCATac                                                           >  1:1177270/1‑139 (MQ=255)
                                                                                    gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCATCATAGCTATAAg                                                    <  1:473125/140‑1 (MQ=255)
                                                                                     attTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGa                                                   <  1:2180245/140‑1 (MQ=255)
                                                                                      tttCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGAt                                                  <  1:835109/140‑1 (MQ=255)
                                                                                       ttCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                 >  2:1403959/1‑140 (MQ=255)
                                                                                       ttCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                 >  1:7074/1‑140 (MQ=255)
                                                                                         cATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGa                                               <  2:133942/140‑1 (MQ=255)
                                                                                               taataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAACCCGGTATAAAATCCTCATAGGTATAAGATTGATGGAAc                                         >  2:934074/1‑140 (MQ=255)
                                                                                                aataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACt                                        >  2:1999975/1‑140 (MQ=255)
                                                                                                aataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACt                                        <  1:285081/140‑1 (MQ=255)
                                                                                                 ataaAAATATACTTTTTTTCTGTCTGCTCTTTAGCGTGCTAACAAAAACTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTg                                       <  2:845765/140‑1 (MQ=255)
                                                                                                  taaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGa                                      <  1:1403959/140‑1 (MQ=255)
                                                                                                   aaaaaTATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGaa                                     <  2:1362430/140‑1 (MQ=255)
                                                                                                   aaaaaTATACTTTTTCTCTGGCAGCTCTTTACCTTGCTATCAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGaa                                     <  2:1027520/140‑1 (MQ=255)
                                                                                                         ataCTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCatat                               >  1:2122867/1‑140 (MQ=255)
                                                                                                          tgcTTTTTCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCTATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATAtt                              >  1:2034415/3‑140 (MQ=255)
                                                                                                               tttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTg                                       >  2:1949075/1‑126 (MQ=255)
                                                                                                               tttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTg                                       <  1:1949075/126‑1 (MQ=255)
                                                                                                               tttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgc                         >  2:2439959/1‑140 (MQ=255)
                                                                                                                ttCTCTGGCTGCTCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                        >  2:836686/1‑140 (MQ=255)
                                                                                                                       gctgctCTTTAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGa                 <  1:1622437/140‑1 (MQ=255)
                                                                                                                               ttAGCTTGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtc         >  2:1246699/1‑140 (MQ=255)
                                                                                                                                  gCTTGCTAACAAAAATCATAAATTATGCCGTCTATCCGTAGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTcc      <  2:1120761/140‑1 (MQ=255)
                                                                                                                                   cTTGCTAACAAAAATTGTAAAATATGCCGTATTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt     >  1:243735/1‑140 (MQ=255)
                                                                                                                                     tGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGt   <  2:280201/140‑1 (MQ=255)
                                                                                                                                     tGCTAACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACAGAACCATCTTGACGTGTTGCGGAGCCTTTTCTCCTGt   >  1:1001937/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
GGGAACAGACCAGTTTTCACCTCGTTTATCATGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTC  >  NC_000913/4578411‑4578684

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