Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947322,1420TC36.6% 54.9 / 64.2 71intergenic (+427/+31)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base C (12/14);  total (33/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 = 6.18e-02
Rejected as consensus: Frequency below/above cutoff threshold.

TGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGCA  >  NC_002947/321994‑322291
                                                                                                                                                    |                                                                                                                                                       
tGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                                        <  1:1322113/150‑1 (MQ=255)
 ggAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcg                                                                                                                                                       <  1:860887/150‑1 (MQ=255)
  gatgCATGAGGGCTAGGACGACCGACCCATGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGCCGCTCATGGCGGTTTGGCACTGCTCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                                      <  1:327681/147‑1 (MQ=255)
    aGCATCAGGGTCAGGACGACCGACCCACGACGTCTTCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGa                                                                                                                                                    <  1:1310440/150‑1 (MQ=255)
     gCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                                   <  1:197006/150‑1 (MQ=255)
       aTCAGGGTCAGGACGACCGGCCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAag‑‑a                                                                                                                                               <  1:851515/150‑1 (MQ=255)
       aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAag‑‑a                                                                                                                                               >  1:1333319/1‑150 (MQ=255)
       aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAag‑‑a                                                                                                                                               <  1:258588/150‑1 (MQ=255)
       aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAag‑‑a                                                                                                                                               >  1:955463/1‑150 (MQ=255)
        tCAGGGTCAGGACGACCGACCCACGACGNNTGCTTGCCGGGCGANANCNCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAag‑‑ag                                                                                                                                              >  1:825805/1‑150 (MQ=255)
          aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGc                                                                                                                                            <  2:955463/150‑1 (MQ=255)
            ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCGGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAg                                                                                                                                          <  2:669440/150‑1 (MQ=255)
                aGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGCAgg                                                                                                                                         >  2:182866/1‑147 (MQ=255)
                aGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGCAgg                                                                                                                                         <  1:182866/147‑1 (MQ=255)
                             cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTaaaa                                                                                                                         >  1:670345/1‑150 (MQ=255)
                                      tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                                   <  2:25146/117‑1 (MQ=255)
                                      tGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                                   >  1:25146/1‑117 (MQ=255)
                                          cgCCGGGCGATACCCCAAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTACACCTGTACCTGCCTCTTCGCCGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAg                                                                                                            <  2:825805/149‑1 (MQ=255)
                                              gggCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                                        <  2:1333319/150‑1 (MQ=255)
                                               ggCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAcc                                                                                                       >  2:628216/1‑150 (MQ=255)
                                                        cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCa                                                                                                                                           <  1:444609/105‑1 (MQ=255)
                                                        cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCa                                                                                                                                           >  2:444609/1‑105 (MQ=255)
                                                        cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAag‑‑a                                                                                                                                               <  2:853889/101‑1 (MQ=255)
                                                        cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAag‑‑a                                                                                                                                               >  1:853759/1‑101 (MQ=255)
                                                        cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAag‑‑a                                                                                                                                               >  1:853889/1‑101 (MQ=255)
                                                        cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAag‑‑a                                                                                                                                               <  2:853759/101‑1 (MQ=255)
                                                         ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACcaca                                                                                             >  2:1129374/1‑150 (MQ=255)
                                                         ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGc                                                                                                                                            <  1:655600/103‑1 (MQ=255)
                                                         ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGc                                                                                                                                            >  2:655600/1‑103 (MQ=255)
                                                           aaGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGCAGCCGGGCCAGT‑‑GGGGAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAg                                                                                           >  2:775723/1‑150 (MQ=255)
                                                              ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCa                                                                                                                      >  2:763072/1‑120 (MQ=255)
                                                              ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCa                                                                                                                      <  1:763072/120‑1 (MQ=255)
                                                                      aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAgg                                                                                                                                         >  2:1324226/1‑93 (MQ=255)
                                                                      aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAgg                                                                                                                                         <  1:1324226/93‑1 (MQ=255)
                                                                               gagaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTa                                                                                       <  1:333634/134‑1 (MQ=255)
                                                                               gagaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTa                                                                                       >  2:333634/1‑134 (MQ=255)
                                                                                     cAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGACTACGCGGATCGCCCACCTGGATCt                                                                 >  2:1038186/1‑150 (MQ=255)
                                                                                        gcgcTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgcg                                                              >  1:88618/1‑150 (MQ=255)
                                                                                                               aTCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtag                                       >  1:1372858/1‑150 (MQ=255)
                                                                                                                tCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTaaaa                                                                                                                         >  2:791204/1‑67 (MQ=255)
                                                                                                                tCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTaaaa                                                                                                                         <  1:791204/67‑1 (MQ=255)
                                                                                                                   aTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTTCGCCGCACTGTTCGCCGATCTGTAGTAgg                                   <  2:1410105/150‑1 (MQ=255)
                                                                                                                          cTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTGACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGATCTGTAGTAGGTCTTGtt                            >  1:298506/1‑149 (MQ=255)
                                                                                                                           tGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                           <  1:628216/150‑1 (MQ=255)
                                                                                                                             tACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGc                         <  2:1372858/150‑1 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                                    >  2:991071/1‑88 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                                    <  1:991071/88‑1 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGAt                                                                                >  1:753443/1‑92 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGAt                                                                                <  2:753443/92‑1 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt                      >  1:950365/1‑150 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                               <  2:85713/141‑1 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                               >  1:85713/1‑141 (MQ=255)
                                                                                                                                cTGCCTCTTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt                      >  1:647625/1‑150 (MQ=255)
                                                                                                                                 tGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGa                                                                                 <  2:314435/90‑1 (MQ=255)
                                                                                                                                 tGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGa                                                                                 >  1:314435/1‑90 (MQ=255)
                                                                                                                                   ccTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGtt                   <  2:950365/150‑1 (MQ=255)
                                                                                                                                    ctctTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTc                                 >  1:970277/1‑135 (MQ=255)
                                                                                                                                    ctctTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTc                                 <  2:970277/135‑1 (MQ=255)
                                                                                                                                    ctctTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtagta                                     >  2:560936/1‑131 (MQ=255)
                                                                                                                                    ctctTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtagta                                     <  1:560936/131‑1 (MQ=255)
                                                                                                                                     tctTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                                        >  1:865076/1‑127 (MQ=255)
                                                                                                                                     tctTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                                        <  2:865076/127‑1 (MQ=255)
                                                                                                                                     tctTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                               <  1:414780/136‑1 (MQ=255)
                                                                                                                                     tctTCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCGACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                               >  2:414780/1‑136 (MQ=255)
                                                                                                                                      ctTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt                <  2:670345/150‑1 (MQ=255)
                                                                                                                                       ttCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                            >  2:1108972/1‑121 (MQ=255)
                                                                                                                                       ttCGCGAGCATGCCCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                            <  1:1108972/121‑1 (MQ=255)
                                                                                                                                         cgcgGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTCTCAGCAATGTTGCTCaa             <  1:50665/150‑1 (MQ=255)
                                                                                                                                          gagggCCTGCCAGCGAAG‑‑TCGCAGGTACAGGCAACTTAAAATCAGAACCTGTAGCTACCCGTAACCACACGGAAACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCAGTAGTAGGTCTTGTCAGCAATGTTGCTCact            <  2:298506/148‑3 (MQ=255)
                                                                                                                                           cgGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATg           <  2:779730/150‑1 (MQ=255)
                                                                                                                                                aTGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgccgcc      >  2:446080/1‑150 (MQ=255)
                                                                                                                                                    tCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGca  >  2:371970/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
TGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAG‑‑AGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGCA  >  NC_002947/321994‑322291

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