Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,783 T→C 5.5% V64A (GTC→GCC)  ygbL → putative class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7830TC 5.5% 262.0 / 5.1 92V64A (GTC→GCC) ygbLputative class II aldolase
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base C (3/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.91e-01

ACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCT  >  NC_000913/2863638‑2863930
                                                                                                                                                 |                                                                                                                                                   
aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                                  >  2:423189/1‑149 (MQ=255)
aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                                  >  1:459518/1‑149 (MQ=255)
aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                                  >  1:803793/1‑149 (MQ=255)
aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACCGGAATTTACTTGCGCCACCGACAGGTTCATCCCGGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                                  >  2:787288/1‑149 (MQ=255)
   cgGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCcgctgc                                                                                                                                               <  2:754771/149‑1 (MQ=255)
     gATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAAGCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCgg                                                                                                                                             <  2:74026/149‑1 (MQ=255)
       ttGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGTCTGGCAATCTGTCGATGCTTTTACTTGACGGTAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGATAGCCGCTGCGGAt                                                                                                                                           <  2:190640/149‑1 (MQ=255)
       ttGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGATCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGAt                                                                                                                                           <  2:110851/149‑1 (MQ=255)
       ttGCCAGTTCATTCTTGCAGCGCGGCTATGCAACCGGTGCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGAt                                                                                                                                           <  1:1000407/149‑1 (MQ=255)
         gCCAGTTCAATCTGTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCAAGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATgg                                                                                                                                         <  1:1096445/149‑1 (MQ=255)
          ccAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTTCTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                        >  2:1084861/1‑149 (MQ=255)
          ccAGTTCACTCTTTCAGCGCGGCTATGCAACCGGAGCGGCTGGCAAGCTGTCGCGGCTTGTACCTGCCGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                                                                        <  1:785380/149‑1 (MQ=255)
            aGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGa                                                                                                                                      >  2:138598/1‑149 (MQ=255)
            aGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGa                                                                                                                                      >  1:622214/1‑149 (MQ=255)
             gTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGaa                                                                                                                                     >  1:413283/1‑149 (MQ=255)
              ttCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAAt                                                                                                                                    <  2:411902/149‑1 (MQ=255)
               tCATTCTTGCAGCGCGGCTAAGCAACCGGTTCGGCTGGCAATCTGTCGATGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATg                                                                                                                                   <  1:327081/149‑1 (MQ=255)
                 aTTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGt                                                                                                                                 >  2:1131073/1‑149 (MQ=255)
                  ttCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGtt                                                                                                                                >  2:776035/1‑149 (MQ=255)
                    cTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGGTGCTTTTACCTGACGGGAATTTACTGGCGACACAGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGAGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTaa                                                                                                                              >  2:1165240/1‑149 (MQ=255)
                    cTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTaa                                                                                                                              >  2:8921/1‑149 (MQ=255)
                        cAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtgg                                                                                                                          <  1:1109783/149‑1 (MQ=255)
                          gcgcGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACAACGACAGGTTCATGACTGGGCAATCTCGATCCGCAGCGGCTTTCAAAAGTCGCCGCGGATGGCGAATGGTTAAgtggtg                                                                                                                        >  2:468384/1‑149 (MQ=255)
                              ggCAATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCGTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACaa                                                                                                                    <  2:118281/149‑1 (MQ=255)
                                     cAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGCCACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCg                                                                                                             >  1:884800/1‑149 (MQ=255)
                                       aCCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGaa                                                                                                           <  2:402277/149‑1 (MQ=255)
                                         cGGTTCGGGTGGCAATCTGTCGCTGCTTTTACCTGTCGGGAATTTACTGGCGACACCGAGAGGTTCATGCCTGGGCAATCTCTATCCGCAGCGGCTTTCCAAAGTCTCCGCGGATGGCGTATGGTTAAGTGGTGACAAAGGCTCGAAag                                                                                                         <  1:132244/149‑1 (MQ=255)
                                         cGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACCCCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAag                                                                                                         <  1:1064625/149‑1 (MQ=255)
                                         cGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAag                                                                                                         >  1:135584/1‑149 (MQ=255)
                                          ggTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAaga                                                                                                        <  1:910543/149‑1 (MQ=255)
                                           gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACCCCGACAGGTTCAAGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAACTCGCCGCGGATGGAGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                       >  2:1047329/1‑149 (MQ=255)
                                           gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCTAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                       <  1:16146/149‑1 (MQ=255)
                                           gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                       >  1:404967/1‑149 (MQ=255)
                                           gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                       <  1:632691/149‑1 (MQ=255)
                                           gTTCGGCTGGCAATCTGGAGCTGCTTTTACCGGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCTCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                       <  1:645391/149‑1 (MQ=255)
                                             tCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGt                                                                                                     <  1:7176/149‑1 (MQ=255)
                                             tCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGt                                                                                                     <  1:762096/149‑1 (MQ=255)
                                                 cTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGctc                                                                                                 >  1:562711/1‑149 (MQ=255)
                                                   ggCAAGCTGGCGCGGCTTTTACCTGACGGGAATGTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGGGGGGACAAACCCTCGAAAGAGGTGCTCtt                                                                                               <  1:1144224/149‑1 (MQ=255)
                                                           gTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTgg                                                                                       >  2:1130405/1‑149 (MQ=255)
                                                            tCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                                                                      >  1:1160887/1‑149 (MQ=255)
                                                            tCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                                                                      <  1:95737/149‑1 (MQ=255)
                                                            tCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                                                                      >  1:285842/1‑149 (MQ=255)
                                                             cGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgcg                                                                                     <  1:581234/149‑1 (MQ=255)
                                                             cGCTGCTTTTACCTGACGGCAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCTAATGGTTAAGTGGTGACTAACGCTCGTAAGAGGTGCTCTTTCATCTGgcg                                                                                     <  2:558523/149‑1 (MQ=255)
                                                               ctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCt                                                                                   >  1:979803/1‑149 (MQ=255)
                                                                  ctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTa                                                                                >  2:725236/1‑149 (MQ=255)
                                                                       aCCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAACCCCTCGAAAGAGGTGCTCTTTAATCTGGCGCTGTATCGca                                                                           >  1:227103/1‑149 (MQ=255)
                                                                         cTGACGGGAATTTACTGTCGACACCGCAAGGTTAATTCCTGGTAAATATATATCCGCAGCGGCTTTACAAAGTCGCCGCGGATCGCGAATGGTTAAGTGGTGACAAACGCTCGAAAGAGGTGCACTTTCAAATGGCGCTGTATCGcatc                                                                         >  2:623570/1‑147 (MQ=255)
                                                                               ggAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCACTGTATCGCAACAATCcg                                                                   >  2:419264/1‑149 (MQ=255)
                                                                                gAATTTACTGGCGACACCGACAAGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCcgc                                                                  >  2:678986/1‑149 (MQ=255)
                                                                                 aaTTTACTGGCGACACCGAAAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCcgcg                                                                 <  2:1045542/149‑1 (MQ=255)
                                                                                  aTTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCgcgc                                                                <  1:26748/149‑1 (MQ=255)
                                                                                          gCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc                                                        >  2:94638/1‑149 (MQ=255)
                                                                                          gCGACACCGACAGGTTCATGCCTGGGCAATCTCGAGCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc                                                        <  2:80068/149‑1 (MQ=255)
                                                                                             acacCGTCAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggt                                                     >  2:413457/1‑149 (MQ=255)
                                                                                              cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAACCCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtg                                                    <  1:1165240/149‑1 (MQ=255)
                                                                                              cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtg                                                    <  1:898702/149‑1 (MQ=255)
                                                                                                ccGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTTTAAAGCggtggt                                                  <  1:657743/149‑1 (MQ=255)
                                                                                                ccGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCggtggt                                                  >  2:572489/1‑149 (MQ=255)
                                                                                                  gACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGc                                                <  2:803793/149‑1 (MQ=255)
                                                                                                      gggtCATGCCTGGGCAATCTCGATCCGCCGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAttt                                            <  1:1066304/146‑1 (MQ=255)
                                                                                                      ggTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCAttt                                            <  1:160377/149‑1 (MQ=255)
                                                                                                      ggTTCATGCCTGGGCAATCTCGATCCGCAGCGGCGTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGGGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCGGTAAAGCGGTGGTGCAttt                                            <  1:222731/149‑1 (MQ=255)
                                                                                                       gTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGATGTAAAGCGGTGGTGCATTTg                                           >  1:730813/1‑149 (MQ=255)
                                                                                                         tCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgca                                         <  2:562711/149‑1 (MQ=255)
                                                                                                         tCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgca                                         <  2:404967/149‑1 (MQ=255)
                                                                                                           aTGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTTGCGCTGTATCGCCACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcaaa                                       >  1:333514/1‑147 (MQ=255)
                                                                                                            tGCATGGGCAATCTCGATCCGCAGCGGCTTACCAAAGTCGCCGCGGATGGAGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCCCAACAAACCGCGCTGTAAAGCGGTGGTGCATTTgcacag                                      <  1:806155/149‑1 (MQ=255)
                                                                                                              ccTGGGCAATCTCGTTCCGCAGCGGCTTTCCAAAGTCGCCGCGTATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagca                                    <  2:903493/149‑1 (MQ=255)
                                                                                                              ccTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagca                                    <  1:879855/149‑1 (MQ=255)
                                                                                                              ccTGGGCAATCTCGATCCGCAGCGGCCTGCCACAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACACACCCTCGAAAGCGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagca                                    <  2:979519/149‑1 (MQ=255)
                                                                                                               cTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgcacagcac                                   >  1:1005880/1‑149 (MQ=255)
                                                                                                                  ggCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGGTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATg                                <  1:776035/149‑1 (MQ=255)
                                                                                                                     aaTCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTc                             >  1:29216/1‑149 (MQ=255)
                                                                                                                     aaTCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTAATCTGGCGCTGTATCGCAACAATCCCCGATGTTAAGCGGTGGTCCATTTGCACAGCACATGGTc                             >  1:990327/1‑149 (MQ=255)
                                                                                                                       tctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGAGGGCGAATGGTTAAATGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTTTCGCAACAATCCGCGCTGTAAAGCGGTGGTGCTTTTGCACAGCACATGGTCgt                           >  2:100576/1‑148 (MQ=255)
                                                                                                                        ctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGAc                          >  1:276367/1‑149 (MQ=255)
                                                                                                                         gccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAACCCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACg                         >  2:715531/4‑149 (MQ=38)
                                                                                                                           gATCCGGAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                       <  2:285842/149‑1 (MQ=255)
                                                                                                                           gATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGgc                       >  2:671875/1‑149 (MQ=255)
                                                                                                                           caTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGCCGgc                       >  1:1069913/2‑149 (MQ=255)
                                                                                                                                gcggCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCttt                  >  1:166310/4‑149 (MQ=255)
                                                                                                                                   gCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcct               <  1:1130405/149‑1 (MQ=255)
                                                                                                                                   gCGGCTTTACAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAACCCCTCGAAAGCGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAACGCGGTGGTGCATTTGCACAGCACATGGTCGAAGGCGCTTTcct               <  2:779086/149‑1 (MQ=255)
                                                                                                                                     ggCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACATTCCGCGCTGTAAAGTGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcctgc             <  2:413283/149‑1 (MQ=255)
                                                                                                                                     ggCTTTCCAAAGTCGCAGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTcctgc             <  2:569760/149‑1 (MQ=255)
                                                                                                                                           ccAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCaa       >  2:248158/1‑149 (MQ=255)
                                                                                                                                            cAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTATATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAg      >  1:1017457/1‑149 (MQ=255)
                                                                                                                                              aaGGCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCTCAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAggg    <  1:423189/149‑1 (MQ=255)
                                                                                                                                                gTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCt  <  1:631028/149‑1 (MQ=255)
                                                                                                                                                gTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCt  >  2:302924/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAGGGCT  >  NC_000913/2863638‑2863930

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: