Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,778 C→G 5.6% S62S (TCC→TCG) ‡ ygbL → putative class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7780CG 5.6% 246.9 / 6.4 90S62S (TCC→TCG) ‡ygbLputative class II aldolase
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base G (3/2);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

GAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAG  >  NC_000913/2863632‑2863926
                                                                                                                                                  |                                                                                                                                                    
gagaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaa                                                                                                                                                    >  2:910543/1‑149 (MQ=255)
 tgaggACGCGGATTGCCAGTTCATTCTTTCTGCGCGGCTATGCAACCGTTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaaa                                                                                                                                                   <  2:253695/145‑1 (MQ=255)
 agaTGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaaa                                                                                                                                                   >  2:631028/1‑149 (MQ=255)
      aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                              >  2:423189/1‑149 (MQ=255)
      aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                              >  1:803793/1‑149 (MQ=255)
      aCGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgc                                                                                                                                              >  1:459518/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                                                                                                                                  >  1:622214/1‑149 (MQ=255)
                  aGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGa                                                                                                                                  >  2:138598/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:632691/149‑1 (MQ=255)
                                                 gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAagag                                                                                                   >  1:404967/1‑149 (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:285842/1‑149 (MQ=255)
                                                                  tCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                                                                  <  1:95737/149‑1 (MQ=255)
                                                                  tCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                                                                  >  1:1160887/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)
                                                                                                                                                  |                                                                                                                                                    
GAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGCCTGCAAG  >  NC_000913/2863632‑2863926

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: ‑