Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,6040AT32.6% 148.8 / 11.1 96intergenic (‑94/‑281)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (42/18);  new base T (0/29);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 3.15e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624467‑624737
                                                                                                                                         |                                                                                                                                     
ccTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTAACAAATCAGCTTCATGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                     <  1:1332450/140‑1 (MQ=255)
    ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGTGCGTAATGACAAATTCCACAAAGCGCACAATCAGTCCCCTCGCCACGTTGTGGAGAGGGTTAGGGTGa                                                                                                                                 >  2:930895/1‑140 (MQ=255)
    ttAATAGAAGGGCGTTCCGGTAGAGCGGGGCGAGTCTCACAACTCAGCTTCCTGTTATTATTAAGGTTAGGGGCGTAAGGACTAGTTCGACCTAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTCGGGTgt                                                                                                                                 >  2:928290/1‑139 (MQ=255)
         agtcgggCGTTGCGGTAGAGGGGGGCGTGTCACACAAATCTGATTCCTGTTTTTAATATGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGGTGTGGGTTTGGGTGAGGGGa                                                                                                                            <  2:910142/136‑1 (MQ=255)
         aGATGGGCGTTGCGGTTGTGCGGGGCTTGTCTCACAAATCAGCTTCCTGTTTTTAATATGTTTAAGGGCGTATTGAAAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGGTTTGGGTTAGGGGa                                                                                                                            <  2:1498099/140‑1 (MQ=37)
         aGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                            >  1:883487/1‑140 (MQ=255)
         aGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGGGAAGTGCGTACTGACAAATTCGACAAAGCGCACTATCCGACCCCTCGCCCCTATGGGGAGCGGGTAAGGGTGAGGGGa                                                                                                                            >  1:1808479/1‑140 (MQ=255)
              ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                       >  1:681529/1‑140 (MQ=255)
              ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                       >  1:264506/1‑140 (MQ=255)
                cGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAAACCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccc                                                                                                                     >  1:1425543/1‑139 (MQ=255)
                           tggggggggTGTCTCACAACTCTGCTGCGTGTTTTTAATATGGTTTGGGGGTTAATGACAAATTCGACAACGCGCACAATCGGTCCCCGCGCCGGTTGGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCg                                                                                                          <  1:1073972/137‑1 (MQ=255)
                           tgCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATATGGTTATGGGCGAAATGACAAATTCGACAAAGCGCACAATCCGACCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                          <  1:815512/139‑1 (MQ=38)
                           cgCGGGGCGTGTCTCACAAATCCGCTTCCTGTTCTTAATACGGTTACGGGCGTAATGACAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCg                                                                                                          <  2:1375231/139‑1 (MQ=37)
                           agCGGGGCTAGTCTCACAAATCTGCTTCATTTTATTAATACGTATAAGGGCTTAATCACAAATTCGACAACGCGCACATTCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                          <  1:369413/140‑1 (MQ=25)
                           agCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                          <  2:148067/140‑1 (MQ=255)
                             cGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTAGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGaa                                                                                                        >  1:1121040/1‑140 (MQ=255)
                             cGGGGCGAGTCTCACAAATCAGCTTCCGGTTATTAATAAGGTTAAGGGCGTAATGTCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGaa                                                                                                        >  1:1120635/1‑140 (MQ=255)
                               gggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACa                                                                                                      >  2:812583/1‑140 (MQ=255)
                               gggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACa                                                                                                      >  1:427224/1‑140 (MQ=255)
                                ggCGTGTCTCAAAAATCAGCTTCCTGTTATTAATATGGTTATGGGCGTACTGAAAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACAt                                                                                                     <  1:499779/140‑1 (MQ=37)
                                ggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAt                                                                                                     >  1:295129/1‑140 (MQ=255)
                                     gTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAcc                                                                                                >  2:1139384/1‑140 (MQ=255)
                                       ctcACAATTCAGCTTCCTGTTTTTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:872039/140‑1 (MQ=255)
                                       ctcACAATTCAGCTTCCTGTTATTAATAAGTTTAAGGGCGTAATGACAATTTCGTCAAAGGGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:735373/140‑1 (MQ=37)
                                       cacaaaaaTCAGATTCCTTTTATTAATACGGTTTTGGGGTTATTGACAATTTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTGGGGGTGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:1518094/135‑1 (MQ=255)
                                          acaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAc                                                                                           >  2:1783851/1‑140 (MQ=255)
                                            aaaaCAGTTCCCTGTTATTAATACGGTTACGGGCTTAATGAAAAATACGCAAAAGCCCACAAACCGCCCCCCCCCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                         <  1:301016/140‑1 (MQ=21)
                                            aaaTCAGCTTCCTGATATTAATATGGTTACGGGCGTAATCACAATTTCGCCAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGTTTTGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                         <  1:671575/140‑1 (MQ=25)
                                             aaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTTATGACAAATTCGAAAATGCGCAAAATCCGTCCCCTCGCCCCTTTGGGGTGAGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAcccc                                                                                        <  2:1003643/140‑1 (MQ=37)
                                                   cTTCCTGTTATTAATACGGTTACGGGCGTAATGACAATTTCTACAAAGCCCACAATCCGTCCCCTCGCCCATTTGGGGTGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:429916/140‑1 (MQ=25)
                                                   cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGACCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  >  1:718304/1‑140 (MQ=255)
                                                    ttCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                 >  1:1297380/1‑140 (MQ=255)
                                                     tCCTGTTATTAATAAGTTTATGGGCGTATTCACAAATTCGCCAAAGCGCACAACCCGTCCCCGCGCCCATTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  2:394175/140‑1 (MQ=25)
                                                     tCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  2:394725/140‑1 (MQ=255)
                                                      ccTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAGTTACCCTCAGCCCAGCCCTCAc                                                                               >  2:1434375/1‑140 (MQ=255)
                                                        tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAccc                                                                             >  1:236676/1‑140 (MQ=255)
                                                                   cggTTATGGGCTTAATCAAAATTTCCCAAACGCGCACAACCGGACCCCTCGCCCCTTTGGGTTGAGGTTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                  <  1:142220/139‑1 (MQ=255)
                                                                    ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  1:135420/1‑140 (MQ=255)
                                                                    ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  2:312436/1‑140 (MQ=255)
                                                                    ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  2:69182/1‑140 (MQ=255)
                                                                       tAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgggg                                                              >  1:1775016/1‑140 (MQ=255)
                                                                       tAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgggg                                                              >  1:327069/1‑140 (MQ=255)
                                                                             cttAATGACAAATTCGAAAACGCCAACAATCCGTCCCCACGCCCCTTGGGGTTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  1:1308315/138‑1 (MQ=255)
                                                                             cttAATGAAAAATTCTAAAAAGCGCACAATCCGGCCCCCCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  1:1042292/138‑1 (MQ=255)
                                                                             cttAATCACAAATCCGACAACGCGCACAACCCGTCCCCCCGCCCCTTTGGGGAGCGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  1:1323244/138‑1 (MQ=255)
                                                                             cGTAATGACAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  2:246662/140‑1 (MQ=255)
                                                                             cGTAATGAAAAATTCGAAAATGCGAACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  1:906387/140‑1 (MQ=25)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGTTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:312356/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  2:226050/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  2:1403024/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  2:1611997/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:1068440/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:1068131/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:731692/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:1814300/1‑140 (MQ=255)
                                                                               tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:403712/1‑140 (MQ=255)
                                                                                  ttaCAAATTCTAAAACGCGCACAATCCGTCCCCTCGCCCGTTTGGGGTGTGGTTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1189399/138‑1 (MQ=255)
                                                                                  tGACAAATTCGACACAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGTTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1075456/140‑1 (MQ=25)
                                                                                  tGACAAATTCCACAAAGCGCACAATCGGTCCCCTCGCCCGTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1256302/140‑1 (MQ=255)
                                                                                  gtaaaaaTTCGTCTAAGCGCACAATCCGGCGCCTCGCCCCTTTGGGTAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTAACCCCAGCCCTCACCCTGGAAGGGTGAGGGGGCAGATCGgcg                                                   <  2:1668770/136‑1 (MQ=21)
                                                                                      aatttCGACAACGCGCACAATCCGCCCCCTCGCCCTTTGGGGGTGTGGTTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                               <  1:359226/137‑1 (MQ=255)
                                                                                      aatttCGAAAAAGGGCAAAATCCGCCCCCTCGCCCCTTTGGGGAGTGTGTTTGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                               <  1:158172/137‑1 (MQ=25)
                                                                                      aaaTTCTAAAAAGCGCACAATCCTTCCCCTCGCCCTTTTGGGTTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                               <  1:158198/140‑1 (MQ=255)
                                                                                      aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                               >  1:258899/1‑140 (MQ=255)
                                                                                      aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                               >  2:1705898/1‑140 (MQ=255)
                                                                                      aaaTTCGAAAAAGCGAAAAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                               <  1:952075/140‑1 (MQ=25)
                                                                                         ttCGACAACGCGCACAATCCGACCCCTCGCCCCTTTGGGGAGTGGTTTCGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  1:1371829/140‑1 (MQ=25)
                                                                                         ttCGAAAAAGCGCACAATCCGTCCCCTCCCCCCTTTGGGGAGCGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  1:1651820/140‑1 (MQ=255)
                                                                                          tctaCAAAGCGCACAATCCGGCCCCCCGCCCCTTTGGGGTGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                           <  1:1062432/137‑1 (MQ=255)
                                                                                              cAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  1:449669/140‑1 (MQ=255)
                                                                                              cAAAGCGCACAATCCGTCCCGTCGCCCCTTGGGGTAGAGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  1:1618135/140‑1 (MQ=255)
                                                                                              cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  2:326167/140‑1 (MQ=255)
                                                                                              cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  2:1696680/140‑1 (MQ=255)
                                                                                              cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  2:1756217/140‑1 (MQ=255)
                                                                                              cAAAGCGCAAAACCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                       <  1:451152/140‑1 (MQ=25)
                                                                                                        aaTCCGTCCCCTCGGCCCTTTGGGGAGAGGGTTAGGGTGATGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTAt                             <  2:901329/140‑1 (MQ=255)
                                                                                                           ccGTCCCCTCGCCCCTTTGGGTAGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGGGAACAATAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACTTTGCGCTTTTGCg                          <  2:1280269/140‑1 (MQ=25)
                                                                                                             gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                 >  1:547028/1‑115 (MQ=255)
                                                                                                             gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                 <  2:547028/115‑1 (MQ=255)
                                                                                                             gTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGTCCTCACCCGGGAAGGGGGAGGGGGCAGAACGGCGCAGGACATCACAGTGCGCTTATGCGaa                        >  1:896277/1‑140 (MQ=25)
                                                                                                              tCCCCTCGCCCCTTGGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCGGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAAt                       >  2:488880/1‑140 (MQ=255)
                                                                                                                cgcTCGCCCGTTTGGGGAGCGGGTTCGGGTGAGGGGCACAGCCTGCACTGGTGCGTACATTAACCCTCACCCCAGCCCTCACCCTGGTAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     <  2:619795/138‑1 (MQ=255)
                                                                                                                cccTCGCCCATTTGGGGAGAGGTTTAGGTTGAGGGTAACATCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAAAGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     <  2:619812/140‑1 (MQ=255)
                                                                                                                  cTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                   >  1:1705592/1‑140 (MQ=255)
                                                                                                                  cTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                   >  1:1706183/1‑140 (MQ=255)
                                                                                                                     gCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaa                >  1:689957/1‑140 (MQ=255)
                                                                                                                        ccTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaataa             <  1:1614645/140‑1 (MQ=255)
                                                                                                                           ttGGGGAGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGc          <  1:1530855/140‑1 (MQ=255)
                                                                                                                           ttGGGGAGCGGGTTTGGGGGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGc          <  1:381560/140‑1 (MQ=255)
                                                                                                                           ttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGc          <  2:427224/140‑1 (MQ=255)
                                                                                                                           ttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGc          <  1:381480/140‑1 (MQ=255)
                                                                                                                             ggggTGTGGGTTGGGGTGTGGGGCACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCtt        <  1:74876/140‑1 (MQ=255)
                                                                                                                               ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtct      <  2:584075/140‑1 (MQ=255)
                                                                                                                                   agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGTGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:725017/1‑140 (MQ=255)
                                                                                                                                   agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:913874/1‑140 (MQ=255)
                                                                                                                                   agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:845551/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                     
CCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624467‑624737

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