Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,628,6250AC31.8% 110.2 / 29.4 63V105G (GTC→GGC) pdeFcyclic di‑GMP phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (39/4);  new base C (0/20);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 7.88e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628487‑2628764
                                                                                                                                           |                                                                                                                                           
cTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAc                                                                                                                                            >  1:1818254/1‑140 (MQ=255)
 tCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGAcc                                                                                                                                           <  2:730870/140‑1 (MQ=255)
     cAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGt                                                                                                                                       >  1:1270897/1‑140 (MQ=255)
      aatAACCCAACAAAAGCAGCAAACGGGAAAAGTATAAGAAACTGCGTTGAAGAAAACACACTCTGCTAAAATAAACGTTGGGAAATTAAACGGG‑TATCACCGTGAGAACCGTTCTTAGCCGGGGGAAGAAAGCCCCTGTa                                                                                                                                      <  2:1067883/138‑1 (MQ=255)
        tAATCCTACAAATGCAGCACACTGGAAAAGTATCAGAAACAGCGATGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAAc                                                                                                                                    >  1:1700156/1‑140 (MQ=255)
        tAACCCTACAAATGCAGCAAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTTGAAAATTAAACGGG‑TACCCCCATGTGAAACGTTGTTCCGCCGGGGAGAAAAGCCCCTGTAAc                                                                                                                                    <  1:506333/140‑1 (MQ=255)
          aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCg                                                                                                                                  >  2:89422/1‑140 (MQ=255)
          aTCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGATGTTACGCCGGGGAGTAAAGACCCTGTAACCg                                                                                                                                  >  1:127099/1‑140 (MQ=255)
          aTCCTAAAAATGCAGCAAACGGGAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCTGGGAAATTAACCGGG‑TATCCCCATGTAAAACTTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCg                                                                                                                                  <  2:224010/140‑1 (MQ=255)
           tCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTTGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGc                                                                                                                                 >  2:506333/1‑140 (MQ=255)
             cTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGccc                                                                                                                               >  1:1643956/1‑140 (MQ=255)
               aCAAATGCAGCAAACTGGAGAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCa                                                                                                                             >  1:1127528/1‑140 (MQ=255)
                aaaaTGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCGGGGAAATTAAGCGGG‑TATCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCaa                                                                                                                            <  2:1736842/139‑1 (MQ=255)
                  aaTGCAGAAAACTGAAAAAGTATCAAAAACAGCGTGGCAGAAAACAAAATCTGCAAAAATATCCGCGGGAAAATAAACCGGG‑TATCCCCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  2:225503/140‑1 (MQ=255)
                  aaTGCAGAAAACTGAAAAAGTATAAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGAAAATAAAACGGG‑TATCACCATGTGAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACa                                                                                                                          <  2:1458109/140‑1 (MQ=255)
                    tGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACaga                                                                                                                        >  1:1647355/1‑140 (MQ=255)
                     gcagcaAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                       <  1:1073473/140‑1 (MQ=255)
                     gcagcaAACGGGAAAAGTATCAAAAACAGCGTGGCAGGAAACACAATCTGCCAGAATATACGCGGGGAAATTAAACGGG‑TATCCCCATGTAAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACagag                                                                                                                       <  2:91297/140‑1 (MQ=255)
                     gcagaaAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAAAATATCCGGGGAAAAATTAAACGGG‑TACCCCCATGTAAACCTTTTTTACCCGGGGGAGAAAACCCCCTGTAACCGCCCCAACagag                                                                                                                       <  2:1224705/140‑1 (MQ=255)
                      cagcaAACTGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAAAATATACGCTGGAAAATTAAACGGG‑TATCACCATGTAAAACGTTGTTACGCCGGGGATAAAAGCCCCTGTAACCGCCCCAACAgaga                                                                                                                      <  2:177109/140‑1 (MQ=255)
                      cagaaAACTGGAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATACGCTGGGAAATTAAACGGG‑TATCCCCATGTGAAACGTTGTTACGCCGGGGAGAAAAGCCCCTGTAACCGCCCCAACAgaga                                                                                                                      <  1:7823/140‑1 (MQ=255)
                        gcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGaat                                                                                                                    >  1:467990/1‑140 (MQ=255)
                              ggaaaaaGTAACAAAAACAGCGTTGCAGAAAACAAAATCGCCCAAAATATACGCGGGAAAATAAAACGGG‑TACCACCATGGAAACCGTTGTTCCGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGaatgataat                                                                                                              <  2:1647355/137‑1 (MQ=255)
                                      tATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGAtaataa                                                                                                      >  2:350861/1‑140 (MQ=255)
                                        tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATaaaa                                                                                                    >  2:1089885/1‑140 (MQ=255)
                                            aagcggCGTGGCAGAAAACACAATCTGCAAAAATATCCGCGGGAAAATTAGGGGGT‑AATCACTGTGAAAATCGTTTTCGGGGGGGGGAGAAAACCCCGTTAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  2:1818254/135‑1 (MQ=255)
                                            aaaCAGCGTTGCAGAAAACAAAATCTGCCAGAATATACGCTGGGAAATAAAACGGG‑TATCACCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  2:1127528/140‑1 (MQ=255)
                                            aaaCAGCGTGGCAGAAAACACATTCTGCAAGAATATCCGC‑GGGAAAATAAACGGGTTACCACCATGTAAAAGTTTTTTACGCGGGGGAAAAAAGCCCCTTTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  2:901843/140‑1 (MQ=255)
                                              aCAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTAAAACGTTTTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                                              <  2:1394005/140‑1 (MQ=255)
                                                aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  1:1561780/1‑140 (MQ=255)
                                                aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt                                                                                            >  2:1502837/1‑140 (MQ=255)
                                                  cGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                                          >  1:894167/1‑140 (MQ=255)
                                                      gCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                >  2:1072/1‑130 (MQ=255)
                                                      gCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                                                <  1:1072/130‑1 (MQ=255)
                                                               acaATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGAAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc                                                                             >  1:830470/1‑140 (MQ=255)
                                                                   tCTGCCAGAATATACGCGGGAAAATAAAAGGGG‑TACCCCCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa                                                                         <  2:225301/140‑1 (MQ=255)
                                                                        cAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                    >  1:1654248/1‑140 (MQ=255)
                                                                        cAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGAGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                    >  1:402423/1‑140 (MQ=255)
                                                                        aaGAATATCCGCGGGAAAATAAAACGGG‑TATCACCATGTAAAACGTTGTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                                                    <  2:160546/139‑1 (MQ=255)
                                                                         aGAATATCCGCTGGAAAATAAACCGGG‑TATCACCATGTAAAACGTTTTTACGCGGGGGAGAAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                                   <  2:1007362/140‑1 (MQ=255)
                                                                           aaTATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                                 >  2:827844/1‑140 (MQ=255)
                                                                               taCGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                                             >  2:1389091/1‑140 (MQ=255)
                                                                                 cGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                                           >  2:342383/1‑140 (MQ=255)
                                                                                     gggAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  2:1305907/1‑140 (MQ=255)
                                                                                     gggAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  1:330568/1‑140 (MQ=255)
                                                                                     gggAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                                                       >  1:942459/1‑140 (MQ=255)
                                                                                      ggAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt                                                      >  1:1678358/1‑140 (MQ=255)
                                                                                          atcaaaGGGG‑TATCGCAATGGGAACCTTGGTTCCGCGGGGGGGAAAAGACCTTCTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                                                  <  1:1526183/137‑1 (MQ=255)
                                                                                           ttAAACGGG‑TATCACCATGTGAAACGTTGTTACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCgg                                                 <  2:1181044/140‑1 (MQ=255)
                                                                                                 ggg‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaa                                           <  1:89422/140‑1 (MQ=255)
                                                                                                  gg‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaaa                                          >  1:514961/1‑140 (MQ=255)
                                                                                                   g‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                                         >  1:542634/1‑140 (MQ=255)
                                                                                                        gaCCATGTGAAACTTTGTTACGCCGGGGAGTAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                                     <  2:1373473/139‑1 (MQ=255)
                                                                                                              gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               >  2:630388/1‑140 (MQ=255)
                                                                                                              gtaaagcGTTGTTGCGCGGGGGAGAAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg                               <  2:1657750/134‑1 (MQ=255)
                                                                                                               tgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa                              >  2:826193/1‑140 (MQ=255)
                                                                                                                 aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc                            >  1:1440134/1‑140 (MQ=255)
                                                                                                                    cGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcccc                         >  1:1373473/1‑139 (MQ=255)
                                                                                                                      ttgttCCGCGGGGGAGAAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       <  2:894167/140‑1 (MQ=255)
                                                                                                                      ttgttACGCGGGGGAGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       <  1:207975/140‑1 (MQ=255)
                                                                                                                      ttgttACGCGGGGGAGAAAAGACCCCGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       <  2:330568/140‑1 (MQ=255)
                                                                                                                      ttgtgACGCGGGGAATTAAAGCCCTTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       <  2:942459/140‑1 (MQ=255)
                                                                                                                           aCGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATg                  >  2:694389/1‑140 (MQ=255)
                                                                                                                              ccGGGAAGTAAAGCCCCTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATgccg               <  2:1561780/140‑1 (MQ=255)
                                                                                                                                 gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  2:1001424/1‑140 (MQ=255)
                                                                                                                                 gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  1:383020/1‑140 (MQ=255)
                                                                                                                                 gggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  2:1218267/1‑140 (MQ=255)
                                                                                                                                 gggAGTAAAACCCGTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCCTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCCGACAAAAATCCCTAAGGCGATCCCTGGCAATGCGGccc            <  1:1009665/140‑1 (MQ=255)
                                                                                                                                  ggAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCa           <  1:1305907/140‑1 (MQ=255)
                                                                                                                                   gaaaaaaGCCCTTGAAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAg          <  2:823181/136‑1 (MQ=255)
                                                                                                                                    aGTAAAGCCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       <  2:129272/126‑1 (MQ=255)
                                                                                                                                    aGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTgg                       >  1:129272/1‑126 (MQ=255)
                                                                                                                                        aaGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa     >  2:565386/1‑140 (MQ=255)
                                                                                                                                         aGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaa    >  1:823181/1‑140 (MQ=255)
                                                                                                                                           aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            >  1:972719/1‑130 (MQ=255)
                                                                                                                                           aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc            <  2:972719/130‑1 (MQ=255)
                                                                                                                                           aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  1:600089/1‑140 (MQ=255)
                                                                                                                                           aCCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaaaaa  >  1:866755/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CTCGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAAAAA  >  NC_000913/2628487‑2628764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑