Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,145 T→C 5.5% L35P (CTC→CCC)  cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,1450TC 5.5% 425.5 / 9.5 128L35P (CTC→CCC) cydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (65/56);  new base C (5/2);  total (70/58)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01

Only 100 of 132 total aligned reads displayed.

TCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCC  >  NC_000913/772997‑773293
                                                                                                                                                    |                                                                                                                                                    
 cTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTc                                                                                                                                                     <  3:166857/149‑1 (MQ=255)
  ttCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCa                                                                                                                                                    >  1:581189/1‑149 (MQ=255)
   tCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCAc                                                                                                                                                   >  2:437823/1‑149 (MQ=255)
     cACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCAccc                                                                                                                                                 <  4:271142/149‑1 (MQ=255)
     cACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAccc                                                                                                                                                 <  3:417836/149‑1 (MQ=255)
      aCGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCg                                                                                                                                                >  1:296261/1‑149 (MQ=255)
      aCGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAcccc                                                                                                                                                <  4:424649/149‑2 (MQ=255)
        gACTACTCAGCCAGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTCGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGtt                                                                                                                                              >  4:580090/1‑149 (MQ=255)
          ctactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTc                                                                                                                                            <  1:119375/149‑1 (MQ=255)
            actCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCt                                                                                                                                          <  4:727436/149‑1 (MQ=255)
            actCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCt                                                                                                                                          <  2:80924/149‑1 (MQ=255)
            actCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCt                                                                                                                                          <  1:597309/149‑1 (MQ=255)
             ctCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATAATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTc                                                                                                                                         <  4:118983/149‑1 (MQ=255)
              tCAGCCGGCTGGCTAAGAAAGGAGCCGTCAAATGATCGATTATGAAGGATTGCGTTTTATCTAGCGGTTGCTGGTTGGCGTTCTGCTGATTGGTGTTGCAGTCACTGCCGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCg                                                                                                                                        >  4:796006/1‑149 (MQ=255)
              tCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCg                                                                                                                                        >  4:903821/1‑149 (MQ=255)
              tCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCg                                                                                                                                        >  1:719018/1‑149 (MQ=255)
               cAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCgg                                                                                                                                       <  3:731917/149‑1 (MQ=255)
                aGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                      <  3:653186/149‑1 (MQ=255)
                 gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                                     <  3:75612/149‑1 (MQ=255)
                 gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                                     <  3:264862/149‑1 (MQ=255)
                 gCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTc                                                                                                                                     <  4:679608/149‑1 (MQ=255)
                   cGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGt                                                                                                                                   <  4:497484/149‑1 (MQ=255)
                   cGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGt                                                                                                                                   <  4:657144/149‑1 (MQ=255)
                     gCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGtt                                                                                                                                              <  1:774143/136‑1 (MQ=255)
                     gCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGtt                                                                                                                                              >  2:774143/1‑136 (MQ=255)
                     gCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTaa                                                                                                                                 <  2:284463/149‑1 (MQ=255)
                     gCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAAc                                                                                                                                <  3:97874/149‑1 (MQ=255)
                           tAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTc                                                                                                                                                     <  2:839461/123‑1 (MQ=255)
                           tAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTc                                                                                                                                                     >  1:839461/1‑123 (MQ=255)
                           tAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGacac                                                                                                                           >  4:895814/1‑149 (MQ=255)
                              gACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGa                                                                                                                        <  2:789179/149‑1 (MQ=255)
                                 aGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCg                                                                                                                     >  1:364777/1‑149 (MQ=255)
                                     gtcgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGa                                                                                                                 >  3:110817/1‑149 (MQ=255)
                                     gtcgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGa                                                                                                                 <  3:1075/149‑1 (MQ=255)
                                        gtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                              >  3:641402/1‑149 (MQ=255)
                                        gtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGATGCAGGTTGGCGTTCTGCTGCTTGGTTTTGCAGTCACTGACGGTTTCGACCTGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                              >  2:814458/1‑149 (MQ=255)
                                         tcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTa                                                                                                             >  4:8269/1‑149 (MQ=255)
                                          cAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTAt                                                                                                            <  1:183245/149‑1 (MQ=255)
                                          cAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTAt                                                                                                            <  2:484472/149‑1 (MQ=255)
                                               gATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTa                                                                                                       >  3:594021/1‑149 (MQ=255)
                                                aTCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTaa                                                                                                      >  1:763166/1‑149 (MQ=255)
                                                  cGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACt                                                                                                    >  4:431694/1‑149 (MQ=255)
                                                   gATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTaa                                                                                                                                 >  2:454384/1‑119 (MQ=255)
                                                   gATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTaa                                                                                                                                 <  1:454384/119‑1 (MQ=255)
                                                    aTTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                              <  4:89087/137‑1 (MQ=255)
                                                       aTGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCAtt                                                                                               <  3:676363/149‑1 (MQ=255)
                                                        tGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGAt                                                                                                         >  3:422722/1‑138 (MQ=255)
                                                         gAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGc                                                                                             >  4:43740/1‑149 (MQ=255)
                                                         gAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGc                                                                                             >  1:157200/1‑149 (MQ=255)
                                                           aGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGcac                                                                                           <  4:731233/149‑1 (MQ=255)
                                                               ttGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACcaca                                                                                       >  1:700849/1‑149 (MQ=255)
                                                                  cGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGCTTAACTCCATTGCACCACACTg                                                                                    <  1:567685/149‑1 (MQ=255)
                                                                   gTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTgg                                                                                   >  1:651440/1‑149 (MQ=255)
                                                                     tttATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGa                                                                                 <  4:490940/149‑1 (MQ=255)
                                                                               ggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAg                                                                       <  4:766232/149‑1 (MQ=255)
                                                                                  tgctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                                    >  2:622438/1‑149 (MQ=255)
                                                                                    ctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTg                                                                  <  2:296261/149‑1 (MQ=255)
                                                                                     tgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTgg                                                                 <  1:837334/149‑1 (MQ=255)
                                                                                       gTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCt                                                               >  2:70456/1‑149 (MQ=255)
                                                                                       gTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCt                                                               <  1:581380/149‑1 (MQ=255)
                                                                                         tGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGa                                                             <  3:895814/149‑1 (MQ=255)
                                                                                          ggCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGAt                                                            <  4:352180/149‑1 (MQ=255)
                                                                                          ggCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGAt                                                            >  3:591984/1‑149 (MQ=255)
                                                                                             gTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCAc                                                         >  2:38428/1‑149 (MQ=255)
                                                                                             gTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTATCAcc                                                        >  2:627721/1‑149 (MQ=255)
                                                                                              ttCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCAcc                                                        >  4:204199/1‑149 (MQ=255)
                                                                                              ttCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCAcc                                                        >  3:92012/1‑149 (MQ=255)
                                                                                                 tgctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcgcg                                                     >  1:23973/1‑149 (MQ=255)
                                                                                                      aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                                >  2:421738/1‑149 (MQ=255)
                                                                                                      aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                                >  3:781131/1‑149 (MQ=255)
                                                                                                      aTTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                                >  1:18985/1‑149 (MQ=255)
                                                                                                       ttGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGggcggc                                               >  3:129350/1‑149 (MQ=255)
                                                                                                       ttGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGggcggc                                               >  4:190588/1‑149 (MQ=255)
                                                                                                         ggTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTTGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGgcgc                                             >  1:374896/1‑149 (MQ=255)
                                                                                                          gTTTTGCAGTCACTGACGGCTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCa                                            >  1:703821/1‑149 (MQ=255)
                                                                                                           ttttGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGTATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCAc                                           >  4:604062/1‑149 (MQ=255)
                                                                                                                  gTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtgc                                    <  1:814458/149‑1 (MQ=255)
                                                                                                                   tCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtgc                                    <  1:388385/148‑1 (MQ=255)
                                                                                                                     aCTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtgctgc                                 >  4:376093/1‑149 (MQ=255)
                                                                                                                         aCGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTgg                             >  2:355136/1‑149 (MQ=255)
                                                                                                                          cGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGc                            <  3:376093/149‑1 (MQ=255)
                                                                                                                            gTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCg                          >  1:470826/1‑149 (MQ=255)
                                                                                                                              ttCGACATGGGGGTGGGCATGCTCACCCGTTTCCTTGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGAt                        <  2:374896/149‑1 (MQ=255)
                                                                                                                               tCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATg                       <  3:263326/149‑1 (MQ=255)
                                                                                                                                  aCATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTc                    <  1:421738/149‑1 (MQ=255)
                                                                                                                                   cATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGAc                                                                                <  2:323264/88‑1 (MQ=255)
                                                                                                                                   cATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGAc                                                                                >  1:323264/1‑88 (MQ=255)
                                                                                                                                   cATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCt                   <  1:610870/149‑1 (MQ=255)
                                                                                                                                    aTGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTa                  >  2:149776/1‑149 (MQ=255)
                                                                                                                                     tGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTAt                 <  2:157200/149‑1 (MQ=255)
                                                                                                                                      acgggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATg                >  1:851539/3‑149 (MQ=255)
                                                                                                                                      acgggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATg                >  2:107365/3‑149 (MQ=255)
                                                                                                                                        gggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGcc              >  4:625106/1‑149 (MQ=255)
                                                                                                                                        gggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGcc              >  4:286234/1‑149 (MQ=255)
                                                                                                                                        gggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGcc              >  3:600009/1‑149 (MQ=255)
                                                                                                                                              gCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCGATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCg        <  4:591984/149‑1 (MQ=255)
                                                                                                                                                  gCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTtct    >  4:591235/1‑149 (MQ=255)
                                                                                                                                                    tCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCt                               <  1:707432/120‑1 (MQ=255)
                                                                                                                                                    tCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCt                               >  2:707432/1‑120 (MQ=255)
                                                                                                                                                    tCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTcc  >  4:81366/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCC  >  NC_000913/772997‑773293

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