Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 641,355 G→A 6.7% intergenic (+37/+84) 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,3550GA 6.7% 321.3 / 9.1 104intergenic (+37/+84)ahpF/uspGalkyl hydroperoxide reductase, AhpF component/universal stress protein G
Reads supporting (aligned to +/- strand):  ref base G (53/44);  new base A (2/5);  total (55/49)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 110 total aligned reads displayed.

AGCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAG  >  NC_000913/641207‑641504
                                                                                                                                                    |                                                                                                                                                     
aGCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCgg                                                                                                                                                      <  2:931798/150‑1 (MQ=255)
 gCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCggg                                                                                                                                                     <  1:1328589/150‑1 (MQ=255)
  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:351432/1‑150 (MQ=255)
         cTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGACGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAAcc                                                                                                                                             >  2:1498413/1‑150 (MQ=255)
           gTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTg                                                                                                                                           <  2:120488/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)
                   gTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCttt                                                                                                                                   <  1:59640/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                                                                                                                            <  1:220429/150‑1 (MQ=255)
                          aCAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                            >  2:119922/1‑150 (MQ=255)
                          aCAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATgg                                                                                                                            <  1:1180591/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)
                                              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                                                                                                     >  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)
                                                                 ccTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCgg                                                                                     <  2:597599/150‑1 (MQ=255)
                                                                  ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCggg                                                                                    <  1:752106/150‑1 (MQ=255)
                                                                  ctctctGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCggg                                                                                    >  1:930683/1‑150 (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                                                                   <  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:839401/1‑150 (MQ=255)
                                                                                         tGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGc                                                             <  1:448318/150‑1 (MQ=255)
                                                                                          gATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGca                                                            <  1:119922/150‑1 (MQ=255)
                                                                                                cACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCa                                                      >  2:802898/1‑150 (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)
                                                                                                            aTAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGAt                                          <  1:294135/150‑1 (MQ=255)
                                                                                                             tAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGc                                                                                       >  2:1072425/1‑104 (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:694555/1‑150 (MQ=255)
                                                                                                                     agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                 >  1:1156678/1‑150 (MQ=255)
                                                                                                                     agTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggc                                 >  1:617189/1‑150 (MQ=255)
                                                                                                                      gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                <  1:960376/150‑1 (MQ=255)
                                                                                                                      gTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATtggcgtggcg                                <  1:1087785/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)
                                                                                                                                  tGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGa                    >  2:209283/1‑150 (MQ=255)
                                                                                                                                      aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                <  2:221550/150‑1 (MQ=255)
                                                                                                                                      aTTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCg                <  2:179044/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)
                                                                                                                                           ttAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGa           >  2:1359275/1‑150 (MQ=255)
                                                                                                                                           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)
                                                                                                                                                    |                                                                                                                                                     
AGCGGGTGACTGTACGACGGTTCCGTACAAGCAGATCATCATCGCCACTGGCGAAGGTGCCAAAGCCTCTCTGAGTGCTTTTGACTACCTGATTCGCACCAAAACTGCATAAGAAGAAGTAAGATTCACCTGCAATTGCTTAGCCGCCGGGGTCAAACCTGGCGGCTTTTTTATGGCATTAAAAAGCCCCGCCGGGATGGCTCCGGGCAGGGCGGGATACTTATTCTGGCAATTAACGCACAACCAGCACCGGCAGATTGGCGTGGCGGATCACGCTCGAGGCGTTAGAACCTAACAG  >  NC_000913/641207‑641504

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: