Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,963 T→A 5.4% Y260N (TAT→AAT)  umuC → translesion error‑prone DNA polymerase V subunit; DNA polymerase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,231,9630TA 5.4% 318.9 / 4.8 106Y260N (TAT→AAT) umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base T (56/44);  new base A (3/3);  total (59/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCG  >  NC_000913/1231817‑1232112
                                                                                                                                                   |                                                                                                                                                     
tCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATc                                                                                                                                                    <  1:402070/150‑1 (MQ=255)
    gTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTc                                                                                                                                                <  1:428551/150‑1 (MQ=255)
       tATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACACTGTTTGCAACTGGAAGAGTTTGTACCGACGAAGCAGGAAATTATGTGTTCCGGCTCGTTTGGTGAACGCATC‑ACGGATTATCCCTCGAt                                                                                                                                             >  2:465014/1‑150 (MQ=255)
        aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGACGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATg                                                                                                                                            <  2:580249/150‑1 (MQ=255)
         tCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGc                                                                                                                                           >  2:430368/1‑150 (MQ=255)
          ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCg                                                                                                                                          >  2:1058468/1‑150 (MQ=255)
           cGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                         >  1:310744/1‑150 (MQ=255)
           cGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                         >  1:110917/1‑150 (MQ=255)
           cGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACACTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                         >  1:170161/1‑150 (MQ=255)
            gTAAACATTTTAATGTCGTTCGAGAAAGAACGGTGAGTGAACTGCGCGGCGAACCCTGTTTGCACCTGGAAGAGTTTGCACCGACGACGCAGGAAATTATCTGTTACCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGc                                                                                                                                        <  2:1139228/150‑1 (MQ=255)
            gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGc                                                                                                                                        <  1:810184/150‑1 (MQ=255)
            gTAAACATTTTAATGTCGTGCTCGAAAGAACGGAGCGTGAACTGCGCGGCGAACCCTGTATGCAACTGGAAGAGTTTGCACGGAGGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCCTC‑ACGGAGTATCCGTCGATGCGGc                                                                                                                                        >  1:336590/1‑150 (MQ=255)
             tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGTCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCa                                                                                                                                       >  1:219927/1‑150 (MQ=255)
             tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgttc                                                                                                                                       <  1:662624/150‑4 (MQ=255)
               aaCATTTTAATGTCGTGCACGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgttcac                                                                                                                                     <  2:164138/150‑6 (MQ=255)
                 cATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGcc                                                                                                                                   >  1:900744/1‑150 (MQ=255)
                 cATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGcc                                                                                                                                   >  1:1231992/1‑150 (MQ=255)
                   ttttAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCAt                                                                                                                                 >  1:458334/1‑150 (MQ=255)
                        aTGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTa                                                                                                                            <  1:1731121/150‑1 (MQ=255)
                          gTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGt                                                                                                                          <  1:747663/150‑1 (MQ=255)
                            cGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                                        <  2:946720/150‑1 (MQ=255)
                            cGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                                        <  1:1768408/150‑1 (MQ=255)
                               gCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGc                                                                                                                     >  1:840901/1‑150 (MQ=255)
                                cTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGATATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                                    >  2:590103/1‑150 (MQ=255)
                                cTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                                    <  2:1094878/150‑1 (MQ=255)
                                  cGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                                  >  1:1387790/1‑150 (MQ=255)
                                     aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCg                                                                                                               <  1:1119144/150‑1 (MQ=255)
                                         aCGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATAATCTGTTCCCGCTCGTGTGGTGAACGCATC‑ACGGATTATCCGGCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcg                                                                                                           <  1:913145/150‑1 (MQ=255)
                                          cGGTGCGTGAACTGCGCGGCGCACCCTGTTTGCAACTGTACGAGTTTGCACCGACGAAGCAGGAAATTATCTGATCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcgg                                                                                                          <  1:338346/150‑1 (MQ=255)
                                          cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcgg                                                                                                          >  1:1712790/1‑150 (MQ=255)
                                               cGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaa                                                                                                     >  2:1578547/1‑150 (MQ=255)
                                               cGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGAATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaa                                                                                                     >  2:225555/1‑150 (MQ=255)
                                                gTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGTATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaa                                                                                                    <  2:1186385/150‑1 (MQ=255)
                                                 tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaa                                                                                                   >  1:1626382/1‑150 (MQ=255)
                                                 tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaa                                                                                                   <  2:151023/150‑1 (MQ=255)
                                                   aaCTGCGCGGCGAACCCTGTTTGCTACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAc                                                                                                 <  2:788796/150‑1 (MQ=255)
                                                    acggcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACt                                                                                                <  1:1724907/147‑1 (MQ=255)
                                                    aCTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACt                                                                                                <  1:1392326/150‑1 (MQ=255)
                                                      tGCGCGGCGAACCCTGTATGCAACAGGAAGAGTTAGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTc                                                                                              <  2:290780/150‑1 (MQ=255)
                                                       gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCg                                                                                             <  2:744805/150‑1 (MQ=255)
                                                       gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATAATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCg                                                                                             <  1:1445213/150‑1 (MQ=255)
                                                        cgcgGCGAACCCTGTTTGCAACTGGCAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGAAGTTACGCTGCCCGTGCGGCGGAAAAACTTCGc                                                                                            <  1:365587/150‑1 (MQ=255)
                                                        cgcgGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGc                                                                                            <  2:1520747/150‑1 (MQ=255)
                                                                   cTGTTAGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATAATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCa                                                                                 <  2:601671/150‑1 (MQ=255)
                                                                    agTTTGCAACTGGAAGAGTTTGCACCGACGAAGAAGGAAAGTATCTGTTCCCGCGCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGAAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTGCGCAGCGAGCATCaa                                                                                >  2:456416/2‑150 (MQ=255)
                                                                     gTTTGCACCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAat                                                                               <  2:750068/150‑1 (MQ=255)
                                                                     gTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAat                                                                               <  2:756600/150‑1 (MQ=255)
                                                                      tttGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAata                                                                              <  2:219927/150‑1 (MQ=255)
                                                                         gCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACTCATC‑ACGGATTATCCGTCGATTCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTg                                                                           <  1:341158/150‑1 (MQ=255)
                                                                         gCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCAGTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTg                                                                           >  1:1401547/1‑150 (MQ=255)
                                                                          cAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                          >  1:62926/1‑150 (MQ=255)
                                                                           aaCTGGAAGATTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                         <  1:1365498/150‑1 (MQ=255)
                                                                           aaCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCACAGCGAGCATCAATATTGTc                                                                         <  2:848658/150‑1 (MQ=255)
                                                                               ggAATAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGTCCTTTTGTAGTTAAGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCAGCAATATTGTCGGtt                                                                     >  1:270620/1‑150 (MQ=255)
                                                                               ggAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                     >  1:710785/1‑150 (MQ=255)
                                                                               ggAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                     >  2:1729295/1‑150 (MQ=255)
                                                                                gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                    >  2:1848419/1‑150 (MQ=255)
                                                                                gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                    >  2:813328/1‑150 (MQ=255)
                                                                                  agagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAt                                                                  >  1:386611/1‑150 (MQ=255)
                                                                                   gagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAtc                                                                 <  2:110917/150‑1 (MQ=255)
                                                                                     gTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGTTCGCATC‑ACGGATTATCGGTAGATGCGGCAGGCCATTTGTAGTTACGCTGCCGGGGCGGCGGAAATCCTTCGCAGCGAGCATCAATATTGTCGGttttagtc                                                               >  2:253138/1‑145 (MQ=255)
                                                                                     gTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAtctc                                                               >  1:59842/1‑150 (MQ=255)
                                                                                     gTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTAGTAGTTACGCTGCCCGGGCGGCGGAAAACCTTCGCAGCGAGCATCAAGATTGTCGGTTTAgctc                                                               >  2:278591/1‑150 (MQ=255)
                                                                                           aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGtt                                                         >  2:1557964/1‑150 (MQ=255)
                                                                                             cgacgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTtta                                                       <  2:1648502/150‑1 (MQ=255)
                                                                                                 gaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                   <  2:577796/150‑1 (MQ=255)
                                                                                                  caGCAGGAAATTATCTGCTCCCGCTCCGTGGGTGAACGCATC‑ACGGAGTATCCGTCGATGCGGCAGGCCATATGTAGGTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGGTTATCTCCACGTTTATTAAg                                                  <  2:850280/149‑1 (MQ=255)
                                                                                                  aaGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAg                                                  <  1:1576894/150‑1 (MQ=255)
                                                                                                      aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGt                                              >  1:838728/1‑150 (MQ=255)
                                                                                                      aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGt                                              >  2:1536183/1‑150 (MQ=255)
                                                                                                       ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                             >  2:640940/1‑150 (MQ=255)
                                                                                                         aaaTTATCTGTTCCCGCTCGTTTTGTGAACGCATC‑AAGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCAc                                           <  2:448321/150‑1 (MQ=255)
                                                                                                             tATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCCCCAtt                                       >  2:1107449/1‑150 (MQ=255)
                                                                                                                  gTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTgcgc                                  <  2:63892/150‑1 (MQ=255)
                                                                                                                   ttCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCt                                 <  1:1508063/150‑1 (MQ=255)
                                                                                                                       cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                             >  1:359278/1‑150 (MQ=255)
                                                                                                                        gCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATg                            <  2:1327335/150‑1 (MQ=255)
                                                                                                                            gTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAcc                        >  1:421287/1‑150 (MQ=255)
                                                                                                                               tGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCtta                     >  2:394109/1‑150 (MQ=255)
                                                                                                                               tGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCtta                     <  1:867045/150‑1 (MQ=255)
                                                                                                                               tGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCtta                     >  1:1067317/1‑150 (MQ=255)
                                                                                                                                ggTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttat                    >  2:1588806/1‑150 (MQ=255)
                                                                                                                                  tGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCttatta                  >  2:1785092/1‑150 (MQ=255)
                                                                                                                                    gccGCATCGACGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACg                >  2:1622705/3‑150 (MQ=255)
                                                                                                                                     ccGCATCGACGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACgg               >  2:1604159/2‑150 (MQ=255)
                                                                                                                                      cGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc              <  1:813328/150‑1 (MQ=255)
                                                                                                                                      cGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGc              >  2:1864015/1‑150 (MQ=255)
                                                                                                                                       gcatcgaCGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTAGTAAGACGTCCCCATTTGCGCTCAATGACCCTTATTACGGCa             >  2:742904/7‑150 (MQ=255)
                                                                                                                                       gCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCa             >  2:48813/1‑150 (MQ=255)
                                                                                                                                        cATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCaa            >  2:1414281/1‑150 (MQ=255)
                                                                                                                                          tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATa          >  1:774094/1‑150 (MQ=255)
                                                                                                                                           c‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTCCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAg         >  1:1437609/1‑150 (MQ=255)
                                                                                                                                              cGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAgcg       <  1:1025482/150‑1 (MQ=255)
                                                                                                                                              cGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAgcg       <  1:1028576/150‑1 (MQ=255)
                                                                                                                                               ggATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCTGTTTATCTCCACGTTTATTAAGACGTCACCATATGCGCTCAATGAACCTTATTACGGCAATAgcgc      >  2:1670066/1‑150 (MQ=255)
                                                                                                                                               ggATTATCCGTCGATGCGGCAGGCCATATGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTAATTAAGACGTCACCATTTGCGCTCAAAGAACCTTATTACGGCAATAgcgc      >  2:1662581/1‑150 (MQ=255)
                                                                                                                                                gATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGcgcg     >  1:39515/1‑150 (MQ=255)
                                                                                                                                                  ttATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTc   <  2:458334/150‑1 (MQ=255)
                                                                                                                                                  ttATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTc   <  2:170161/150‑1 (MQ=255)
                                                                                                                                                   tATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCg  <  2:1787797/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCG  >  NC_000913/1231817‑1232112

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