Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,656,820 C→G 5.9% Q88H (CAG→CAC fdx ← [2Fe‑2S] ferredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,656,8200CG 5.9% 278.1 / 6.0 81Q88H (CAG→CACfdx[2Fe‑2S] ferredoxin
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base G (2/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGT  >  NC_000913/2656672‑2656965
                                                                                                                                                    |                                                                                                                                                 
aGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATAGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTAGACTACTAAATCTTCGTCGGTAACGCGCGcc                                                                                                                                                   >  7:175936/1‑149 (MQ=255)
 gATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCt                                                                                                                                                  <  5:96214/149‑1 (MQ=255)
 gATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCt                                                                                                                                                  <  3:105533/149‑1 (MQ=255)
  aTCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTg                                                                                                                                                 <  4:198012/149‑1 (MQ=255)
   tCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTgg                                                                                                                                                <  7:261580/149‑1 (MQ=255)
     gggATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCa                                                                                                                                              <  1:236314/149‑1 (MQ=255)
     gggATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCa                                                                                                                                              <  7:59313/149‑1 (MQ=255)
      ggATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAg                                                                                                                                             <  7:376901/149‑1 (MQ=255)
         tACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCtt                                                                                                                                          <  3:60224/149‑1 (MQ=255)
            gCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTaaa                                                                                                                                       <  8:233790/149‑1 (MQ=255)
               tCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACgg                                                                                                                                    >  2:34201/1‑149 (MQ=255)
                  tACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCtt                                                                                                                                 <  6:150701/149‑1 (MQ=255)
                   aCAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCttt                                                                                                                                >  7:233559/1‑149 (MQ=255)
                    cAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTc                                                                                                                               <  8:205224/149‑1 (MQ=255)
                     aGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTc                                                                                                                               <  7:315256/148‑1 (MQ=255)
                     aGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTc                                                                                                                               >  8:315256/1‑148 (MQ=255)
                             gCCAATTTCGCGGCTATCGGTCGACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAg                                                                                                                      >  8:73976/1‑149 (MQ=255)
                             gCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCa                                                                                                                       <  8:39444/148‑1 (MQ=255)
                              ccAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGt                                                                                                                     >  3:301752/1‑149 (MQ=255)
                              ccAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGt                                                                                                                     >  3:296808/1‑149 (MQ=255)
                               cAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTc                                                                                                                    >  3:90859/1‑149 (MQ=255)
                                 aTTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         <  6:202713/142‑1 (MQ=255)
                                 aTTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         >  5:202713/1‑142 (MQ=255)
                                  tttCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         <  4:19062/141‑1 (MQ=255)
                                  tttCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTc                                                                                                                         >  3:19062/1‑141 (MQ=255)
                                     cgcgGCTATAGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGATTTGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAg                                                                                                               <  1:64981/148‑1 (MQ=255)
                                      gcgGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAgg                                                                                                              <  4:90859/148‑1 (MQ=255)
                                         gCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGGTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCttt                                                                                                          <  7:353917/149‑1 (MQ=255)
                                          cTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTg                                                                                                         <  7:126949/149‑1 (MQ=255)
                                             tCGGTCCACTTAAGTCCCATACTAACCACTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTc                                                                                                       >  7:78347/1‑148 (MQ=255)
                                                gTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAg                                                                                                    >  8:231285/1‑148 (MQ=255)
                                                gTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAg                                                                                                    >  1:202505/1‑148 (MQ=255)
                                                gTCCACTTAAGACCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGc                                                                                                   <  5:317068/149‑1 (MQ=255)
                                                   cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                <  5:201612/149‑1 (MQ=255)
                                                   cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                <  8:78347/149‑1 (MQ=255)
                                                   cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                >  1:258114/1‑149 (MQ=255)
                                                   cACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATg                                                                                                >  3:131955/1‑149 (MQ=255)
                                                     cTTAAGTCCAATACTAACCTCTGTTAATGCTCAATCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATgtc                                                                                              <  7:341257/149‑1 (MQ=255)
                                                      ttAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATgtcg                                                                                             >  4:168309/1‑149 (MQ=255)
                                                         aGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATgtcgtc                                                                                           <  1:213403/148‑1 (MQ=255)
                                                         aGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCt                                                                                          <  8:233559/149‑1 (MQ=255)
                                                          gTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCtt                                                                                         >  1:319776/1‑149 (MQ=255)
                                                          gTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCtt                                                                                         >  1:185973/1‑149 (MQ=255)
                                                             ccATAGTAACGTCTGTTAATGCTCACGCGTATGGTTGATAGTGTAACGCGGGATTTCGACTACTATATCTTCGTAGGTAACGCGCGCCTGGCAGCTTTAACGTCTTTCCGGCTCCAGTCCCGTGGCTTTGTCCAGCATGTGGTCTTcc                                                                                       <  1:267572/148‑1 (MQ=255)
                                                              cATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTg                                                                                     >  1:231112/1‑149 (MQ=255)
                                                                tACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGct                                                                                   <  2:319776/149‑1 (MQ=255)
                                                                tACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGc                                                                                    >  6:114001/1‑148 (MQ=255)
                                                                tACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGGCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGct                                                                                   >  5:239609/1‑149 (MQ=255)
                                                                     aCCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTg                                                                                <  7:329977/147‑1 (MQ=255)
                                                                     aCCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAg                                                                              <  7:73976/149‑1 (MQ=255)
                                                                        tctGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCtt                                                                           >  8:17077/1‑149 (MQ=255)
                                                                           gTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGa                                                                               >  3:753/1‑142 (MQ=255)
                                                                           gTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGa                                                                               <  4:753/142‑1 (MQ=255)
                                                                           gTTAATGCTCACGCGCATGGTTGATAGCGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGGAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTcc                                                                        >  8:183245/1‑149 (MQ=255)
                                                                              aaTGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCgg                                                                      >  8:50952/1‑148 (MQ=255)
                                                                               aTGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCCGAGCTTTCCGGc                                                                     <  2:258114/148‑1 (MQ=255)
                                                                                tGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAg                                                                   <  4:102969/149‑1 (MQ=255)
                                                                                  cTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGt                                                                  <  1:210620/148‑1 (MQ=255)
                                                                                          cATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCaaaa                                                         <  7:231285/149‑1 (MQ=255)
                                                                                           aTGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAAc                                                        <  5:272953/149‑1 (MQ=255)
                                                                                              gTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCt                                                      >  4:114409/1‑148 (MQ=255)
                                                                                                   tAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCacg                                                 >  3:137922/1‑148 (MQ=255)
                                                                                                        tAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGAt                                           <  1:34201/149‑1 (MQ=255)
                                                                                                           cgcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATg                                          >  8:262929/1‑147 (MQ=255)
                                                                                                           cgcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCa                                        <  2:202505/149‑1 (MQ=255)
                                                                                                            gcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAg                                       <  4:296808/149‑1 (MQ=255)
                                                                                                            gcgGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCGTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAg                                       <  4:301752/149‑1 (MQ=255)
                                                                                                            gcgGGATTTCGACTACTAAATCATGGTCGGTAACGCGCGCGTGGCAGCTTAAACGGATTTCAGGCTACAGTCCCCAGGCTTTGTAAAGAATGGCGTCTGCCTGCTCTGAGCTTTCCGGGAGTGAGTCAAAACCTTCACGAACGATGCAg                                       <  1:58713/149‑1 (MQ=255)
                                                                                                                       actactAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtggt                            <  7:309933/149‑1 (MQ=255)
                                                                                                                       actactAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGCCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAggtggt                            >  7:367850/1‑149 (MQ=255)
                                                                                                                        ctactaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGgtggtg                           >  1:198498/1‑149 (MQ=255)
                                                                                                                           ctaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGc                          <  7:262929/147‑1 (MQ=255)
                                                                                                                           ctaAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGc                          >  7:129176/1‑147 (MQ=255)
                                                                                                                            taAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTc                                                                         >  6:84567/1‑99 (MQ=255)
                                                                                                                            taAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTc                                                                         <  5:84567/99‑1 (MQ=255)
                                                                                                                            taAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTACGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAg                       >  2:183648/1‑149 (MQ=255)
                                                                                                                                tCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGca                     >  2:290196/1‑147 (MQ=255)
                                                                                                                                        ggTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGAttttt           >  6:68387/1‑149 (MQ=255)
                                                                                                                                          tAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGTCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCTGCAGTGAGTCAAAACCTTCACGAAAGATGAAGTGGAAGGTGGTGCAAGCACAGGATTTTTca         >  3:33947/1‑149 (MQ=255)
                                                                                                                                           aaCGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTca         >  2:258510/1‑148 (MQ=255)
                                                                                                                                               cgcgcCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGc    <  2:231112/149‑1 (MQ=255)
                                                                                                                                                 cgcCTGTCAGCTTAAACGGCTTTCCGGCTCGAGTCACCAGGATTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCgt  <  8:367850/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AGATCGGGATACGCATCGTACAGTGCTTCGCCAATTTCGCGGCTATCGGTCCACTTAAGTCCCATACTAACCTCTGTTAATGCTCACGCGCATGGTTGATAGTGTAACGCGGGATTTCGACTACTAAATCTTCGTCGGTAACGCGCGCCTGGCAGCTTAAACGGCTTTCCGGCTCCAGTCCCCAGGCTTTGTCCAGCATGTCGTCTTCCTGCTCTGAGCTTTCCGGCAGTGAGTCAAAACCTTCACGAACGATGCAGTGGCAGGTGGTGCAAGCACAGGATTTTTCACAGGCGT  >  NC_000913/2656672‑2656965

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