Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,970 Δ1 bp 5.6% coding (785/1269 nt) 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,9700C. 5.6% 389.1 / 15.8 108coding (785/1269 nt)umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base C (50/52);  new base . (3/3);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 108 total aligned reads displayed.

ATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGT  >  NC_000913/1231825‑1232114
                                                                                                                                                  |                                                                                                                                                
aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGAt                                                                                                                                               <  2:258644/149‑1 (MQ=255)
aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGAt                                                                                                                                               <  1:861489/149‑1 (MQ=255)
aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGAt                                                                                                                                               <  4:303710/149‑1 (MQ=255)
  ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGc                                                                                                                                             <  2:1397025/149‑1 (MQ=255)
    gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                           >  3:637538/1‑149 (MQ=255)
     tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCGTC‑ACGGATTATCCGTCGATGCGGc                                                                                                                                          >  1:1172470/1‑149 (MQ=255)
     tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgtt                                                                                                                                          <  3:769573/149‑3 (MQ=255)
     tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGc                                                                                                                                          >  1:42777/1‑149 (MQ=255)
       aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAg                                                                                                                                        >  2:481114/1‑149 (MQ=255)
       aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAg                                                                                                                                        <  2:1172470/149‑1 (MQ=255)
            tttAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCAt                                                                                                                                   <  4:637538/149‑1 (MQ=255)
                  gTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAg                                                                                                                             >  2:51554/1‑149 (MQ=255)
                      tGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTAc                                                                                                                         <  2:42777/149‑1 (MQ=255)
                       gCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACg                                                                                                                        >  1:1261016/1‑149 (MQ=255)
                           gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                                    >  2:774918/1‑149 (MQ=255)
                           gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                                    <  2:1261016/149‑1 (MQ=255)
                           gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                                    >  1:1073767/1‑149 (MQ=255)
                             aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGccc                                                                                                                  <  1:481114/149‑1 (MQ=255)
                                aaCGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTa                                                                                                                              <  4:657648/134‑1 (MQ=255)
                                aaCGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTa                                                                                                                              >  3:657648/1‑134 (MQ=255)
                                           aaCTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaaa                                                                                                    <  3:1467089/149‑1 (MQ=255)
                                             cTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACt                                                                                                  >  1:1422095/1‑149 (MQ=255)
                                                  cggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCa                                                                                             <  1:774918/149‑1 (MQ=255)
                                                   ggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAg                                                                                            >  3:133966/1‑149 (MQ=255)
                                                     cgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCg                                                                                          >  4:651783/1‑149 (MQ=255)
                                                        aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCa                                                                                             >  1:318889/1‑143 (MQ=255)
                                                        aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCa                                                                                             <  2:318889/143‑1 (MQ=255)
                                                        aCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGc                                                                                       >  2:809283/1‑149 (MQ=255)
                                                         cccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTAGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                                      >  3:1235819/1‑149 (MQ=255)
                                                          ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGa                                                                                         <  1:973864/145‑1 (MQ=255)
                                                          ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGa                                                                                         >  2:973864/1‑145 (MQ=255)
                                                          ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                     <  2:1422095/149‑1 (MQ=255)
                                                           cTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATc                                                                                    >  1:114426/1‑149 (MQ=255)
                                                            tGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAc                                                                                                   >  4:281291/1‑133 (MQ=255)
                                                            tGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAAc                                                                                                   <  3:281291/133‑1 (MQ=255)
                                                            tGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCa                                                                                   >  2:14472/1‑149 (MQ=255)
                                                             gTTTACAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCaa                                                                                  <  4:1235819/149‑1 (MQ=255)
                                                                   aaCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                            >  1:1434548/1‑149 (MQ=255)
                                                                   aaCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                            <  1:51554/149‑1 (MQ=255)
                                                                   aaCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                            >  1:694086/1‑149 (MQ=255)
                                                                   aaCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGt                                                                            >  1:1439497/1‑149 (MQ=255)
                                                                    aCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                           <  1:809283/149‑1 (MQ=255)
                                                                     cTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCg                                                                          >  1:157542/1‑149 (MQ=255)
                                                                        gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTa                                                                                                                              <  1:863705/94‑1 (MQ=255)
                                                                        gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTa                                                                                                                              >  2:863705/1‑94 (MQ=255)
                                                                        gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                       >  4:804553/1‑149 (MQ=255)
                                                                        gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                       <  2:1073767/149‑1 (MQ=255)
                                                                        gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                       >  1:347284/1‑149 (MQ=255)
                                                                           gagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                      >  1:647770/1‑147 (MQ=255)
                                                                           gagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                      <  2:647770/147‑1 (MQ=255)
                                                                           gagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAt                                                                    <  2:157542/149‑1 (MQ=255)
                                                                           gagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAt                                                                    <  2:1434548/149‑1 (MQ=255)
                                                                             gTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGc                                                                                           >  3:237038/1‑124 (MQ=255)
                                                                              tttGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAtctc                                                                 >  3:66036/1‑149 (MQ=255)
                                                                                 gCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCAc                                                              <  4:66036/149‑1 (MQ=255)
                                                                                 gCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTAGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGCTTTATCTCCAc                                                              <  3:1282434/149‑1 (MQ=255)
                                                                                   aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGt                                                            <  2:694086/149‑1 (MQ=255)
                                                                                      gacgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTtta                                                         <  3:804553/149‑1 (MQ=255)
                                                                                        cgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatt                                                       <  2:770153/149‑1 (MQ=255)
                                                                                         gaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCa                                                                                                                                    >  4:1152362/1‑71 (MQ=255)
                                                                                           aGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAg                                                    <  3:651783/149‑1 (MQ=255)
                                                                                              aGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGgcggcg                                                                                                          <  2:1454627/92‑1 (MQ=255)
                                                                                               ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACtt                                                                                                 <  2:861058/100‑1 (MQ=255)
                                                                                               ggAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACtt                                                                                                 >  1:861058/1‑100 (MQ=255)
                                                                                                 aaaTTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                              <  2:114426/149‑1 (MQ=255)
                                                                                                      aTCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAtt                                         >  1:308505/1‑149 (MQ=255)
                                                                                                       tCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                                      <  1:1407448/71‑1 (MQ=255)
                                                                                                       tCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                                      <  1:662676/71‑1 (MQ=255)
                                                                                                       tCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                                      >  2:662676/1‑71 (MQ=255)
                                                                                                       tCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCAttt                                        >  3:442585/1‑149 (MQ=255)
                                                                                                        cTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATAGTCGGTTTATCTCCACGTTTATTAAGACGTCACCGTTTg                                       >  4:1235342/1‑149 (MQ=255)
                                                                                                              ccGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATAtt                                                                              >  3:647909/1‑104 (MQ=255)
                                                                                                              ccGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATAtt                                                                              <  4:647909/104‑1 (MQ=255)
                                                                                                               cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                     <  1:635927/128‑1 (MQ=255)
                                                                                                               cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                     >  2:635927/1‑128 (MQ=255)
                                                                                                                  tCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGa                             <  1:14472/149‑1 (MQ=255)
                                                                                                                   cGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGaa                            <  2:1439497/149‑1 (MQ=255)
                                                                                                                            aaCGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCg                                                                          <  3:117049/94‑1 (MQ=255)
                                                                                                                            aaCGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCg                                                                          >  4:117049/1‑94 (MQ=255)
                                                                                                                                catcgaCGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                     >  1:326036/6‑111 (MQ=255)
                                                                                                                                catcgaCGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                     <  2:326036/106‑1 (MQ=255)
                                                                                                                                 atcgaCGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                              >  1:1319549/5‑117 (MQ=255)
                                                                                                                                 atcgaCGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                              <  2:1319549/113‑1 (MQ=255)
                                                                                                                                 atcgaCGGATAATCCGT‑GATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCaa              >  4:856510/5‑149 (MQ=255)
                                                                                                                                  tC‑ACGGATTATCTGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                           >  2:14532/1‑135 (MQ=255)
                                                                                                                                  tC‑ACGGATTATCTGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                           <  1:14532/135‑1 (MQ=255)
                                                                                                                                  tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGc                                                                                       <  3:374254/75‑1 (MQ=255)
                                                                                                                                  tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGc                                                                                       >  4:374254/1‑75 (MQ=255)
                                                                                                                                  tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                     <  4:208392/109‑1 (MQ=255)
                                                                                                                                  tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                     >  3:208392/1‑109 (MQ=255)
                                                                                                                                   c‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATa            <  3:1235342/149‑1 (MQ=255)
                                                                                                                                     aCGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCa                                           <  4:396910/117‑1 (MQ=255)
                                                                                                                                     aCGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCa                                           >  3:396910/1‑117 (MQ=255)
                                                                                                                                     aCGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAg           >  1:1263603/1‑149 (MQ=255)
                                                                                                                                     aCGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAg           >  2:223448/1‑149 (MQ=255)
                                                                                                                                      cGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAgc          <  2:308505/149‑1 (MQ=255)
                                                                                                                                             tCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                           >  3:897491/1‑125 (MQ=255)
                                                                                                                                              ccGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGt  >  4:872098/1‑149 (MQ=255)
                                                                                                                                                gTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATc                                                                                    >  2:1092357/1‑65 (MQ=255)
                                                                                                                                                gTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATc                                                                                    <  1:1092357/65‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTCGGT  >  NC_000913/1231825‑1232114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑