Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,541,403 G→A D333N (GAC→AAC)  feoB → Fe(2(+)) transporter FeoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,541,4030GA97.0% 196.9 / NA 69D333N (GAC→AAC) feoBFe(2(+)) transporter FeoB
Reads supporting (aligned to +/- strand):  ref base G (2/0);  major base A (27/38);  minor base C (1/1);  total (30/39)
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

TCTGTGATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGTACGCTTG  >  NC_000913/3541170‑3541649
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ccTGTGATGTTGTAAGCAACACCCTGACGGCGGACCCCTGCCTTTTCACCACTGCTGTAGATAAAATCCTGCTCAACCGTTTCTTCGGTCTGCTTATTTGCCTCTTTGTGAGGTACCTGATGTTCCTTCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTTTTTTACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGTACTTGCTACACGATTCACTTCCCGAACTGGCTGACTttcttc                                                                                                                                                                                                                                       <  2:144734/250‑1 (MQ=255)
   gtgATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGTCTGACTATCTTCCta                                                                                                                                                                                                                                    <  2:293063/251‑2 (MQ=255)
   gtgATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCtt                                                                                                                                                                                                                                        <  1:160626/247‑1 (MQ=255)
     gATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGGCTTACAGCCGCCGTTTTTAGGCGGGTCCCAGGCGGTCTTTTTGCCTGGGATTCAATGGATTGGGTCCACGCCCCACTTCCCGAACTGGGTGACTATTTTCCTccc                                                                                                                                                                                                                                  >  2:360650/1‑249 (MQ=255)
     gATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCg                                                                                                                                                                                                                                   <  1:290237/250‑1 (MQ=255)
     gATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCg                                                                                                                                                                                                                                   <  1:281246/250‑1 (MQ=255)
      aTGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGcc                                                                                                                                                                                                                                 <  1:259704/251‑1 (MQ=255)
           gTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGGTACAGCCGCTGTTTTACAGCGGATCCCGGGCGCTATTTGTGCTTGGGATTCAATGGATTGGCTACACGCTCCACCTCCCGAACTGGCTGACCATCTTTCTCGTCCAggg                                                                                                                                                                                                                            >  1:163839/1‑251 (MQ=255)
               ccAACTCCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTg                                                                                                                                                                                                                        <  2:283083/250‑1 (MQ=255)
                   cacCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTg                                                                                                                                                                                                                    <  1:388324/251‑1 (MQ=255)
                    acCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGTGGTACCGGCCCTTATTTAAACCGGCGCCGTGTTGCTATTTGTTCTTGGGATTCAATGGATCGGCTACACTGCTCACTTTCCGAACCGGCTGACTATTCTCCTCGCCCAGAGCCTgggcgg                                                                                                                                                                                                                   >  2:341440/1‑248 (MQ=255)
                    acCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGGGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTTTGCATGGGATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCCGGGGCTgggggg                                                                                                                                                                                                                   >  1:337908/1‑248 (MQ=255)
                            ggCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGGTACAGCCGCCGTTTGACGTCGGCTCCGTGGCGCCATTTGTGCATGGTATTCAATGGATTGGCTCCACCCTCCACTTCCCGAACTGGCTGACCATCTTCCTCACCCCGGGCCTGGGGGGCGGCAtt                                                                                                                                                                                                            >  2:68355/1‑250 (MQ=255)
                             gCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGCTTGACGTCCGCCCCGTGGGGCTATTTGTGCATTGTATTCAAAGGCTTGGCTACACGCTCCACTTCCCGAACCGGATGACTATCTTCCTCCCCCAGGGCCTggggggggggatta                                                                                                                                                                                                           >  2:155354/1‑239 (MQ=255)
                              cAGAACCCCGCCGTTTCGCCACTGCGGTAGATAAAATCGTGCTAAACCGTTTCCTCGGTCTTCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCTCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGGTTACTATCTTCCTCGCCCAGGGCCTTGGTGGCGGCATTAac                                                                                                                                                                                                         <  1:343156/251‑1 (MQ=255)
                                gAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTTACGTCGGCCCCGGGGCGCTATTTGGGCATGGTATTCAATGGGTTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGGGGCGGCATTAacac                                                                                                                                                                                                       >  1:24659/1‑251 (MQ=255)
                                         cGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTCCAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGGATTTAATGGATTGGCTACCCGCCCCACTTCCCCAACTGGCTGACTATCTTCCTCGCCCAGGGGCCGGGGGGCGGGATTAAAACCGTGCTGcc                                                                                                                                                                                              >  2:322976/1‑251 (MQ=255)
                                              caccacTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTg                                                                                                                                                                                          <  2:57043/250‑1 (MQ=255)
                                                    tGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTTTTTGTGCCTGGTATTCAATGGATTGGCTACACGCTCCCCTTCCCGAACTGGCTGACTATTTTCCTCGCCCCGGGCCTGGGTGGCGGCGTTAAAACCGTGCCGCCACTGGGGCCGc                                                                                                                                                                                   >  2:280896/1‑251 (MQ=255)
                                                     gCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCAGATTTTCCTCTTTGTGATGTACCTGATGGTCCTGCTGGCTATCAACATCGGCGGGGCGGTACAGCCACGGGTTGACGTCGGGTCCGGGGCGCTATTTGGGAATGGGACTCAATGGGATGGGTAACCGCTCCACTTCCCGAACTGGCTGACTATCTTCTTAGACCAAGGCATGGGTGGCGGCGTTTACACCGTGCTGCCCCTGGGGcccgc                                                                                                                                                                                  >  2:145038/1‑248 (MQ=255)
                                                        gTAGATAAAATCGTGCTCAACCGATTACTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACACCCGCTGTTTGACGTCGGCCCCCTGGCGCCATTTGTACCTGGTATTTAATGGATTGGGAACACGCTCCCCTTCCCGAACTGGCTGACTATCTTCCTCGCCCCGGGACCCGGTGGCGGCGTTAAAAACGTGGTGCCACTCGTGccccncac                                                                                                                                                                               >  2:380450/1‑245 (MQ=255)
                                                          aGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATTGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTg                                                                                                                                                                             <  1:375661/251‑1 (MQ=255)
                                                          aGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGAtt                                                                                                                                                                              <  1:322976/250‑1 (MQ=255)
                                                            aTAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCCATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGTCGATATTCCTCGCCCCGGGCCTTGGTGGCGGCATTAACACCGTGCTGCCCCTGGGGCCGCAGAtttgg                                                                                                                                                                           >  1:315302/1‑248 (MQ=255)
                                                                          aaCCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAACCGCTGTTTGACGCCGGCTCCGTGGCGCTTTTTGTGCATGGGATTCAATGGATTGGCTACACCCTCCACTTCCCGGACTGGCTGACTATCTTCCCCGCCCCGGGCCTGGGGGGGGGCATTAACACCGTGCTGCCACTTGTGCCGccgctttgcatggtttaccttgt                                                                                                                                                             >  2:129574/1‑229 (MQ=255)
                                                                            ccGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTcc                                                                                                                                                           <  1:415769/251‑1 (MQ=255)
                                                                                    tCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGGGTTACATCCGCTGTTTGACGGCGGCTCCCTGGCGCCGTTTTTTCCGGGTATTCACTTGGTTGTCTTCACGCTCCGCTTTCCGAACTGGGTGGCTATCTGTCTCTCCCCAGGCCTCGGTGGCGGCCTTAACACCCAGCTTCCCCTGGTGCCCCCGATTGGCATGGTGTACCCTTTTCTCtccttc                                                                                                                                                   >  2:285525/1‑251 (MQ=255)
                                                                                       gTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTTACGTCGGCTCCGTGGCGCTATTTGTGCATGGGATTCAAAGGGTTGGCTTCACGCTCCACTTCCCGAACTGGCTGGCTATCTTCCTCGCCCCGGGCCTGGGTGGGGGCATTTACACCGTGCCTCCACCGGGTCCGCAGATTGGCGTGATGTACCCGTTCCTCCccctcctt                                                                                                                                                >  2:155050/1‑251 (MQ=255)
                                                                                                    ccgctTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGt                                                                                                                                   <  2:94286/248‑1 (MQ=255)
                                                                                                       ctTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACATCCGCTGTTTGACGTCTGCTCCCAGGCGATTTTTTGGGCTGGTATTTCATGGATTGGCTTAACGCTCCCACTCCCCGCCTGGCTGACTTTCTTCCTCGCCCCCGGCCTGGGGGGGGGCATTCACCCCGTGCTTCCCCTGGGACCGCCGATTGGCATGCTGTACCTGTACCCCTCCTTCCTTGGTGGCTCCGGGtat                                                                                                                                 >  2:222033/1‑250 (MQ=255)
                                                                                                            tgATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACCTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAAAACCGTGCTGCCCCTGGGGCCGCCGATTGGCATGATGTACCCGTTCCTCTCCTTCCTTGAGGACTCCGGGGATATTgcg                                                                                                                           >  1:318117/1‑251 (MQ=255)
                                                                                                                gTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGCGGCGCTCTTTTTGCATCGGATTCAATGGATTGGCTACACGCTCCCCTTCCCGAACTGGCTGGCTATCTTCCTCCCCCACGGCCTGGGCGGCGGCATTAACCCCGTGCTGCCAATGGTGCCGCAGGTTGGCATGATGTAACTGTTCCTCTCCTTCCTTGAGGACTACgggggtatggggtgtg                                                                                                                       >  1:71468/1‑238 (MQ=255)
                                                                                                                 tACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTTCAGCCGCTGTTTTACGTCGGCCCCGTGGCGCCATTTGTGCCTGGTATTCAATGGATTTGCTACACCCTCCACTTCCCGCACTGGCTGACTATCTTCCTCGCCCAGTGCCTGGGTGGCGGCCCTAACACCCTTCTGCCTCTGGTGCCGCAGATTGGCCTGATTTACCTGTTCCTCTCCTTCCTTTAGGACTCCGGGCATTTGGCGCGTgg                                                                                                                      >  1:71698/1‑250 (MQ=255)
                                                                                                                 tACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGAGTTAAAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTTTGCATGGTATTCAATGGATTGGGTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCCATGGGGCCGCAGAGTGGCATGATGTACCTGTTCCTCTCCTTCCTTTAGGAATCCGGGGATATGGCGAGTgc                                                                                                                      >  1:325321/1‑251 (MQ=255)
                                                                                                                  aCCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGTCGCGTgc                                                                                                                      <  1:324531/250‑1 (MQ=255)
                                                                                                                   ccTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTgcgg                                                                                                                    <  1:88902/251‑1 (MQ=255)
                                                                                                                      gATGTTCCTGCTGGCTTGCAACATCGGCGGGGCGTTCCATCCGCTGTTTGACGTCGGCCCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACACCGGGTATATGGCGCGTGCGGCGt                                                                                                                 <  2:415317/251‑1 (MQ=255)
                                                                                                                      gATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTgcggc                                                                                                                   <  1:115015/249‑1 (MQ=255)
                                                                                                                      gATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGGCGGCTCCGTGGGGGTATTTGTGCAAGGTATTCAATGGATTGGGTACACGCTCCACTTCCCGAACTGGCTGTCTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTTGCATGATGTACCTGTTCCTCTCCTTCATTGAGGAGTCCGGGTATATGGGGCGGgcggcgg                                                                                                                 >  1:365493/1‑250 (MQ=255)
                                                                                                                         gTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGtt                                                                                                                <  2:215795/249‑1 (MQ=255)
                                                                                                                           tCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTtgt                                                                                                             <  1:30105/250‑1 (MQ=255)
                                                                                                                            ccTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGa                                                                                                           <  1:405601/251‑1 (MQ=255)
                                                                                                                            ccTGCTGGCTATCAACATCGGCGGCGAGTTACTGCCGCTGTTTGACGTCGTCTCCGTGGCGCTATTTGTGCATTGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTTCCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTAGATGGCGCGTGCGGCGTTtgtg                                                                                                            <  2:405423/250‑1 (MQ=255)
                                                                                                                               gctgGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTTACGTCGGCTCCGTTGCGCTATTTTTGCATGGTATTTAATGGATTTGCTACCCGCTCCCCTTTCCGAACTGGCTTACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTTGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGGTTGTGATgg                                                                                                        >  1:348461/1‑251 (MQ=255)
                                                                                                                                 tgGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCCGGGCGCCATTTGTGCAAGGTATTCAATGGATTGGCTACACCCCCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGGGGCATTAACACCGTGCTGCCACTGGGGCCGCAGATTGGCCTGATGTAACTGTTCCTCTCCTTCCTTGAGGACTCCGGGGATATGGTGCGTGGGGCGTTTGTGATggcc                                                                                                      >  1:258445/1‑249 (MQ=255)
                                                                                                                                        cAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTg                                                                                                <  2:286824/250‑1 (MQ=255)
                                                                                                                                         aaCATCGGCGGGGGGTTACAGACGCTTTTTTACGTCGGCCCCCTGGCGATATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGGCTATTTTTCTCGCCCCGGGCCTGGGTGGCGGCATTAACACCCTGCTGCCCCTGGGGCCGCATATTGGCCTGAAGTAACTGTTCCTTTCCTTCCTTGATGACTTCGGGTATATGGCGCGGGCGGCGTTTGTGATGGagcgcctgct                                                                                              >  2:278258/1‑242 (MQ=255)
                                                                                                                                                gcggGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGgcg                                                                                       >  1:65318/1‑251 (MQ=255)
                                                                                                                                                  gtggCGTTACAGCCGCTTTTTGACGTCGGCTCCGGGGCGCTATTGGTGGATGGTATTCAATGGATTGGCTACTCGCTCCACTTCCCGAACTGGCTGACTATCTGCCTCGCCCAGGGCCTGGGTGGCGTCATTAACACCGTGCTGCCACTGGTGATGCAGAGTGTCAGGATGTTCTTGTTCCTCTCCTTCCTTGAGGACTACGGGTATATGGCGCGTGCGGCGTTTGTGATGGGCCGTCTGAttgaggagct                                                                                     <  2:45817/249‑10 (MQ=255)
                                                                                                                                                       gTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCt                                                                                <  1:38352/251‑1 (MQ=255)
                                                                                                                                                         tACAGCCGCGGTTTGACTTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACGCCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGGGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCtt                                                                               <  2:198090/250‑1 (MQ=255)
                                                                                                                                                              ccGCTGTTTGACGTCAGCTCCGTGGCGCTATTTGTGCCTGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCggg                                                                         >  1:117703/1‑251 (MQ=255)
                                                                                                                                                                   gTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGaaa                                                                     <  1:313352/250‑1 (MQ=255)
                                                                                                                                                                     ttGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATc                                                                   <  2:348461/250‑1 (MQ=255)
                                                                                                                                                                              gCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCCCCACGTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGTAAATCCTTTGTGCCg                                                         <  2:292205/251‑1 (MQ=255)
                                                                                                                                                                                     ggCGCTATTTGTGCCTGGTATTCAATGGATTGGCTACACGCTCCCCTTCCCGAACTTGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTTCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGAtc                                                   <  2:254121/250‑1 (MQ=255)
                                                                                                                                                                                      gcgcTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGAtcgt                                                 <  1:278258/251‑1 (MQ=255)
                                                                                                                                                                                               tgtgCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGtcggtttc                                          <  1:274252/250‑1 (MQ=255)
                                                                                                                                                                                               tgtgCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGtcggtttc                                          <  1:395501/250‑1 (MQ=255)
                                                                                                                                                                                                gtgCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTTGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGtcggtttcgg                                        <  2:125292/251‑1 (MQ=255)
                                                                                                                                                                                                gtgCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGtcggtttcgg                                        >  1:128795/1‑251 (MQ=255)
                                                                                                                                                                                                            cAATGGATTGGCTACATGCTCCACTTCCCGCACTGGCTTGGTTTCTTCCGTGCCCAGGGCCTTGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTGGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGGTGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGAGCGTCGGTTTCGGTTgtaccgtacc                            <  2:44889/251‑1 (MQ=255)
                                                                                                                                                                                                                          cacGCTCCACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTCCCGTCGGTAATGGGTg              >  1:151722/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ccACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGATGCaagt         >  1:237435/1‑247 (MQ=255)
                                                                                                                                                                                                                                  ccgcTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGt         <  2:262298/247‑1 (MQ=255)
                                                                                                                                                                                                                                  cACTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACCCCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGTAc       <  2:90765/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTTCCCGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGTACg      <  1:280896/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tCCCGAACTGGCTGACTATCTTCATCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCCTGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGTACGCtt   >  2:18281/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cccGAACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGTACGCTTg  <  1:129574/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
TCTGTGATGTGGTAAGCAACACCCTGACGGCAGAACCCAGCCGTTTCACCACTGCGGTAGATAAAATCGTGCTCAACCGTTTCCTCGGTCTGCCGATTTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACACGCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTGCGGCGTTTGTGATGGACCGTCTGATGCAGGCGCTGGGCTTGCCGGGGAAATCCTTTGTGCCGCTGATCGTCGGTTTCGGTTGTAACGTACCGTCGGTAATGGGTGCACGTACGCTTG  >  NC_000913/3541170‑3541649

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