Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,231,963 T→A 5.2% 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.2% 332.8 / 4.4 96Y260N (TAT→AAT) umuCtranslesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base T (48/43);  new base A (2/3);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTC  >  NC_000913/1231817‑1232111
                                                                                                                                                   |                                                                                                                                                    
tCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTAt                                                                                                                                                    >  3:116843/1‑149 (MQ=255)
    gTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAAGTGCGCGGCGAAAATTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGt                                                                                                                                                <  2:504999/149‑1 (MQ=255)
    gTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGt                                                                                                                                                <  2:854643/149‑1 (MQ=255)
    gTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGt                                                                                                                                                <  3:614024/149‑1 (MQ=255)
     tttATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTc                                                                                                                                               <  4:322888/149‑1 (MQ=255)
       tATCCGTAAACAATTTAATGTCGTGCTCGAAAGAACGGTGCGTCAACTGCGCGGCGAAACCCGTTTGCAACTCGAAGTTTTTGCACCGACTAAGCAGGAAATTATCTGTTCCCGCTCGTCTGGTGAACGCATCGACGGATAATCCGtga                                                                                                                                              <  4:581213/149‑3 (MQ=255)
        aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGAt                                                                                                                                            <  3:706881/149‑1 (MQ=255)
        aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGAt                                                                                                                                            >  4:32991/1‑149 (MQ=255)
         tCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATg                                                                                                                                           >  3:275708/1‑149 (MQ=255)
           cGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTCTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATAATCCGT‑GATGCgt                                                                                                                                        <  4:755816/149‑2 (MQ=255)
            gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                        >  3:3590/1‑149 (MQ=255)
            gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                        >  4:782343/1‑149 (MQ=255)
            gTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGATGAACGCATC‑ACGGATTATCCGTCGATGCgg                                                                                                                                        >  2:21362/1‑149 (MQ=255)
             tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgtt                                                                                                                                       >  1:134638/1‑147 (MQ=255)
              aaaCATTTTAATGTCGTGCTCGAAAGAGCGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCa                                                                                                                                      <  1:380510/149‑1 (MQ=255)
              aaaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCGACGGATAATCCGT‑GATGCgttc                                                                                                                                      <  2:134638/149‑4 (MQ=255)
               aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAg                                                                                                                                     >  1:115350/1‑149 (MQ=255)
               aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAg                                                                                                                                     >  4:27829/1‑149 (MQ=255)
                          gTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAg                                                                                                                          >  1:27374/1‑149 (MQ=255)
                           tCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGt                                                                                                                         >  4:113231/1‑149 (MQ=255)
                             gTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTa                                                                                                                       >  2:414482/1‑149 (MQ=255)
                              tGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTAc                                                                                                                      <  4:401664/149‑1 (MQ=255)
                                 tCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt                                                                                                                   >  4:300356/1‑149 (MQ=255)
                                   gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                                 <  1:564587/149‑1 (MQ=255)
                                   gAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGc                                                                                                                 >  3:698260/1‑149 (MQ=255)
                                    aaaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGcc                                                                                                                <  3:395778/149‑1 (MQ=255)
                                    aaaGAACGGTGCGAGAACTGCGCGGCGAACCCTGTTTGCAACTGGGAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGACGCGGCAGGCCATTTGTAGTTACGCTGcc                                                                                                                >  2:802297/1‑149 (MQ=255)
                                     aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGccc                                                                                                               >  3:762040/1‑149 (MQ=255)
                                          cGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcg                                                                                                          >  3:124570/1‑149 (MQ=255)
                                            gTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcggc                                                                                                        <  4:714498/149‑1 (MQ=255)
                                            gTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGggcggc                                                                                                        <  2:115350/149‑1 (MQ=255)
                                                 tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCggg                                                                                                            >  3:85217/1‑140 (MQ=255)
                                                 tGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCggg                                                                                                            <  4:85217/140‑1 (MQ=255)
                                                   aaCTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaaa                                                                                                 >  1:478100/1‑149 (MQ=255)
                                                   aaCTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGaaaaa                                                                                                 >  1:109856/1‑149 (MQ=255)
                                                       gcgcGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTc                                                                                             <  3:238855/149‑1 (MQ=255)
                                                         gcggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGc                                                                                           >  1:119293/1‑149 (MQ=255)
                                                          cggcgAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCa                                                                                          <  4:275708/149‑1 (MQ=255)
                                                               aaCCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGAAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAg                                                                                     <  3:142114/149‑1 (MQ=255)
                                                                 cccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                                   <  2:564553/149‑1 (MQ=255)
                                                                 cccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCa                                                                                   >  1:265408/1‑149 (MQ=255)
                                                                  ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                  >  2:691181/1‑149 (MQ=255)
                                                                  ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                  >  4:85017/1‑149 (MQ=255)
                                                                  ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                  >  3:110274/1‑149 (MQ=255)
                                                                  ccTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCAt                                                                                  >  4:25641/1‑149 (MQ=255)
                                                                    tGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCa                                                                                >  3:110083/1‑149 (MQ=255)
                                                                      tttGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAat                                                                              <  2:109856/149‑1 (MQ=255)
                                                                         gCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATAtt                                                                           <  3:236158/149‑1 (MQ=255)
                                                                            aCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                        <  3:32991/149‑1 (MQ=255)
                                                                            aCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                        >  2:42292/1‑149 (MQ=255)
                                                                            aCTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                        >  4:853092/1‑149 (MQ=255)
                                                                                gtagagTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCAGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGAAGCGAGCATCAATATTGTCgggt                                                                    >  4:270472/3‑147 (MQ=255)
                                                                                gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                    <  3:658790/149‑1 (MQ=255)
                                                                                gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                    <  2:314711/149‑1 (MQ=255)
                                                                                gAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGtt                                                                    >  2:581932/1‑149 (MQ=255)
                                                                                 aaGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                   <  3:27829/149‑1 (MQ=255)
                                                                                 aaGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGttt                                                                   <  2:43101/149‑1 (MQ=255)
                                                                                    agTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAtc                                                                <  3:271186/149‑1 (MQ=255)
                                                                                      tttGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTAtctc                                                              >  2:649474/1‑149 (MQ=255)
                                                                                       ttGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTcc                                                             <  4:678827/149‑1 (MQ=255)
                                                                                          cACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACg                                                          >  3:869401/1‑149 (MQ=255)
                                                                                           aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                        <  3:249669/134‑1 (MQ=255)
                                                                                           aCCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTc                                                                        >  4:249669/1‑134 (MQ=255)
                                                                                                cgaAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTttatt                                                    <  2:27374/149‑1 (MQ=255)
                                                                                                   aGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                  >  2:375311/1‑148 (MQ=255)
                                                                                                   aGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                  <  1:375311/148‑1 (MQ=255)
                                                                                                    gCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGa                                                <  2:549894/149‑1 (MQ=255)
                                                                                                     cAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGAc                                               <  1:539606/149‑1 (MQ=255)
                                                                                                        gAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTc                                            >  3:456098/1‑149 (MQ=255)
                                                                                                         aaaTTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCa                                           <  2:265408/149‑1 (MQ=255)
                                                                                                          aaTTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCAc                                          <  4:124570/149‑1 (MQ=255)
                                                                                                               tCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGCCTGGCCCGTGGTAATTACGGTGCCCCGGAGGCGGAAACGCTACGCAGCGAGCATCAATATTGTCGGTTTGTCGCCGCGCATATTAAGACGTCACCAtta                                     >  2:203803/1‑148 (MQ=255)
                                                                                                                 tGTTCACGCTCGTTTGGAGAACGAATC‑ACGGATTATCTGTCGATGCGGCTGGCCATTTGTAGTTACGCTGCCGGGGCGGCGGAAAAACTTCGCAGCGCGCATCGATATTTTCGGTTTATCTCCACGTTTATTAATACGTCACCAGTTgc                                   >  2:590814/1‑149 (MQ=255)
                                                                                                                      ccGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTATGACGTCACCATTTGCGCTCa                              >  3:389326/1‑149 (MQ=255)
                                                                                                                       cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTGTTAAGACGTCACCATTTGCGCTCaa                             <  2:119293/149‑1 (MQ=255)
                                                                                                                       cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCaa                             >  4:139327/1‑149 (MQ=255)
                                                                                                                       cGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCaa                             <  2:478100/149‑1 (MQ=255)
                                                                                                                        gCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAAt                            >  1:181052/1‑149 (MQ=255)
                                                                                                                           cGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGaa                         >  3:682439/1‑149 (MQ=255)
                                                                                                                            gTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                        <  3:382532/149‑1 (MQ=255)
                                                                                                                            gTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                        <  4:110083/149‑1 (MQ=255)
                                                                                                                            gTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                        <  4:110274/149‑1 (MQ=255)
                                                                                                                            gTTTGGTGAACGCATC‑ACGGATTATCCGGCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAAc                        <  2:712676/149‑1 (MQ=255)
                                                                                                                                  tGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTatg                                                                                                                     >  3:663877/1‑48 (MQ=255)
                                                                                                                                  tGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACg                                                                                                                     <  4:663877/50‑1 (MQ=255)
                                                                                                                                        cATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCa            >  3:549768/1‑149 (MQ=255)
                                                                                                                                         atcgaCGGATAATCCGT‑GATGCGTCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCaa           >  1:185283/5‑149 (MQ=255)
                                                                                                                                         aTC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCaa           <  1:649474/149‑1 (MQ=255)
                                                                                                                                          tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt          >  3:173603/1‑149 (MQ=255)
                                                                                                                                          tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt          <  4:762040/149‑1 (MQ=255)
                                                                                                                                          tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt          >  4:830797/1‑149 (MQ=255)
                                                                                                                                          tC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAAt          >  2:438815/1‑149 (MQ=255)
                                                                                                                                           c‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATa         <  4:682439/149‑1 (MQ=255)
                                                                                                                                                 aTTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                  >  1:855609/1‑103 (MQ=255)
                                                                                                                                                 aTTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTaa                                                  <  2:855609/103‑1 (MQ=255)
                                                                                                                                                 aTTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGcgcg    >  4:388605/1‑149 (MQ=255)
                                                                                                                                                   tATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTc  >  4:151945/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATC‑ACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCCCGGGCGGCGGAAAAACTTCGCAGCGAGCATCAATATTGTCGGTTTATCTCCACGTTTATTAAGACGTCACCATTTGCGCTCAATGAACCTTATTACGGCAATAGCGCGTC  >  NC_000913/1231817‑1232111

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