Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,714,238 G→A 7.3% intergenic (+8/+201) srmB → / ← yfiE ATP‑dependent RNA helicase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,714,2380GA 7.3% 258.9 / 8.4 82intergenic (+8/+201)srmB/yfiEATP‑dependent RNA helicase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/43);  new base A (3/3);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

GCTAAACGTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCG  >  NC_000913/2714093‑2714378
                                                                                                                                                 |                                                                                                                                            
gCTAAACGTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCa                                                                                                                                           <  1:140021/149‑1 (MQ=255)
gCTAAACGTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCa                                                                                                                                           <  1:8654/149‑1 (MQ=255)
gCTAAACGTGCTGAGAAGAAAAAAGATAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGTGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGCGCCACCGCAAACGACAGAAGAGAAATCTCAATGCCa                                                                                                                                           >  4:134063/1‑149 (MQ=255)
 cTAAACGTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAg                                                                                                                                          <  2:151159/149‑1 (MQ=255)
       gTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTaa                                                                                                                                    >  1:4879/1‑149 (MQ=255)
         gCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCTACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGc                                                                                                                                  <  4:826687/149‑1 (MQ=255)
           tGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCt                                                                                                                                >  2:96499/1‑149 (MQ=255)
           tGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCt                                                                                                                                >  2:32085/1‑149 (MQ=255)
             agaagaAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTgg                                                                                                                              <  4:320677/149‑1 (MQ=255)
               aagaaAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTa                                                                                                                            >  4:813831/1‑149 (MQ=255)
               aagaaAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTa                                                                                                                            <  4:735423/149‑1 (MQ=255)
                   aaaaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaa                                                                                                                        >  1:337141/1‑149 (MQ=255)
                   aaaaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaa                                                                                                                        <  1:50202/149‑1 (MQ=255)
                   aaaaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAttga                                                                                                                        <  1:32085/149‑3 (MQ=255)
                   aaaaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAttga                                                                                                                        <  1:96499/149‑3 (MQ=255)
                     aacgCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCt                                                                                                                                >  4:284379/4‑139 (MQ=255)
                     aaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCt                                                                                                                                <  3:284379/139‑1 (MQ=255)
                     aaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaaaa                                                                                                                      >  1:104388/1‑149 (MQ=255)
                     aaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaaaa                                                                                                                      <  3:130330/149‑1 (MQ=255)
                     aaaGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTaaaa                                                                                                                      >  4:315320/1‑149 (MQ=255)
                       aGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGt                                                                                                                    >  3:199208/1‑149 (MQ=255)
                         cTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTgg                                                                                                                              >  1:184630/1‑137 (MQ=255)
                         cTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTgg                                                                                                                              <  2:184630/137‑1 (MQ=255)
                         cTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGc                                                                                                                  >  4:436637/1‑149 (MQ=255)
                         cTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGc                                                                                                                  >  4:383973/1‑149 (MQ=255)
                          tAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCt                                                                                                                 >  2:604592/1‑149 (MQ=255)
                           aaagagaaagagAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGAAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTg                                                                                                                <  3:863021/149‑1 (MQ=255)
                            aagagaaagagaAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAGTGCCAGGTTTAAGCCTGGTATTAAAAGTGCTgg                                                                                                               >  4:572500/1‑149 (MQ=255)
                            aagagaaagagaAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTgg                                                                                                               >  2:751556/1‑149 (MQ=255)
                            aagagaaagagaAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTgg                                                                                                               <  3:180745/149‑1 (MQ=255)
                            aagagaaagagaAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTgg                                                                                                               <  1:121619/149‑1 (MQ=255)
                             agagaaagagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGc                                                                                                              >  1:534947/1‑149 (MQ=255)
                               agaaagagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCaa                                                                                                            <  4:45191/149‑1 (MQ=255)
                                 aaagagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAAc                                                                                                          <  3:545981/149‑1 (MQ=255)
                                 aaagagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAAc                                                                                                          <  1:256307/149‑1 (MQ=255)
                                 aaagagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAAc                                                                                                          >  1:766992/1‑149 (MQ=255)
                                   agagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGc                                                                                                        >  4:53019/1‑149 (MQ=255)
                                   agagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAttgaga                                                                                                                      <  2:321376/135‑5 (MQ=255)
                                   agagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATACCAGGCTTAAACCTGGCAttgaga                                                                                                                      >  1:321376/1‑131 (MQ=255)
                                    gagaaGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCa                                                                                                       >  4:819911/1‑149 (MQ=255)
                                           cgcgGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGATTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGc                                                                                                <  2:769189/149‑1 (MQ=255)
                                               ggTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGa                                                                                            <  2:269475/149‑1 (MQ=255)
                                                   aaaaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTaa                                                                                                                                    <  4:647560/105‑1 (MQ=255)
                                                   aaaaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTaa                                                                                                                                    >  3:647560/1‑105 (MQ=255)
                                                    aaaaaCGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATgcgc                                                                                       <  2:174386/149‑1 (MQ=255)
                                                          gCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCt                                                                                 >  3:459625/1‑149 (MQ=255)
                                                           ccATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCtt                                                                                >  2:713303/1‑149 (MQ=255)
                                                                  gACACCAAAAATAGTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGc                                                                         <  4:435413/149‑1 (MQ=255)
                                                                        aaaaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCTCCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                   >  3:805653/1‑149 (MQ=255)
                                                                        aaaaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                   >  1:670546/1‑149 (MQ=255)
                                                                        aaaaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGt                                                                   >  1:250128/1‑149 (MQ=255)
                                                                         aaaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGTCAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTg                                                                  <  4:230862/149‑1 (MQ=255)
                                                                          aaaTATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTgg                                                                 <  1:220148/149‑1 (MQ=255)
                                                                                tGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATg                                                           <  3:819911/149‑1 (MQ=255)
                                                                                    aaGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAat                                                       <  1:778470/149‑1 (MQ=255)
                                                                                     aGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAata                                                      <  1:591634/149‑1 (MQ=255)
                                                                                      gcgcCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTAt                                                                              <  4:125986/124‑1 (MQ=255)
                                                                                      gcgcCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTAt                                                                              >  3:125986/1‑124 (MQ=255)
                                                                                       cgcCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAatata                                                    <  3:568236/149‑1 (MQ=255)
                                                                                          cGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTg                                                 >  4:631214/1‑149 (MQ=255)
                                                                                                ccAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGTTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTg                                           <  2:766992/149‑1 (MQ=255)
                                                                                                 cAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGc                                          >  3:843982/1‑149 (MQ=255)
                                                                                                      ggAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGAt                                     <  4:843982/149‑1 (MQ=255)
                                                                                                                 gCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                          <  2:337141/149‑1 (MQ=255)
                                                                                                                 gCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCg                          <  3:295077/149‑1 (MQ=255)
                                                                                                                  ccACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCgg                         >  3:435467/1‑149 (MQ=255)
                                                                                                                   cACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGa                        <  2:276888/149‑1 (MQ=255)
                                                                                                                     ccGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATa                      >  2:304305/1‑149 (MQ=255)
                                                                                                                      cGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                     <  2:643497/149‑1 (MQ=255)
                                                                                                                      cGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATaa                     <  3:427278/149‑1 (MQ=255)
                                                                                                                       gCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAg                    <  3:636818/149‑1 (MQ=255)
                                                                                                                        cAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAgg                   <  2:1793/149‑1 (MQ=255)
                                                                                                                         aaaCGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGc                  <  4:435467/149‑1 (MQ=255)
                                                                                                                           aCGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGt                <  2:343592/149‑1 (MQ=255)
                                                                                                                            cGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGtt               <  1:713303/149‑1 (MQ=255)
                                                                                                                                agaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcg           >  1:800115/1‑149 (MQ=255)
                                                                                                                                 gaagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgc          >  4:202752/1‑149 (MQ=255)
                                                                                                                                  aagaGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgcc         >  2:598617/1‑149 (MQ=255)
                                                                                                                                    gagTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgc          >  4:800090/1‑146 (MQ=255)
                                                                                                                                    gagTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCAcgc          <  3:800090/146‑1 (MQ=255)
                                                                                                                                       tAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATc    <  3:53019/149‑1 (MQ=255)
                                                                                                                                         aTCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg  <  1:604592/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GCTAAACGTGCTGAGAAGAAAAAAGCTAAAGAGAAAGAGAAGCCGCGGGTGAAAAAACGCCATCGCGACACCAAAAATATTGGTAAGCGCCGTAAACCAAGCGGAACGGGCGTGCCACCGCAAACGACAGAAGAGTAATCTCAATGCCAGGTTTAAGCCTGGTATTAAAAGTGCTGGCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGCTCATGCAATATATTGAATTTGCACGATCTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCG  >  NC_000913/2714093‑2714378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: