Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,814,762 A→C K299T (AAA→ACA)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,814,7620AC98.6% 228.9 / NA 72K299T (AAA→ACA) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (34/37);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTC  >  NC_000913/1814618‑1814907
                                                                                                                                                |                                                                                                                                                 
acaACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTc                                                                                                                                             <  2:178735/151‑1 (MQ=255)
 caACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCa                                                                                                                                            <  2:843974/151‑1 (MQ=255)
    cGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGCAAAGCCTCACTc                                                                                                                                         >  2:406606/1‑151 (MQ=255)
     gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACt                                                                                                                                          >  1:663161/1‑149 (MQ=255)
     gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTc                                                                                                                                         >  2:119374/1‑150 (MQ=255)
     gTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCg                                                                                                                                        <  1:156267/151‑1 (MQ=255)
        aTGTGGGCGATGTCGGATCGCGTCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGttt                                                                                                                                     <  2:180689/151‑1 (MQ=255)
              gCGATGTCGGATCGCTGCATCCCCCGCTGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATg                                                                                                                               <  2:1144372/151‑1 (MQ=255)
                gATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGaa                                                                                                                             <  1:805681/151‑1 (MQ=255)
                    tCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTGCGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGGTTGGGATGAAGca                                                                                                                          >  2:1202276/1‑150 (MQ=255)
                     cGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGcaca                                                                                                                        <  2:843246/151‑1 (MQ=255)
                      ggATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGaa                                                                                                                             >  2:1201254/1‑145 (MQ=255)
                       gATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGACAGCCTCACACGGTTGGGATGAAACACaa                                                                                                                       >  2:1034806/1‑150 (MQ=255)
                            cgGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACCCCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACt                                                                                                                  <  1:214812/150‑1 (MQ=255)
                            cgGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTc                                                                                                                 <  1:735419/151‑1 (MQ=255)
                               cATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCAcc                                                                                                              >  2:1201477/1‑151 (MQ=255)
                                 tCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCg                                                                                                             <  2:37453/150‑1 (MQ=255)
                                    cccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACg                                                                                                         >  1:283327/1‑151 (MQ=255)
                                    cccGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGGACAGCCTCCCTCGTTTGGGATGAAGCACAAAAACTCACCGGACg                                                                                                         >  2:859629/1‑151 (MQ=255)
                                         aGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGAcc                                                                                                    >  1:389747/1‑151 (MQ=255)
                                             aCCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCg                                                                                                  <  2:479504/149‑1 (MQ=255)
                                                gCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGGGACCCGGACtt                                                                                             >  2:770670/1‑151 (MQ=255)
                                                   ccATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTcc                                                                                           <  2:671156/150‑1 (MQ=255)
                                                       ggAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAccgc                                                                                      <  1:828520/151‑1 (MQ=255)
                                                       ggAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCAccg                                                                                       <  2:600307/150‑1 (MQ=255)
                                                           gttttCGTTATTCACCCCTCCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACTTTTCCACTGGAACCCACTGGCAGTTACGGCCTCACTCTTTTGGGATGAAGCACAAAACCTCACCGGACGTGACCCTGCTTTCCACcgccgct                                                                                  <  2:1041410/148‑2 (MQ=255)
                                                            gCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAAATCACCGGACGTGACCCGGACTTCCACCGCcgcg                                                                                  >  2:822529/1‑150 (MQ=255)
                                                             cTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACcgccgc                                                                                   >  1:460118/1‑148 (MQ=255)
                                                              ttCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGGAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAg                                                                                >  1:267073/1‑150 (MQ=255)
                                                              ttCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGGCGTGACCCGGACTTCCACCGCCGCGAg                                                                                >  2:290532/1‑150 (MQ=255)
                                                                 ggTATTCCCCCCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTTCGTTTCCACTGGAAACCACTGGAAGGTACAGCCTCACTCGTTTGGGATGAGGCACAAAAACCCCCCGGCCGTGACCCGGCCTTCCACCGCCGCGAGTtgt                                                                            <  2:1018965/151‑1 (MQ=255)
                                                                           ccTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGGTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCCCCGCCGAGAGGTGTGGGaaaac                                                                    >  2:512888/1‑146 (MQ=255)
                                                                            cTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCAtt                                                                 >  1:1139165/1‑151 (MQ=255)
                                                                                    tGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAgg                                                         >  1:346314/1‑151 (MQ=255)
                                                                                                 agggaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCCGACGCGACCCGGACTTCCACCGCCGCGAGGTGTGGGAAGCCATTGAAGCAGGCGATTTTCCgg                                              >  2:32810/1‑149 (MQ=255)
                                                                                                  gggaaggCAACTTTTGTACGTTTCCACTGGAAGCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGATTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAAt                                           <  2:485502/151‑1 (MQ=255)
                                                                                                  gggaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGGTGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGCGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGaa                                            >  2:789387/1‑150 (MQ=255)
                                                                                                   ggaaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATa                                          <  1:711156/151‑1 (MQ=255)
                                                                                                     aaggCAACGTTTGTACGTTTCCACTGGCAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACCCACCGGACGAGACCCGGGCTTCCACCGCCGCGGGTTGTGGGAAGCCATTGAAGCAGGCGAATTTCCGGAATACg                                        >  2:336858/1‑151 (MQ=255)
                                                                                                     aaggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACg                                        <  1:876738/151‑1 (MQ=255)
                                                                                                      aggCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATAc                                         >  2:946747/1‑149 (MQ=255)
                                                                                                         cAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAAc                                     >  1:729091/1‑150 (MQ=255)
                                                                                                          aaCGTTTGTGCGTTTCCATTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGAGCCCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTg                                   <  2:772597/151‑1 (MQ=255)
                                                                                                          aaCGTTTGTACGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACt                                    <  2:224183/150‑1 (MQ=255)
                                                                                                            cGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTggg                                 >  2:798321/1‑151 (MQ=255)
                                                                                                            cGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTggg                                 <  1:870870/151‑1 (MQ=255)
                                                                                                            cGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCCCACTCGTTTGGGATGAAGCCCAAAAACTCACCGGACGTGGCCCGGACTTCCACCGCCCCGGGTTGTGGGAAGCCATTGGAACAAGCGATTTTCCGGGAAACCAACTggg                                 >  2:713891/1‑151 (MQ=255)
                                                                                                              tttGTACGTTTCCACTGGAACCCACGGGCAGGTACAGCCTCCCGCGTTTGGGATGAAGCACAAAACCTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTggg                                 <  2:615797/149‑1 (MQ=255)
                                                                                                              tttGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGCGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCt                               >  1:310221/1‑151 (MQ=255)
                                                                                                               ttttACGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCCCCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCt                               <  2:662081/147‑1 (MQ=255)
                                                                                                                tGTACGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAACCTCACCGGACGTGACCCGGACTTCCCCCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCtt                              <  2:180704/150‑1 (MQ=255)
                                                                                                                   aCGTTTCCACTGGAAACCACTGGCAGGTACAGCCTCCCTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTcc                            >  2:291918/1‑149 (MQ=255)
                                                                                                                    cGTTTCCACTGGAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGTACGTGACCGGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAg                          <  2:709229/150‑1 (MQ=255)
                                                                                                                     gTTTCCACTGGAACCCACTGGCAGGTACAGCCTCCCTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAg                          <  2:101298/149‑1 (MQ=255)
                                                                                                                     gTTTCCACTGGAAACCACTGGCAGGTACAGCCTCTCTGGTTTGGGATGAAGCGCAAAATCTCGCCGGACGTGACCCGGACTTCCACCCCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCa                           <  1:223755/148‑1 (MQ=255)
                                                                                                                     gTTTCCACTGGAAACCACTGGCAGGGACAGCCTCACCCGTTTGGGATGAAGCCCAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGGGGTGTGGGAAGCCCTTGAAGCAGGCGATTTTCCGGGATACGAACTGGGCTTCCAg                          >  1:758284/1‑149 (MQ=255)
                                                                                                                      tttCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTTTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGtt                        >  1:538896/1‑150 (MQ=255)
                                                                                                                      tttCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGtt                        <  1:824322/150‑1 (MQ=255)
                                                                                                                      tttCCACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGtt                        <  1:693861/150‑1 (MQ=255)
                                                                                                                         ccACTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGGCATTGAAGCAGGCGATTTTCCGGAATACGAAATGGGGTTCCAGTTGAtt                    >  2:583569/1‑151 (MQ=255)
                                                                                                                           aCTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTcc                  <  1:1192722/151‑1 (MQ=255)
                                                                                                                            cTGGAAACCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTcc                  >  2:1148224/1‑150 (MQ=255)
                                                                                                                              ggAACCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAACCTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTg                <  2:408182/150‑1 (MQ=255)
                                                                                                                                aaaCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTgaag             <  1:1202353/151‑1 (MQ=255)
                                                                                                                                aaaCCACTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTgaag             <  1:1045709/151‑1 (MQ=255)
                                                                                                                                aaaCCACTGGCAGGAACAGCCTCACTCCTTTGGGATGAATCACAAAAACTCACCGGGAGTGACCCGGACTTCCACCGCCGCGACATGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGGTTCCAGATGATTCCTg                >  2:447412/1‑148 (MQ=255)
                                                                                                                                 aaCCACTGGCAGGTACAGCCTCCCTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTgaaga            <  2:351556/151‑1 (MQ=255)
                                                                                                                                    cACTGGCAGGTACAGCCTCCCTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAagaaga         <  1:628431/151‑1 (MQ=255)
                                                                                                                                     aCTGGCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGAt        <  2:926137/151‑1 (MQ=255)
                                                                                                                                        ggCAGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGATGATTCCTGAAGAAGATGa      >  1:1063718/1‑150 (MQ=255)
                                                                                                                                           aGGTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAAtt   >  1:530435/1‑150 (MQ=255)
                                                                                                                                            ggTACAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCCTCCAGTTGATTCCTGAAGAAGATGAATTc  >  1:665291/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACAACGTGATGTGGGCGATGTCGGATCGCGGCATCCCCCGCAGTTACCGCACCATGGAAGGCTTCGGTATTCACACCTTCCGCCTGATTAATGCCGAAGGGAAGGCAACGTTTGTACGTTTCCACTGGAAACCACTGGCAGGTAAAGCCTCACTCGTTTGGGATGAAGCACAAAAACTCACCGGACGTGACCCGGACTTCCACCGCCGCGAGTTGTGGGAAGCCATTGAAGCAGGCGATTTTCCGGAATACGAACTGGGCTTCCAGTTGATTCCTGAAGAAGATGAATTC  >  NC_000913/1814618‑1814907

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