Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,026,683 G→A 9.1% D445D (GAC→GAT) ‡ dgcQ ← putative diguanylate cyclase DgcQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,026,6830GA 9.1% 192.8 / 7.4 66D445D (GAC→GAT) ‡dgcQputative diguanylate cyclase DgcQ
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base A (2/4);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

CCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCA  >  NC_000913/2026535‑2026832
                                                                                                                                                    |                                                                                                                                                     
ccGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGt                                                                                                                                                      <  1:758174/150‑1 (MQ=255)
 cGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTc                                                                                                                                                     <  1:169261/150‑1 (MQ=255)
  gCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCa                                                                                                                                                    <  1:784668/150‑1 (MQ=255)
     gCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTa                                                                                                                                                 >  1:296439/1‑150 (MQ=255)
        tCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATc                                                                                                                                              <  1:954402/150‑1 (MQ=255)
           gTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCt                                                                                                                                           <  1:473552/150‑1 (MQ=255)
             cAGCCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCttt                                                                                                                                         >  1:432734/1‑150 (MQ=255)
             cAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGTCCAAAGCGGTCATTAATCGCttt                                                                                                                                         <  1:643273/150‑1 (MQ=255)
                aCTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGATTAATGACCGCtttggt                                                                                                                                      <  1:81321/150‑4 (MQ=255)
                   cgcgCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCGTGCGCACGCAAGGAACTGCTAATTAATCCGGCATCATGAGAAAGAACACGGTCGCCCGCCTGATGAGCAAAGCGATTAATGACCGCtttggtcat                                                                                                                                   <  2:286678/150‑7 (MQ=255)
                   cgcgCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATg                                                                                                                                   >  2:858869/1‑150 (MQ=255)
                      gcCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTc                                                                                                                                <  2:1097265/150‑1 (MQ=255)
                       ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGATTAATGACCGCtttggt                                                                                                                                      <  1:198645/143‑4 (MQ=255)
                       ccTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGATTAATGACCGCtttggt                                                                                                                                      >  2:198645/1‑140 (MQ=255)
                        cTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCaa                                                                                                                              >  1:450503/1‑150 (MQ=255)
                                aTCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGAcc                                                                                                                      >  2:456801/1‑150 (MQ=255)
                                  cacacaAAACTCGTCACCACAGACCCGCGCGGCAACGTGCTGCGCACGGAAGGAACTGCTATTTAATCCGGCAGCATGAGAATGAACACTGTCGCCCGCCTGATGACCTAAGCGGTCATTAATCGGTTTAAAATGGTGAAGATCGACGTg                                                                                                                    <  2:447369/150‑1 (MQ=255)
                                  cacacaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTAAAGATCGACCTg                                                                                                                    >  2:716203/1‑150 (MQ=255)
                                   acacaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGa                                                                                                                   >  2:361506/1‑150 (MQ=255)
                                   acacaAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGa                                                                                                                   >  1:1217292/1‑150 (MQ=255)
                                      caAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGATTAATGACCGCtttgg                                                                                                                                       <  2:91078/127‑3 (MQ=255)
                                      caAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGATTAATGACCGCtttgg                                                                                                                                       >  1:91078/1‑125 (MQ=255)
                                            tcctcACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGaa                                                                                                          <  1:502556/150‑1 (MQ=255)
                                                 accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAgg                                                                                                     <  2:764072/150‑1 (MQ=255)
                                                 accaccGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAgg                                                                                                     <  2:1273830/150‑1 (MQ=255)
                                                    accGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATg                                                                                                  <  2:332730/150‑1 (MQ=255)
                                                       gACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTg                                                                                               <  2:1340065/150‑1 (MQ=255)
                                                         cccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGt                                                                                             >  2:102550/1‑150 (MQ=255)
                                                          ccgcccgGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTg                                                                                            <  1:827087/150‑1 (MQ=255)
                                                                ggCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTg                                                                                      <  2:880092/150‑1 (MQ=255)
                                                                 gCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGa                                                                                     >  2:1228578/1‑150 (MQ=255)
                                                                 gCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGa                                                                                     >  2:1027487/1‑150 (MQ=255)
                                                                  cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                    >  1:810925/1‑150 (MQ=255)
                                                                  cAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGac                                                                                    <  1:456801/150‑1 (MQ=255)
                                                                       tCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAt                                                                               <  1:1228578/150‑1 (MQ=255)
                                                                        ccTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAtt                                                                              >  1:1299885/1‑150 (MQ=255)
                                                                         cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAttt                                                                             <  2:553593/150‑1 (MQ=255)
                                                                         cTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAAttt                                                                             >  1:947013/1‑150 (MQ=255)
                                                                            cgcacgcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCTCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagc                                                                          <  2:153587/150‑1 (MQ=255)
                                                                                 gcaAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTagcgagcg                                                                     >  2:434280/1‑150 (MQ=255)
                                                                                      gAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACgg                                                                <  1:306934/150‑1 (MQ=255)
                                                                                       aaCTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACggg                                                               <  2:1088023/150‑1 (MQ=255)
                                                                                          tGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCtt                                                            >  1:407098/1‑150 (MQ=255)
                                                                                               attaatCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCg                                                       <  2:1237776/150‑1 (MQ=255)
                                                                                                       ggCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTgc                                               >  1:614865/1‑150 (MQ=255)
                                                                                                           gcaTGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCa                                           >  1:905317/1‑150 (MQ=255)
                                                                                                              tGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGCCGGGCTTTTTCGAACAGTGCGCCACgg                                        >  1:728180/1‑150 (MQ=255)
                                                                                                                 gaaagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTta                                     <  2:450503/150‑1 (MQ=255)
                                                                                                                 gaaagaaCACGGGCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTta                                     <  1:361506/150‑1 (MQ=255)
                                                                                                                   aagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtata                                   <  1:1027487/150‑1 (MQ=255)
                                                                                                                   aagaaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtata                                   >  2:1287395/1‑150 (MQ=255)
                                                                                                                     caaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTtatata                                 <  1:438894/149‑1 (MQ=255)
                                                                                                                      aaCACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATaa                                >  2:230084/1‑150 (MQ=255)
                                                                                                                        cacGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAc                              >  1:233441/1‑150 (MQ=255)
                                                                                                                           ggTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAAcgcg                           <  1:167827/150‑1 (MQ=255)
                                                                                                                                 ccgccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAgg                     <  1:858869/150‑1 (MQ=255)
                                                                                                                                  cgccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgt                    >  2:1161424/1‑150 (MQ=255)
                                                                                                                                   gccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtg                   >  2:669717/1‑150 (MQ=255)
                                                                                                                                    ccTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGgtgt                  <  2:947013/150‑1 (MQ=255)
                                                                                                                                      tgatgaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCg                >  1:241211/1‑150 (MQ=255)
                                                                                                                                       gatgaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGt               <  1:669717/150‑1 (MQ=255)
                                                                                                                                        atgaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTg              <  2:432734/150‑1 (MQ=255)
                                                                                                                                         tgaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCTAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGc             <  1:308197/150‑1 (MQ=255)
                                                                                                                                          gaCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGcc            <  2:722968/150‑1 (MQ=255)
                                                                                                                                           aCCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCa           >  1:700601/1‑150 (MQ=255)
                                                                                                                                               aaGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGcc       <  2:296439/150‑1 (MQ=255)
                                                                                                                                                   ggTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGcc   >  1:1181339/1‑150 (MQ=255)
                                                                                                                                                    gTCATTAATCGCTTAAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCa  >  2:822678/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGCAGCCTCCGTCAGACTCGCGCCTGGCAGAATCACACAAAACTCCTCACCACCGACCCGCCCGGCAACGTCCTGCGCACGCAAGGAACTGCTAATTAATCCGGCAGCATGAGAAAGAACACGGTCGCCCGCCTGATGACCAAAGCGGTCATTAATCGCTTTAAAATGGTCAAGATCGACCTGAATGACAGAAAAAGGATGTTGGTGCGTCTGACACAATTTAGCGAGCGGACGGGCTTTTTCGAACAGTGCGCCACGGTTATATAAACGCGTTAAGGTGTCGTGCCACGCCTGCCA  >  NC_000913/2026535‑2026832

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: ‑