Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,885,293 A→C 5.5% N73T (AAC→ACC)  yidC → membrane protein insertase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,885,2930AC 5.5% 425.3 / 14.5 129N73T (AAC→ACC) yidCmembrane protein insertase
Reads supporting (aligned to +/- strand):  ref base A (54/67);  new base C (4/3);  total (58/71)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

Only 100 of 130 total aligned reads displayed.

AGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAC  >  NC_000913/3885149‑3885441
                                                                                                                                                |                                                                                                                                                    
aGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACAcc                                                                                                                                                  >  3:606954/1‑149 (MQ=255)
   aTAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgt                                                                                                                                               <  4:881984/149‑1 (MQ=255)
      aaaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggt                                                                                                                                            >  4:111009/1‑149 (MQ=255)
       aaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtg                                                                                                                                           >  3:860198/1‑149 (MQ=255)
       aaaCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTCGATCTGACCATCAACACCCgtggtg                                                                                                                                           >  3:809305/1‑149 (MQ=255)
         aCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                         <  1:663049/149‑1 (MQ=255)
         aCCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCgtggtggt                                                                                                                                         <  1:552237/149‑1 (MQ=255)
           ccGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGa                                                                                                                                       >  3:759814/1‑149 (MQ=255)
            cGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGAt                                                                                                                                      >  3:193241/1‑149 (MQ=255)
             gCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATg                                                                                                                                     >  3:203614/1‑149 (MQ=255)
              cAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGt                                                                                                                                    <  2:507991/149‑1 (MQ=255)
              cAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGt                                                                                                                                    <  2:832652/149‑1 (MQ=255)
               aaCCTCAGGCCGAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAGCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                   <  2:461889/149‑1 (MQ=255)
                 ccTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGa                                                                                                                                 <  3:152297/149‑1 (MQ=255)
                   tCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCACCACCACGGGTGTTGATggtcaga                                                                                                                               <  3:643063/149‑7 (MQ=255)
                     aGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCaa                                                                                                                             >  1:210274/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)
                                    cGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCt                                                                                                              >  1:223966/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                                                                                                         >  2:396796/1‑149 (MQ=255)
                                         aCAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTAccc                                                                                                         <  1:828578/149‑1 (MQ=255)
                                          cAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCg                                                                                                        >  2:976027/1‑149 (MQ=255)
                                           aaCGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                   <  2:355806/121‑1 (MQ=255)
                                           aaCGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGtt                                                                                                                                   >  1:355806/1‑121 (MQ=255)
                                           aaCGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                       >  1:461945/1‑149 (MQ=255)
                                           aaCGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGa                                                                                                       >  3:554080/1‑149 (MQ=255)
                                             cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                     <  2:970156/149‑1 (MQ=255)
                                             cGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                     <  3:947332/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)
                                                          gggTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcc                                                                                                                  <  3:50202/123‑1 (MQ=255)
                                                          gggTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGctgcc                                                                                                                  >  4:50202/1‑123 (MQ=255)
                                                          gggTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTa                                                                                        >  1:594344/1‑149 (MQ=255)
                                                                     ccgACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGATTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACTCGAAAGAGCTGAACTCTACCCAGCCGTTc                                                                             <  3:897676/149‑1 (MQ=255)
                                                                        aCCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAg                                                                          <  3:363719/149‑1 (MQ=255)
                                                                           aGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGACGTTCCAGCTg                                                                       >  4:374326/1‑149 (MQ=255)
                                                                            gggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTa                                                                                        <  2:56497/131‑1 (MQ=255)
                                                                             ggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                     <  3:374326/149‑1 (MQ=255)
                                                                             ggCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGtt                                                                     <  3:526166/149‑1 (MQ=255)
                                                                              gCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTg                                                                    <  4:726550/149‑1 (MQ=255)
                                                                                gTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCCTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGaaa                                                                                                     >  1:189906/1‑114 (MQ=255)
                                                                                   ccGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAAc                                                               <  2:1334/149‑1 (MQ=255)
                                                                                   ccGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAAc                                                               <  4:759814/149‑1 (MQ=255)
                                                                                        cAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAc                                                          <  2:456416/149‑1 (MQ=255)
                                                                                        cAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCAGCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCAc                                                          <  4:644740/149‑1 (MQ=255)
                                                                                          gTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCg                                                        >  3:589837/1‑149 (MQ=255)
                                                                                              ccAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGt                                                    >  1:653588/1‑149 (MQ=255)
                                                                                               cAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCa                                                           <  4:867292/141‑1 (MQ=255)
                                                                                               cAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCa                                                           >  3:867292/1‑141 (MQ=255)
                                                                                                 tgggAAACGGATCTCGGTTAAGTCCGACGTGCTTGATCTGACCATCAACAGCCGTGGTTGTTATGTTGAGCAAGCTCTGCTGCCTGCTTAGCCAAAAGAGCTAACCTCTACCCAGCCGTTTCAGCTATTGGAAACTTCACCGCAGttta                                                 <  4:647480/148‑1 (MQ=255)
                                                                                                  gggAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttat                                                <  1:74077/149‑1 (MQ=255)
                                                                                                    gAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttattt                                              >  3:401344/1‑149 (MQ=255)
                                                                                                    gAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttattt                                              <  2:461945/149‑1 (MQ=255)
                                                                                                       aCTGGTCTCGGTTAAGAGCGAGGTGCTTGATATGGCCGTCATTACCCGTGGTGATGATGTTGTGTAGGCTCTGCTGCGGTATCATTCGAATGAGCTGATCCATCTCCATCTGTTCCAGCTGTTGGAAACCTCACCGCAGTTTATTTATc                                           <  2:878172/149‑1 (MQ=255)
                                                                                                       aCTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATc                                           <  2:1022222/149‑1 (MQ=255)
                                                                                                        cTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCa                                          <  1:380586/149‑1 (MQ=255)
                                                                                                            tctcGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGca                                      >  1:995992/1‑149 (MQ=255)
                                                                                                              tcGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGcaca                                    >  1:331695/1‑149 (MQ=255)
                                                                                                                ggTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAACCTTCACCGCAGTTTATTTATCAGGCACaga                                  <  2:434948/149‑1 (MQ=255)
                                                                                                                     aGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGt                             <  1:323451/149‑1 (MQ=255)
                                                                                                                     aGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGt                             >  3:934303/1‑149 (MQ=255)
                                                                                                                      gACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAgg                                        >  4:31408/1‑137 (MQ=255)
                                                                                                                      gACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAgg                                        <  3:31408/137‑1 (MQ=255)
                                                                                                                         cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                         <  4:737437/149‑1 (MQ=255)
                                                                                                                         cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                         <  3:127804/149‑1 (MQ=255)
                                                                                                                         cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                         >  3:346531/1‑149 (MQ=255)
                                                                                                                         cGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGa                         <  4:554080/149‑1 (MQ=255)
                                                                                                                               gCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTc                   <  1:577720/149‑1 (MQ=255)
                                                                                                                                 ttGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGt                 >  2:680877/1‑149 (MQ=255)
                                                                                                                                  tGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTATCAGGCACAGAGCGGTCTGACCGGTCGTGa               <  3:413619/149‑1 (MQ=255)
                                                                                                                                   gATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGa               <  4:277257/149‑1 (MQ=255)
                                                                                                                                    aTCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttatt                                               >  2:592336/1‑116 (MQ=255)
                                                                                                                                    aTCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGtttatt                                               <  1:592336/116‑1 (MQ=255)
                                                                                                                                    aTCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGAt              >  1:130144/1‑149 (MQ=255)
                                                                                                                                    aTCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGAt              >  1:393305/1‑149 (MQ=255)
                                                                                                                                      cTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATgg            <  4:65593/149‑1 (MQ=255)
                                                                                                                                       tcaacATCACCACCACGGGTGTTGATGTTGAGGAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc           >  4:898085/5‑149 (MQ=255)
                                                                                                                                       tcaacATCACCACCACGGGTGTTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGc           >  4:562768/5‑149 (MQ=255)
                                                                                                                                        gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc          <  3:316372/149‑1 (MQ=255)
                                                                                                                                        gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc          <  1:499837/149‑1 (MQ=255)
                                                                                                                                        gACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGcc          >  3:663915/1‑149 (MQ=255)
                                                                                                                                         aCCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc         <  1:274250/149‑1 (MQ=255)
                                                                                                                                         aCCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGccc         <  3:287692/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)
                                                                                                                                           cATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCgg       >  1:659023/1‑149 (MQ=255)
                                                                                                                                             tCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATTAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGAt     >  2:704250/1‑149 (MQ=255)
                                                                                                                                                acacCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAc  >  2:312286/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGGATAAAAACCCGCAACCTCAGGCCCAACAGACCACGCAGACAACGACCACCGCAGCGGGTAGCGCCGCCGACCAGGGCGTACCGGCCAGTGGCCAGGGGAAACTGATCTCGGTTAAGACCGACGTGCTTGATCTGACCATCAACACCCGTGGTGGTGATGTTGAGCAAGCTCTGCTGCCTGCTTACCCGAAAGAGCTGAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAAC  >  NC_000913/3885149‑3885441

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: