Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,710 C→T 7.1% intergenic (+27/+18) fadH → / ← higA 2,4‑dienoyl‑CoA reductase/antitoxin/DNA‑binding transcriptional repressor HigA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,7100CT 7.1% 257.0 / 6.9 85intergenic (+27/+18)fadH/higA2,4‑dienoyl‑CoA reductase/antitoxin/DNA‑binding transcriptional repressor HigA
Reads supporting (aligned to +/- strand):  ref base C (44/34);  new base T (2/4);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATC  >  NC_000913/3233562‑3233856
                                                                                                                                                    |                                                                                                                                                  
tGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACg                                                                                                                                                   <  2:639616/150‑1 (MQ=255)
tGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACg                                                                                                                                                   <  1:1318768/150‑1 (MQ=255)
 ggCGCAACCGCTGATTGAGAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCTATGTGGCTATGGAGCTGGACGCACGACTGGCAATTGCCCAGGGAACACGGCTGGCCCTGGAGATTTACATCGTTATTTGAAGCGCCGGATGCAACGc                                                                                                                                                  >  2:888892/1‑150 (MQ=255)
     cAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATcc                                                                                                                                              <  1:846705/150‑1 (MQ=255)
       aCCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCgg                                                                                                                                            >  1:965769/1‑150 (MQ=255)
        ccGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGc                                                                                                                                           <  1:723544/150‑1 (MQ=255)
         cGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCa                                                                                                                                          <  2:588676/150‑1 (MQ=255)
         cGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATCGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCa                                                                                                                                          <  2:722835/150‑1 (MQ=255)
           cTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGTGCCGGATGCGTTGCATCCGGcgct                                                                                                                                        <  1:253459/150‑4 (MQ=255)
           cTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGcgct                                                                                                                                        <  1:424331/150‑4 (MQ=255)
           cTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCGTTGCATCCGGcgct                                                                                                                                        <  2:302655/150‑4 (MQ=255)
                    aGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGAAATTGCCCAGGGAACACGGCGGGCGCTGGAGATTTAAATCGTTATTTGACGCGCCGGATGCAACGCATCCGGCACGTCGCGTTaa                                                                                                                               <  1:903489/150‑1 (MQ=255)
                     gCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAAt                                                                                                                              >  2:216294/1‑150 (MQ=255)
                     gCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAAt                                                                                                                              >  1:1338917/1‑150 (MQ=255)
                        ggAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCaa                                                                                                                           >  2:269511/1‑150 (MQ=255)
                          aaaaCGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATa                                                                                                                         <  1:907367/150‑1 (MQ=255)
                          aaaaCGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTCCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATa                                                                                                                         >  2:327612/1‑150 (MQ=255)
                                    aTTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCgg                                                                                                               >  1:168029/1‑150 (MQ=255)
                                     tttAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAAGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCggg                                                                                                              >  1:1302683/1‑150 (MQ=255)
                                       tAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGgag                                                                                                            >  1:546890/1‑150 (MQ=255)
                                       tAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGAACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGCCACGTAGCGTTAATCAATAAACAAGGCGGgag                                                                                                            >  2:1264996/1‑150 (MQ=255)
                                            ggcggcggcGATGTGGCTATGGAGCGGCACGCACGACGGGAAATTGCCCAGGCAACACGGCGGGCGCGGGATATTTAAATCGTTATTTCACGCGCCGTATGCACCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATg                                                                                                       <  2:2554/150‑1 (MQ=255)
                                            ggcggcTGCGATGTGGCTATGGAGCTGGACGCACGAGGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTTAAGCGCGGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATg                                                                                                       <  2:885032/150‑1 (MQ=255)
                                              cggcTGCGATGTGGCTATGGCGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGGGCGGTAGATTTAAATCGTTATTTCAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGcc                                                                                                     <  2:282065/150‑1 (MQ=255)
                                                gctgctATGTGGCTATGTAGCTGTACGCACTACGGGAAATTGCCCTGGTAACACGGCTGGCGCTGGATATTTAAATCTTTTTTTGAAGCGCCGTATGAAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGa                                                                                                   <  1:284174/150‑1 (MQ=255)
                                                  tGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAAt                                                                                                 >  1:1223217/1‑150 (MQ=255)
                                                    cGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGATATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcg                                                                                               <  2:995182/150‑1 (MQ=255)
                                                    cGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcg                                                                                               >  2:186035/1‑150 (MQ=255)
                                                      aggtgGCTATGGAGCGGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCGGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATcgcg                                                                                             <  1:1101644/148‑1 (MQ=255)
                                                       tgtgGCTATGCAGCTGGACGCACGACGGCCAATTGCCCAGGGAACCCGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGt                                                                                            <  1:34542/150‑1 (MQ=255)
                                                       ggtgGCTTTGTTGCTGGACGCACGAGGGGAAATTGCCCTGGTAACAGGGCGGGCGGTGTATATTTAAATGGTTATTTGATGCGCGGGATGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGt                                                                                            <  1:1014829/149‑1 (MQ=255)
                                                             tATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCaa                                                                                      >  2:1112351/1‑150 (MQ=255)
                                                              aTGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAt                                                                                     >  1:777051/1‑150 (MQ=255)
                                                               tGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATATGAAGCGCCGGATGCAACGCATCCCGCACGTCGCGTTAATCAATAAACAAGCCGGGAGAAATGCCGAATCGCCTTGCCAAtt                                                                                    >  2:751958/1‑150 (MQ=255)
                                                                ggAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAttt                                                                                   >  2:135852/1‑150 (MQ=255)
                                                                   gCTGGCCGCACGCCGGCCAATGGCCCCGGGAACCCGGCGGGCGCTGCTTTTTAAAATCGTTATTTGAAGCGCGGTATGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAAtttttt                                                                                <  1:966817/150‑1 (MQ=255)
                                                                     tGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAg                                                                              <  2:898344/150‑1 (MQ=255)
                                                                     tGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAg                                                                              <  2:1338917/150‑1 (MQ=255)
                                                                      ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATGAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                             >  1:183865/1‑150 (MQ=255)
                                                                      ggACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGc                                                                             >  2:170236/1‑150 (MQ=255)
                                                                       gACGAACGACGGGCAATTGCCCAGGTAACAGTGCTGGCGGTGGATATTTAAATCGTTATTTGAACCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCg                                                                            <  2:557766/150‑1 (MQ=255)
                                                                           cACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGtt                                                                        >  1:937986/1‑150 (MQ=255)
                                                                             cgacgGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTcc                                                                      >  2:1487303/1‑150 (MQ=255)
                                                                               acgGGCAATTGCCCAGGTAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAg                                                                    <  1:94915/150‑1 (MQ=255)
                                                                                 gggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCg                                                                  >  1:1533966/1‑150 (MQ=255)
                                                                                  ggCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGt                                                                 >  1:96900/1‑150 (MQ=255)
                                                                                    cAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTa                                                               >  1:351027/1‑150 (MQ=255)
                                                                                        tGCCCTGGTAACACGGCTGGCGCTGGAGATTTAAATCTTTATTTGAAGCGCCGTATGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAAttt                                                           <  1:1069894/150‑1 (MQ=255)
                                                                                          cccAGGGAACCCGGCTGGCGCTGGATATTTAAATCGTTATTTGAAGCGGGGTTTGCAAGGCATCCGGCACGTCGCGTTAATCAATAAACAATGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         <  1:216294/150‑1 (MQ=255)
                                                                                          cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAGCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         >  2:438098/1‑150 (MQ=255)
                                                                                          cccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         >  2:758029/1‑150 (MQ=255)
                                                                                          cccAGGGAACACGGCGGCCGCTGGAGATTTAAATCGTTATTTGACGCGCCGGATGAAACGCAGCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTcc                                                         <  1:785638/150‑1 (MQ=255)
                                                                                           ccAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCct                                                        >  2:4721/1‑150 (MQ=255)
                                                                                              gggAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCtt                                                     >  2:785118/1‑150 (MQ=255)
                                                                                                   cacGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGc                                                >  2:335460/1‑150 (MQ=255)
                                                                                                    acGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCt                                               <  1:917652/150‑1 (MQ=255)
                                                                                                     cGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTc                                              >  1:737667/1‑150 (MQ=255)
                                                                                                       gcgggCGCTGGAGATTTAAATGTTTTTTTCAAGCGCCGTATCAAACCCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCaa                                            <  1:1010292/147‑1 (MQ=255)
                                                                                                          ggCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAc                                         >  1:770018/1‑150 (MQ=255)
                                                                                                            cgcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgc                                       >  2:718060/1‑150 (MQ=255)
                                                                                                             gcTGGCGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      <  2:275002/150‑1 (MQ=255)
                                                                                                             gcTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                      <  2:1165498/150‑1 (MQ=255)
                                                                                                               tgcagaTTTAAATCGTTATTTGACGCGCGGGATGCACCGCACCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                    <  1:202200/147‑1 (MQ=255)
                                                                                                               tGGAGATTAAAAGCGTTATTTGACGCGCCGTAGGCAACGCATCCGTCACGTCGCGTTAATCAATAAACAAGGCTGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTTTTCCAGCGTTAATTTCGTGTTTCCGCTCATAACGGGTg                                    <  2:73207/150‑1 (MQ=255)
                                                                                                                  agaTTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGaca                                 >  2:113520/1‑150 (MQ=255)
                                                                                                                   gaTTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGacac                                >  1:1473680/1‑150 (MQ=255)
                                                                                                                      ttAATTCGTTTTTTGAAGCGCCGTATGCAACGCATCGGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCAt                             <  2:1223217/150‑1 (MQ=255)
                                                                                                                        aaaTCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAg                           >  2:1435756/1‑150 (MQ=255)
                                                                                                                             gTTATTTGAAGCGCCAGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTa                      >  1:1478827/1‑150 (MQ=255)
                                                                                                                              ttATTTTAAGCGCCGGTTGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAATGCGGGAGAAATGCCGAATCGCGTTGGCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTAc                     <  2:1170929/150‑1 (MQ=255)
                                                                                                                               tATTTCAAGCGCCGTATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACt                    <  1:1210082/150‑1 (MQ=255)
                                                                                                                                aTTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTg                   >  2:965474/1‑150 (MQ=255)
                                                                                                                                aTTTGAAGCGCCGGATGCAACGCATCCGGCACGGCGCGTTAATCAATAAACCAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCGCAAAACGCGTGACACCATAGATTTACTg                   <  2:1060850/150‑1 (MQ=255)
                                                                                                                                    gAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCaa               >  2:1032038/1‑150 (MQ=255)
                                                                                                                                      cgtgcCGGATGCGTTGCATCCGGCGCGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAAtt             >  1:699148/4‑150 (MQ=255)
                                                                                                                                         gcgggttGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGCTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTcc          <  1:722220/144‑1 (MQ=255)
                                                                                                                                         gcCGGATGCGTTGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTcc          >  1:459514/1‑150 (MQ=255)
                                                                                                                                         gcCGGATGCAGCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTcc          <  1:438098/150‑1 (MQ=255)
                                                                                                                                            ggATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc       >  2:567260/1‑150 (MQ=255)
                                                                                                                                            ggATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc       >  2:972352/1‑150 (MQ=255)
                                                                                                                                            ggATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc       >  2:1324925/1‑150 (MQ=255)
                                                                                                                                               tGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGa    >  2:38196/1‑150 (MQ=255)
                                                                                                                                               tGCAACCCAACCTGCACGAAGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATAGCGTTGCCAATATATTAGAGTGTTCCAGCGGTAATTTCCTCTTCCCGCTCAAAACGCGTGAAACCATAGATTTACTGACAATTTCCGGCAGa    >  2:1252337/1‑150 (MQ=255)
                                                                                                                                                 cAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATc  >  1:683369/1‑150 (MQ=255)
                                                                                                                                                 cAACGCATCCGGCACGTCGCGTTAATAAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATc  >  2:101585/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCGCAACCGCTGATTGATAGCGGGAAAACGGTGCATTTAATTGGCGGCTGCGATGTGGCTATGGAGCTGGACGCACGACGGGCAATTGCCCAGGGAACACGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATC  >  NC_000913/3233562‑3233856

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