Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,772 G→A 5.6% R60R (CGG→CGA)  ygbL → putative class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7720GA 5.6% 258.2 / 4.2 89R60R (CGG→CGA) ygbLputative class II aldolase
Reads supporting (aligned to +/- strand):  ref base G (41/43);  new base A (3/2);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGC  >  NC_000913/2863627‑2863919
                                                                                                                                                 |                                                                                                                                                   
gagaGGAGATGACGCGGATTGCCAGTTCATTCTTACAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCTCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGGTCATGCCTGGGCAATCTCGATCCGCAGCGGCtt                                                                                                                                                  >  1:787241/1‑149 (MQ=255)
     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:762096/149‑1 (MQ=255)
                                                        tCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGt                                                                                          <  1:7176/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:404967/149‑1 (MQ=255)
                                                                                                                    tCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTgca                              <  2:562711/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)
                                                                                                                                                 |                                                                                                                                                   
GAGAGGAGATGACGCGGATTGCCAGTTCATTCTTTCAGCGCGGCTATGCAACCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGCGGTGGTGCATTTGCACAGCACATGGTCGACGGCGCTTTCCTGC  >  NC_000913/2863627‑2863919

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