Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,578,192 G→T 5.6% Q383K (CAA→AAA)  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,1920GT 5.6% 351.1 / 4.4 120Q383K (CAA→AAA) mcrB5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Reads supporting (aligned to +/- strand):  ref base G (60/52);  new base T (6/2);  total (66/54)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

TCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGT  >  NC_000913/4578044‑4578339
                                                                                                                                                    |                                                                                                                                                   
tCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTTTTCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCCTTTTTTg                                                                                                                                                     <  2:1626079/149‑1 (MQ=255)
 tcggtatcaCAATTTGATTAAGCCATTGGGTATCCGGAGAGGGGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTgg                                                                                                                                                    <  1:116794/144‑1 (MQ=255)
 ccGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTgg                                                                                                                                                    <  2:1460685/149‑1 (MQ=255)
  tgTCTTCACAATTTCATGAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGc                                                                                                                                                   <  2:1810312/148‑1 (MQ=255)
    tcatcaCAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCAt                                                                                                                                                 >  2:1529054/1‑149 (MQ=255)
      atcaCAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATaa                                                                                                                                               >  2:807772/1‑149 (MQ=255)
       tcaCAATTTCATTAAGCCATTGCGTATCCGGAGAGGAGACATATTCCAACCTTCAGCAGATGTAACTCTGCTCCGGCCGGCGTCATCTCCATCGGAATGTGCCCTGTCTGCTTATGTCCTGGTGCAACTCGTTCATTGTTTGTCAtgat                                                                                                                                              >  2:1787484/1‑146 (MQ=255)
        cacaATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAg                                                                                                                                             <  1:987946/149‑1 (MQ=255)
         acaATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAgg                                                                                                                                            >  2:1343244/1‑148 (MQ=255)
         acaATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGa                                                                                                                                            <  2:1115818/149‑1 (MQ=255)
         acaATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCGCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAgg                                                                                                                                            >  1:1357022/1‑148 (MQ=255)
             tttCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCGTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAAACTca                                                                                                                                        >  1:1304964/1‑149 (MQ=255)
              ttCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTcaa                                                                                                                                       <  2:1606535/149‑1 (MQ=255)
               tCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTcaac                                                                                                                                      <  1:1020526/149‑1 (MQ=255)
                      gCCATTGCATATCCGGAGAGGGGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATATAGACTCAACAAATGaa                                                                                                                               <  1:1132282/149‑1 (MQ=255)
                       ccATTGCGTAACCGGAGAGGCGCCATCTTCCAACCCACCGCCGAAGTAGCTATGCCCAACGCGGACTCCCTTCCCAACGATAACGGCCGCCTCCCTGATTTCCTGGTTCCCCTCCTTCATGTTCTGGCATAAAGACTCAACaaaataag                                                                                                                              >  2:314421/1‑144 (MQ=255)
                          ttGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAACTGACGGtt                                                                                                                           >  2:842883/1‑149 (MQ=255)
                                 tCCGGAGAGGTGCCATCTTTCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGTCATAAAGACTCAACAAATGAAGGGTCTGCttt                                                                                                                    >  2:210964/1‑149 (MQ=255)
                                 tCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCttt                                                                                                                    >  1:930182/1‑149 (MQ=255)
                                 tCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGAGAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAAGGAAGGTTCTTCttt                                                                                                                    >  2:1319090/1‑149 (MQ=255)
                                        agGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAGGGTTCTGCTTTTttcttc                                                                                                             >  1:247197/1‑149 (MQ=255)
                                         ggTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCa                                                                                                            <  1:1515092/149‑1 (MQ=255)
                                           tGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGt                                                                                                          >  1:1649924/1‑149 (MQ=255)
                                           ggCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCGTAAAGACTCCACAAATGAAGTTTCTGCTTTTTTATTCAGt                                                                                                          >  2:1524342/2‑149 (MQ=255)
                                             ccATCTTCCAACCCACGGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTaa                                                                                                        >  2:1466220/1‑149 (MQ=255)
                                                tCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGCTTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTaaaaa                                                                                                     >  1:551240/1‑149 (MQ=255)
                                                tCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATGAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCCGTaaaaa                                                                                                     >  2:557492/1‑149 (MQ=255)
                                                 cTTCCAACCCACAGCAGAAGTAACTATGCCCATTGCGGAATCCTTTCCCAAGTAGAGTGGCCTCTTTGCTGAGTGCCGGGTTCATGTTGTTCATTCTTTGGCATAAAGACTCAACAAATTAAGGTTCTGCTTTTTTATTCAGTaaaaaa                                                                                                    <  2:994680/149‑1 (MQ=255)
                                                      aaCCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCACCAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCg                                                                                               >  1:1228243/1‑149 (MQ=255)
                                                          cacaGCAGAAGTAACTATGCCCAAAGCGGAATCCTTTCCCAAGGATAGGGGCCTCTTGGCTGCTTTCCTGGTTCGACTCGTTCCTTTTTTGGACTAAGGACTCCACAAATGAAGGTTCTGCTGTTTTATTCAGTAAAAAATTCCGGAAc                                                                                           >  2:899344/1‑149 (MQ=255)
                                                              gcagAAGTAACTATGCCCCCTGCCGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACtgtg                                                                                       <  1:1768093/149‑1 (MQ=255)
                                                              gcagAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACtgtg                                                                                       <  2:1683418/149‑1 (MQ=255)
                                                               cagAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTgtgg                                                                                      >  1:1731613/1‑149 (MQ=255)
                                                                agAAGTAACTATGCCCAATGCGGTATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAATACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTgtggt                                                                                     <  2:1137777/149‑1 (MQ=255)
                                                                agAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTgtggt                                                                                     <  1:1338980/149‑1 (MQ=255)
                                                                  aaGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCGTAAAGACTCAACAAATGGAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGgtgt                                                                                   >  2:1148501/1‑149 (MQ=255)
                                                                  aaGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTGCATTTTTTGGCATAACGGTCCCACAAATGACGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGgtgt                                                                                   >  1:1500737/1‑149 (MQ=255)
                                                                   aGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTa                                                                                  <  1:619200/149‑1 (MQ=255)
                                                                       aCTATTCCCAATGCTGAATCCTTTCCCAAGTATAGTGGCCTCTTTGCTGATTTCCTGCTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCaa                                                                              <  2:1572949/149‑1 (MQ=255)
                                                                       aCTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCaa                                                                              <  2:152372/149‑1 (MQ=255)
                                                                         tATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCaaaa                                                                            >  1:1427556/1‑149 (MQ=255)
                                                                         tATGCCCAATGCGGAATCCTTTCCCAAGGATAGGGGCCTCTTTGCTGATTTCCTGGTTCCACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGGGCCTCaaac                                                                            >  1:1487980/1‑148 (MQ=255)
                                                                          aTGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAAc                                                                           <  1:445183/149‑1 (MQ=255)
                                                                            gCCCAATGCGCAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCt                                                                         <  1:844805/149‑1 (MQ=255)
                                                                              ccAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTgg                                                                       <  1:1500881/149‑1 (MQ=255)
                                                                               cAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTCGTTTGTTCAGTAAAAAATTCCGGAACTGTGGTGTATCAACACCTGGc                                                                      >  1:141321/1‑149 (MQ=255)
                                                                               cAATGCGGAATCCTTTCCCAAGGATAGTGGCCGCTTTGCTGATTGCCTGGTTCAACTCGATCATTTTTTTTTTTTTTTTCTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAAATTCGGAACTGTGGTGTATCAAACCCTgtc                                                                      >  2:1335307/1‑147 (MQ=255)
                                                                                   gCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTAAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCaa                                                                  <  1:430040/149‑1 (MQ=255)
                                                                                      gAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAatat                                                               <  2:1419710/149‑1 (MQ=255)
                                                                                       aaTCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATc                                                              <  1:704079/149‑1 (MQ=255)
                                                                                        aTCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAATACCCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCt                                                             <  1:1211375/149‑1 (MQ=255)
                                                                                         tCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTa                                                            >  1:1042702/1‑149 (MQ=255)
                                                                                            tttCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTTTTTTTTTTTTTCAACAAATGCAGGTTCTGATTTTTTATTCAGTAAATGATTCCGGAACTGTGGTGTTTCAAATCGTTGCCCAATACCCATgg                                                         >  2:1678807/1‑148 (MQ=255)
                                                                                            tttCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAACTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAAGATCTATga                                                         >  2:466982/1‑149 (MQ=255)
                                                                                            tttCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATgg                                                         >  1:462610/1‑148 (MQ=255)
                                                                                            tgtCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATga                                                         >  1:1071456/3‑149 (MQ=255)
                                                                                            attCCCAAGGATAGTGGCCTTTTTGCTGATTTCCTGGGTCAACTCGTTCATTCTTTGGCATAAAGACTCAACAAATGAAGGGTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTAGCAAAACCTGGCTCAATATCTATga                                                         <  2:764932/148‑1 (MQ=255)
                                                                                             ttCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTGGCATGAAGACTCAACAAATGAAGGTTCTGATTTTTTATTCAGTAAAAAATTCCGGAACTGGGGTGTTTCAAAACCTCGCTCAAAATCTATgaaa                                                       >  2:604518/1‑149 (MQ=255)
                                                                                              tCCCAAGGAAAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATgaaa                                                       <  1:589452/149‑1 (MQ=255)
                                                                                                  aaGGATAGTGTCCTCTTTGCTGATTGCCTGGTTCATCTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATgaaagaaa                                                   <  1:337325/149‑1 (MQ=255)
                                                                                                  aaGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAATAATTGCGGAACTGTGGTGTATCAAAACCTGGATCAATAACTATgaaagaaa                                                   >  2:1537619/1‑149 (MQ=255)
                                                                                                  aaGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATgaaagaaa                                                   >  2:1232083/1‑149 (MQ=255)
                                                                                                    ggATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTTTTTTTATTTCTCAACAAATGAAGGCTCTGCTTTTTTATTCAGTGAAAAATGCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCCATTAAAGCAAAt                                                 >  2:269770/1‑149 (MQ=255)
                                                                                                    ggATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTAGGAAAGAAAAt                                                 >  2:691661/1‑149 (MQ=255)
                                                                                                        aGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCACGATCTATGACAGAAAATCgcc                                             >  1:1750083/1‑147 (MQ=255)
                                                                                                        aGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGAATAAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAGCCTGGCTCCATATCTATGAAAGAAAATGGtg                                             >  1:1750092/1‑148 (MQ=255)
                                                                                                        aGTGGCCTCTTTGCTGATTCCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCCGTAGTGGCTTCTTGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTc                                             >  1:1318860/1‑149 (MQ=255)
                                                                                                           ggCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGc                                          <  1:1386712/149‑1 (MQ=255)
                                                                                                              ctctTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGTCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAAACTGTCTCAATATCTATGAAAGAAAATCGTCTGCGTa                                       >  2:1021972/1‑149 (MQ=255)
                                                                                                               tctTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTTTTTTTTTTTTTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAATAAGTCCGGAACTGTGGTGTTTCCAAACCTTGCTCAATATCTATGAAAGCAAATCGTCTGCGTAg                                      >  2:1623863/1‑149 (MQ=255)
                                                                                                                ctTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAgg                                     >  2:809562/1‑149 (MQ=255)
                                                                                                                 tttGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTACAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAggg                                    >  1:726417/1‑149 (MQ=255)
                                                                                                                 tttGATGATTTCCTCGTTCATCTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAggg                                    <  2:453999/149‑1 (MQ=255)
                                                                                                                  ttGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGc                                   <  1:693835/149‑1 (MQ=255)
                                                                                                                   tGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAGGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGACCTGTGGTGCACCAGAACCTGGCGCAATATCTATGAAAGAAAATCGTCTTGGTAGGGCa                                  >  1:592279/1‑149 (MQ=255)
                                                                                                                    gCTGATTTCCTGGTTCAACTCGTTCATTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATa                                <  1:296125/149‑1 (MQ=255)
                                                                                                                     cttaTTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGTCTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATa                                <  2:250325/146‑1 (MQ=255)
                                                                                                                     cTGGTTTCCTGGTTCAACTCGTTCATTTTTTGGCGTAAAGACCCAACAAATGAAGGTTCTGCTTTTCTATTGAGTAAACAAATCCGGACCTGTGGTGTATCAAATCCTGGCCCAATATCTATGAAAGAAAATCGTCTGCGTAGAGCATa                                <  2:615633/149‑1 (MQ=255)
                                                                                                                     cTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATa                                <  2:679307/149‑1 (MQ=255)
                                                                                                                     cTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATa                                <  1:807772/149‑1 (MQ=255)
                                                                                                                     cTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTCTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATa                                <  1:795023/149‑1 (MQ=255)
                                                                                                                      tGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAg                               <  2:1104698/149‑1 (MQ=255)
                                                                                                                      tGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAAACATCTCCGTAGCGCATAg                               >  2:170158/1‑149 (MQ=255)
                                                                                                                       gagttCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCTTTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTAGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGCGCGTAGt                              >  2:1020395/4‑149 (MQ=255)
                                                                                                                               ggTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCTTCTGCGTAGGGCATAGTCAACAACg                      <  1:1011469/149‑1 (MQ=255)
                                                                                                                                gTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTCAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATAGCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACgg                     >  2:45312/1‑149 (MQ=255)
                                                                                                                                 ttCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGc                    <  1:616476/149‑1 (MQ=255)
                                                                                                                                   cAACGCGTTCATTTTTTTGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGcaa                  <  2:1381614/149‑3 (MQ=255)
                                                                                                                                     aCTCGTTCATTTTTTGGCATAAAGACTCGGCAAATGAAGGTTCTGCTTTTTTATTCAGTAAAACATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGGAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCaga                >  1:1178730/1‑149 (MQ=255)
                                                                                                                                     aCTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCaga                >  1:749510/1‑149 (MQ=255)
                                                                                                                                     aCTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCaga                >  1:47097/1‑149 (MQ=255)
                                                                                                                                          ttCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCg           >  1:1432288/1‑149 (MQ=255)
                                                                                                                                          ttCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCAGGAAATGTGGAGTATCAAACCCTGGCTCAATATCTATGACAGAAAATCGTCTGAGTAGGGCAAAGTCGACAACGGCCagagagga           >  2:863479/1‑147 (MQ=255)
                                                                                                                                            cATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTAAAAACCGGCCAGAGAGCGAt         >  1:683037/1‑149 (MQ=255)
                                                                                                                                             gttttttGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATc        <  2:930182/148‑1 (MQ=255)
                                                                                                                                             aTTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTTTTTTTTTTTAAATTCCGGAACTGTGGTGTATTAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTATGCCTTATTTTACAACGGCCAGAGAGCGATc        >  2:443176/1‑149 (MQ=255)
                                                                                                                                             aTTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATc        <  1:280116/149‑1 (MQ=255)
                                                                                                                                             aTTTCTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATc        <  1:1044899/149‑1 (MQ=255)
                                                                                                                                               tttttGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCgg      >  2:1467773/1‑149 (MQ=255)
                                                                                                                                                   tGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGt  <  1:949276/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCCGTCATCACAATTTCATTAAGCCATTGCGTATCCGGAGAGGTGCCATCTTCCAACCCACAGCAGAAGTAACTATGCCCAATGCGGAATCCTTTCCCAAGGATAGTGGCCTCTTTGCTGATTTCCTGGTTCAACTCGTTCATTTTTTGGCATAAAGACTCAACAAATGAAGGTTCTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGT  >  NC_000913/4578044‑4578339

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