Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,885,300 T→G 5.8% R75R (CGT→CGG) ‡ yidC → membrane protein insertase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,885,3000TG 5.8% 414.3 / 12.9 127R75R (CGT→CGG) ‡yidCmembrane protein insertase
Reads supporting (aligned to +/- strand):  ref base T (55/64);  new base G (5/3);  total (60/67)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 135 total aligned reads displayed.

ATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGC  >  NC_000913/3885152‑3885446
                                                                                                                                                    |                                                                                                                                                  
aTAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgt                                                                                                                                                    <  4:881984/149‑1 (MQ=255)
   aaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggt                                                                                                                                                 <  4:703016/149‑1 (MQ=255)
    aaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtg                                                                                                                                                >  3:860198/1‑149 (MQ=255)
      aCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                              <  1:663049/149‑1 (MQ=255)
      aCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                              <  1:552237/149‑1 (MQ=255)
       cccGCAACCTCAGGCTCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggtg                                                                                                                                             <  3:994601/149‑1 (MQ=255)
        ccGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGa                                                                                                                                            >  3:759814/1‑149 (MQ=255)
         cGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGAt                                                                                                                                           >  3:193241/1‑149 (MQ=255)
          gCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCACCACCACGGGTGTTGATg                                                                                                                                          <  2:428764/149‑1 (MQ=255)
          gCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATg                                                                                                                                          >  3:203614/1‑149 (MQ=255)
           cAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGt                                                                                                                                         <  2:832652/149‑1 (MQ=255)
           cAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGt                                                                                                                                         <  2:507991/149‑1 (MQ=255)
            aaCCTCAGGCCGAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAGCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                        <  2:461889/149‑1 (MQ=255)
              ccTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGa                                                                                                                                      <  3:152297/149‑1 (MQ=255)
                  aGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCaa                                                                                                                                  >  1:210274/1‑149 (MQ=255)
                      ccAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGctc                                                                                                                              >  1:434948/1‑149 (MQ=255)
                       cAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGGGGTGATGTTGAGCAAGctct                                                                                                                             >  4:366508/1‑149 (MQ=255)
                        aaCAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAGACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctg                                                                                                                            >  3:51112/1‑149 (MQ=255)
                        aaCAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctg                                                                                                                            >  4:333955/1‑149 (MQ=255)
                         aCAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTctgc                                                                                                                           <  3:311060/149‑1 (MQ=255)
                            gACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCtgctgc                                                                                                                        <  4:25984/149‑1 (MQ=255)
                                 cGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCt                                                                                                                   >  3:12347/1‑149 (MQ=255)
                                  gCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCACCACCACGGGTGTTGATg                                                                                                                                          <  2:805736/125‑1 (MQ=255)
                                  gCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCACCACCACGGGTGTTGATg                                                                                                                                          >  1:805736/1‑125 (MQ=255)
                                    aGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAc                                                                                                                <  3:129396/149‑1 (MQ=255)
                                      aCAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAccc                                                                                                              <  1:828578/149‑1 (MQ=255)
                                      aCAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAccc                                                                                                              >  2:396796/1‑149 (MQ=255)
                                        aaCGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                        >  1:355806/1‑121 (MQ=255)
                                        aaCGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                            >  1:461945/1‑149 (MQ=255)
                                          cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                          <  2:970156/149‑1 (MQ=255)
                                          cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                          <  3:947332/149‑1 (MQ=255)
                                            accaccGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAaga                                                                                                        >  1:723949/1‑149 (MQ=255)
                                               accGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCt                                                                                                     >  3:597278/1‑149 (MQ=255)
                                               accGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCt                                                                                                     <  3:868855/149‑1 (MQ=255)
                                               accGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCt                                                                                                     <  2:543410/149‑1 (MQ=255)
                                                ccGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTg                                                                                                    <  4:431011/149‑1 (MQ=255)
                                                  gCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGaa                                                                                                  <  3:111009/149‑1 (MQ=255)
                                                  gCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGaa                                                                                                  <  4:446731/149‑1 (MQ=255)
                                                   cAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAAc                                                                                                 <  4:597278/149‑1 (MQ=255)
                                                       gggTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTa                                                                                             >  1:594344/1‑149 (MQ=255)
                                                                     aCCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAg                                                                               <  3:363719/149‑1 (MQ=255)
                                                                         gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTa                                                                                             <  2:56497/131‑1 (MQ=255)
                                                                         gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTa                                                                                             >  1:56497/1‑131 (MQ=255)
                                                                         gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCCGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGt                                                                           <  3:720460/149‑1 (MQ=255)
                                                                          ggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                          <  1:337602/149‑1 (MQ=255)
                                                                          ggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                          <  3:526166/149‑1 (MQ=255)
                                                                             gTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCCTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                          >  1:189906/1‑114 (MQ=255)
                                                                             gTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCCTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                          <  2:189906/114‑1 (MQ=255)
                                                                                ccGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAAc                                                                    <  4:759814/149‑1 (MQ=255)
                                                                                ccGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAAc                                                                    <  2:1334/149‑1 (MQ=255)
                                                                                  ggCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACtt                                                                  <  4:515511/149‑1 (MQ=255)
                                                                                     cAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAc                                                               <  2:456416/149‑1 (MQ=255)
                                                                                       gTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCg                                                             >  3:589837/1‑149 (MQ=255)
                                                                                           ccAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGt                                                         >  1:653588/1‑149 (MQ=255)
                                                                                            cAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCa                                                                >  3:867292/1‑141 (MQ=255)
                                                                                              tgggAAACGGATCTCGGTTAAGTCCGACGTGCTTGATCTGACCATCAACAGCCGTGGTTGTTATGTTGAGCAAGCTCTGCTGCCTGCTTAGCCAAAAGAGCTAACCTCTACCCAGCCGTTTCAGCTATTGGAAACTTCACCGCAGttta                                                      <  4:647480/148‑1 (MQ=255)
                                                                                                 gAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttattt                                                   <  2:461945/149‑1 (MQ=255)
                                                                                                 gAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttattt                                                   >  3:401344/1‑149 (MQ=255)
                                                                                                     cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCa                                               <  1:380586/149‑1 (MQ=255)
                                                                                                         tctcGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGca                                           >  3:515863/1‑149 (MQ=255)
                                                                                                         tctcGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGca                                           >  1:995992/1‑149 (MQ=255)
                                                                                                           tcGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTATCAGGCACag                                        >  4:413619/1‑149 (MQ=255)
                                                                                                             ggTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAACCTTCACCGCAGTTTATTTATCAGGCACaga                                       <  2:434948/149‑1 (MQ=255)
                                                                                                             ggTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACaga                                       <  1:251155/149‑1 (MQ=255)
                                                                                                                  aGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGt                                  >  3:934303/1‑149 (MQ=255)
                                                                                                                  aGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGt                                  <  1:323451/149‑1 (MQ=255)
                                                                                                                   gACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAgg                                             >  4:31408/1‑137 (MQ=255)
                                                                                                                   gACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAgg                                             <  3:31408/137‑1 (MQ=255)
                                                                                                                      cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                              >  3:346531/1‑149 (MQ=255)
                                                                                                                      cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                              <  4:737437/149‑1 (MQ=255)
                                                                                                                      cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                              <  4:554080/149‑1 (MQ=255)
                                                                                                                            gCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTc                        <  1:577720/149‑1 (MQ=255)
                                                                                                                            gCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTc                        <  3:390414/149‑1 (MQ=255)
                                                                                                                              ttGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGt                      >  2:680877/1‑149 (MQ=255)
                                                                                                                              ttGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGt                      >  3:113876/1‑149 (MQ=255)
                                                                                                                               tGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTATCAGGCACAGAGCGGTCTGACCGGTCGTGa                    <  3:413619/149‑1 (MQ=255)
                                                                                                                                gATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGa                    >  1:636444/1‑149 (MQ=255)
                                                                                                                                gATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGa                    <  4:277257/149‑1 (MQ=255)
                                                                                                                                 aTCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGAt                   >  1:130144/1‑149 (MQ=255)
                                                                                                                                   cTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATgg                 <  4:65593/149‑1 (MQ=255)
                                                                                                                                    tcaacATCACCACCACGGGTGTTGATGTTGAGGAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                >  4:898085/5‑149 (MQ=255)
                                                                                                                                    tGACCATCAACACCCGGGGTGGTGATGTTTAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc                >  4:812107/1‑149 (MQ=255)
                                                                                                                                     gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc               <  3:316372/149‑1 (MQ=255)
                                                                                                                                     gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc               >  3:663915/1‑149 (MQ=255)
                                                                                                                                     caacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc               >  4:581700/4‑149 (MQ=255)
                                                                                                                                      aCCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc              <  1:274250/149‑1 (MQ=255)
                                                                                                                                      aCCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc              <  3:287692/149‑1 (MQ=255)
                                                                                                                                      aCCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc              <  3:149030/149‑1 (MQ=255)
                                                                                                                                        cATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGttta                                                      >  2:752115/1‑107 (MQ=255)
                                                                                                                                        cATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGttta                                                      <  1:752115/107‑1 (MQ=255)
                                                                                                                                        cATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCgg            >  4:546953/1‑149 (MQ=255)
                                                                                                                                        cATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCgg            >  1:497337/1‑149 (MQ=255)
                                                                                                                                         aTCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                             <  2:940965/83‑1 (MQ=255)
                                                                                                                                         aTCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCt                                                                             >  1:940965/1‑83 (MQ=255)
                                                                                                                                          tCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATTAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGAt          >  2:704250/1‑149 (MQ=255)
                                                                                                                                             acacCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAc       >  2:312286/1‑149 (MQ=255)
                                                                                                                                               acCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGcc                                                                                      >  4:583338/1‑68 (MQ=255)
                                                                                                                                               acCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGcc                                                                                      <  3:583338/68‑1 (MQ=255)
                                                                                                                                               acCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAccc     >  2:896921/1‑149 (MQ=255)
                                                                                                                                                  cGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGc  >  1:97227/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGC  >  NC_000913/3885152‑3885446

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