Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 641,358 G→T 6.6% intergenic (+40/+81) ahpF → / ← uspG alkyl hydroperoxide reductase, AhpF component/universal stress protein G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913641,3580GT 6.6% 355.2 / 9.5 106intergenic (+40/+81)ahpF/uspGalkyl hydroperoxide reductase, AhpF component/universal stress protein G
Reads supporting (aligned to +/- strand):  ref base G (55/44);  new base T (2/5);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 110 total aligned reads displayed.

CGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATG  >  NC_000913/641209‑641507
                                                                                                                                                     |                                                                                                                                                     
cGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCgggg                                                                                                                                                       <  1:795702/150‑1 (MQ=255)
  ggTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTc                                                                                                                                                     >  2:1106777/1‑150 (MQ=255)
  ggTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTc                                                                                                                                                     >  1:1265213/1‑150 (MQ=255)
       cTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAAcc                                                                                                                                                >  2:668342/1‑150 (MQ=255)
       cTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAAcc                                                                                                                                                >  2:617715/1‑150 (MQ=255)
       cTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAAcc                                                                                                                                                >  2:351432/1‑150 (MQ=255)
       cTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGACGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAAcc                                                                                                                                                >  2:1498413/1‑150 (MQ=255)
         gTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                              <  2:120488/150‑1 (MQ=255)
            cgacgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGGGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTggcg                                                                                                                                           <  1:851011/150‑1 (MQ=255)
              acgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCggcggc                                                                                                                                         <  1:1020972/150‑1 (MQ=255)
              acgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCggcggc                                                                                                                                         <  2:1001243/150‑1 (MQ=255)
              acgGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCAGGTTTGACCCCggcggc                                                                                                                                         <  1:1061306/150‑1 (MQ=255)
                    ccGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCtttttt                                                                                                                                   >  2:1128598/1‑150 (MQ=255)
                       tACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATg                                                                                                                                <  1:1226616/150‑1 (MQ=255)
                       tACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTCTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATg                                                                                                                                >  1:907579/1‑150 (MQ=255)
                       tACAAGCACATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTAtt                                                                                                                                >  2:1460263/1‑149 (MQ=255)
                        aCAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                               <  2:1401308/150‑1 (MQ=255)
                        aCAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                               >  2:119922/1‑150 (MQ=255)
                        aCAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                               <  1:1180591/150‑1 (MQ=255)
                        aCAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                               <  1:220429/150‑1 (MQ=255)
                             cAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTa                                                                                                                          <  2:195218/150‑1 (MQ=255)
                                atcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaaaa                                                                                                                       <  2:887127/150‑1 (MQ=255)
                                 tcatcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTaaaaa                                                                                                                      >  1:143471/1‑150 (MQ=255)
                                  catcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAg                                                                                                                     >  2:945829/1‑150 (MQ=255)
                                    tcatcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGcc                                                                                                                   >  2:1453002/1‑150 (MQ=255)
                                      atcGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGcccc                                                                                                                 >  1:179044/1‑150 (MQ=255)
                                          ccACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCccgccg                                                                                                             >  2:148412/1‑150 (MQ=255)
                                           cACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCgg                                                                                                            <  2:219232/150‑1 (MQ=255)
                                            aCTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCggg                                                                                                           >  1:1004236/1‑150 (MQ=255)
                                            aCTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCggg                                                                                                           <  1:1056055/150‑1 (MQ=255)
                                            aCTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCAGCCggg                                                                                                           <  2:276104/150‑1 (MQ=255)
                                              tGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGTGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGAt                                                                                                         <  1:493993/150‑1 (MQ=255)
                                               ggCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATg                                                                                                        <  1:606845/150‑1 (MQ=255)
                                               ggCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATg                                                                                                        >  2:142787/1‑150 (MQ=255)
                                                 cGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCTCCAAAACTGCATTAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGc                                                                                                      <  2:1052518/150‑1 (MQ=255)
                                                    aGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTcc                                                                                                   >  1:644107/1‑150 (MQ=255)
                                                     ggTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCg                                                                                                  >  2:319279/1‑150 (MQ=255)
                                                       tGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCggg                                                                                                >  2:197930/1‑150 (MQ=255)
                                                          cAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAg                                                                                             >  1:474026/1‑150 (MQ=255)
                                                          cAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGCATGGCTCCGGGCAg                                                                                             >  2:769054/1‑150 (MQ=255)
                                                                ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCggg                                                                                       >  1:930683/1‑150 (MQ=255)
                                                                ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCggg                                                                                       <  1:752106/150‑1 (MQ=255)
                                                                    ctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATAc                                                                                   <  1:148412/150‑1 (MQ=255)
                                                                     tGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACt                                                                                  <  1:653848/150‑1 (MQ=255)
                                                                     tGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACt                                                                                  >  1:743657/1‑150 (MQ=255)
                                                                          gCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTc                                                                             <  1:1510679/150‑1 (MQ=255)
                                                                           cTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCt                                                                            <  2:521139/150‑1 (MQ=255)
                                                                                 aCTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAAt                                                                      >  2:1087785/1‑150 (MQ=255)
                                                                                 aCTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAAt                                                                      <  1:1128598/150‑1 (MQ=255)
                                                                                  cTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACTTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAAtt                                                                     >  1:111324/1‑150 (MQ=255)
                                                                                    aCCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCTATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                                   >  2:208800/1‑150 (MQ=255)
                                                                                    aCCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTaa                                                                   >  1:424851/1‑150 (MQ=255)
                                                                                     ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                                  >  1:1082700/1‑150 (MQ=255)
                                                                                     ccTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                                  <  1:351432/150‑1 (MQ=255)
                                                                                     ccTGATTCACACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAAc                                                                  >  2:1325699/1‑150 (MQ=255)
                                                                                      cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACg                                                                 >  1:1380880/1‑150 (MQ=255)
                                                                                      cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACg                                                                 <  1:945829/150‑1 (MQ=255)
                                                                                      cTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACg                                                                 <  2:1265213/150‑1 (MQ=255)
                                                                                       tGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGc                                                                <  1:448318/150‑1 (MQ=255)
                                                                                        gATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGca                                                               <  1:119922/150‑1 (MQ=255)
                                                                                              cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCa                                                         >  2:802898/1‑150 (MQ=255)
                                                                                              cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCa                                                         <  1:319279/150‑1 (MQ=255)
                                                                                                ccAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGc                                                       >  2:434778/1‑150 (MQ=255)
                                                                                                ccAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGc                                                       >  1:1348711/1‑150 (MQ=255)
                                                                                                 cAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCAGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGCGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACACCCCGcc                                                      >  1:386480/1‑149 (MQ=255)
                                                                                                   aaaCTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAcc                                                    >  1:590161/1‑150 (MQ=255)
                                                                                                       tGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCa                                                <  1:887146/150‑1 (MQ=255)
                                                                                                          aTAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGAt                                             <  2:907579/150‑1 (MQ=255)
                                                                                                           tAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGc                                                                                          <  1:1072425/104‑1 (MQ=255)
                                                                                                              gaagaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggc                                         >  1:6479/1‑150 (MQ=255)
                                                                                                               aagaagTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcg                                        <  2:424851/150‑1 (MQ=255)
                                                                                                                   agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                    >  1:617189/1‑150 (MQ=255)
                                                                                                                   agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                    >  1:1156678/1‑150 (MQ=255)
                                                                                                                   agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                    >  1:694555/1‑150 (MQ=255)
                                                                                                                    gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                   <  1:1087785/150‑1 (MQ=255)
                                                                                                                    gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                   <  1:960376/150‑1 (MQ=255)
                                                                                                                      aaGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGa                                 >  1:934356/1‑150 (MQ=255)
                                                                                                                      aaGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGa                                 <  1:802898/150‑1 (MQ=255)
                                                                                                                       aGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGAt                                >  1:903939/1‑150 (MQ=255)
                                                                                                                        gATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATc                               >  2:314653/1‑150 (MQ=255)
                                                                                                                            cACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGc                           >  1:575106/1‑150 (MQ=255)
                                                                                                                             aCCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCt                          <  1:1325699/150‑1 (MQ=255)
                                                                                                                              ccTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTc                         <  2:1190582/150‑1 (MQ=255)
                                                                                                                              ccTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTc                         >  1:164327/1‑150 (MQ=255)
                                                                                                                               ccgCAATTGCCTACCCGGCGGGGGCAAACCCGGCGGGCTTTTTATTGGATTAAAAACCCCCGCCGGGCAGGCCCCGGGGACGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCg                        <  1:197930/148‑1 (MQ=255)
                                                                                                                                    aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                   <  2:179044/150‑1 (MQ=255)
                                                                                                                                    aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                   <  2:221550/150‑1 (MQ=255)
                                                                                                                                    aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                   <  1:1039710/150‑1 (MQ=255)
                                                                                                                                     ttGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGt                  <  1:1498413/150‑1 (MQ=255)
                                                                                                                                       gCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTa                <  2:743657/150‑1 (MQ=255)
                                                                                                                                       aaaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAcc                                                    >  1:576575/5‑114 (MQ=38)
                                                                                                                                       aaaaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCAcc                                                    <  2:576575/110‑1 (MQ=38)
                                                                                                                                         aaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGa              <  1:1189312/148‑1 (MQ=255)
                                                                                                                                         aaaGCCGCCAGGTTTGACCCCGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGa              >  2:1189312/3‑150 (MQ=255)
                                                                                                                                          tAGCCGCCGGGGTCAAACCTGGCGCCTTTTTTATGGCATTAAAAAGCCCCGCAGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGaa             >  2:651199/1‑150 (MQ=255)
                                                                                                                                              cgccgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTa         >  2:818814/1‑150 (MQ=255)
                                                                                                                                                 cgGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACa      >  1:167104/1‑150 (MQ=255)
                                                                                                                                                  ggggTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAg     <  2:111324/150‑1 (MQ=255)
                                                                                                                                                   gggTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGa    >  1:1289604/1‑150 (MQ=255)
                                                                                                                                                     gTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATg  <  1:1477225/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAGATG  >  NC_000913/641209‑641507

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