Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,5980AT34.4% 92.3 / 14.5 65intergenic (‑88/‑287)fepB/entCferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base T (0/22);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 1.16e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624459‑624737
                                                                                                                                           |                                                                                                                                           
gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                            <  2:2101751/141‑1 (MQ=255)
  aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                          <  2:1508141/141‑1 (MQ=255)
   aaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                         >  2:2270057/1‑141 (MQ=255)
         cTGTTAATCGAAGGGCGTTGCGGTCGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                                                                                                                                   >  2:2349583/1‑141 (MQ=255)
                tAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                            >  2:1213373/1‑141 (MQ=255)
                                   cgCGGGGCGATTCACAAAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGa                                                                                                         <  2:2252286/140‑1 (MQ=255)
                                   agCGGGGCTAGTCACACAAATCAGATTCCTGTTATTAATATGTTTAAGGGCGTAATGACAAATACCACAAAGCGCAAAAGCCGGCCCCTCGCCCCTTTGGGGAGTGTTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGa                                                                                                         <  2:259722/141‑1 (MQ=255)
                                   agCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGa                                                                                                         <  2:599486/141‑1 (MQ=255)
                                   agCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGa                                                                                                         <  2:599391/141‑1 (MQ=255)
                                   agCGGGGCGAGTCTAACAAATCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAATGAAAAATTCCACAAAGCGCACAATCCGTCCCCCCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGa                                                                                                         <  2:1147560/141‑1 (MQ=38)
                                     gggggCGAGTCTCACAAATCAGCTTCGTGTTATTAATACGGTTAAGGGCGTAATGACAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGGAGTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAAc                                                                                                       <  2:1810426/140‑1 (MQ=255)
                                     cGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAAc                                                                                                       <  2:569987/141‑1 (MQ=255)
                                      ggggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACa                                                                                                      >  1:1275569/1‑141 (MQ=255)
                                        ggCGAGTCTCACAAAACAGCTTCCTTTTATTAATAAGGTTAAGGGCTTAATGACAAATTCGACAAAGCGCAAAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                    <  2:1521628/141‑1 (MQ=255)
                                              tatcAAAATTCTGCTTCTTTTTTTTAATTAGTTTATGGGCTTATTCAAAATTTCGAAAACGCGCAAAATCGGTCCCCCCGCCCTTTGGGGTTGTGTGTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:2102319/139‑1 (MQ=255)
                                               ctcACAAATCAGCTACCTGTTATTAATAAGGTTACGGGCGTAATCACAAATTCGACAAAGCGCCAAATCCGTCCCCACGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTc                                                                                             <  2:144303/141‑1 (MQ=37)
                                               cccacaAATCCGCTTCCTGTAATTAATACGGTTAAGGGCGTAATGAAAAATTCGACAACGCGCAAAATCCGTCCCCCCGCCCCTTTGGGTATTGGGTTTGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTc                                                                                             <  1:93492/139‑1 (MQ=255)
                                               caaaaaaaTCTGATTCTTTTTTTTAATATGGTTTGGGGCTTATTTAAAATTTCGAAAAAGCAAAAAATCCGTCCCCTCGCCGTTTGGGGGTTTGGGTTTGGTTGGGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTc                                                                                             <  1:1756577/136‑1 (MQ=255)
                                               atcACAAATCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTc                                                                                             <  2:1998583/140‑1 (MQ=255)
                                                 cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAc                                                                                           >  1:625805/1‑141 (MQ=255)
                                                 cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAc                                                                                           >  1:624804/1‑141 (MQ=255)
                                                    aaaTCAGCTTCCTTTTATTAATATGGTTAAGGGCTTATTCAAAAATTCCAAAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGATTGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAcccc                                                                                        <  1:1672721/141‑1 (MQ=25)
                                                    aaaTCAGCTTCCTGTTATTATTAAGGTTACGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGGTGAGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAcccc                                                                                        <  2:2664549/141‑1 (MQ=255)
                                                    aaaTCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAATTACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAcccc                                                                                        <  2:2663875/141‑1 (MQ=255)
                                                         aGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                   >  1:1094041/1‑141 (MQ=255)
                                                         aGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                   >  1:1094047/1‑141 (MQ=255)
                                                         aGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                   >  1:1094307/1‑141 (MQ=255)
                                                         aGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                   >  1:1743851/1‑141 (MQ=255)
                                                           ctaccTGTTATTAATATGGTTACGGGCTTAAACACAAATACCACAAAGCGAACAAACCGACCCCTCGCCCCTTTGGGGAGAGGTTTAGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                 <  2:1738044/138‑1 (MQ=21)
                                                           cTTCCTGTTATTAATAAGTTTAAGGGCGTAATCACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCTAACCCTc                                                                                 <  2:1723876/141‑1 (MQ=37)
                                                           cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACACTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                 <  2:491828/141‑1 (MQ=255)
                                                           cTTCCTGTAATTAATAAGTTTACGGCCGAAAGGACAAATTCCAAAAACCCCACAAGCCGTCCCCACGCCCCTTGGGGTTTTGGGTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                 <  1:480576/141‑1 (MQ=255)
                                                             tCCTTTTATTAATAAGTTTAAGGGCTTAATGAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAc                                                                               <  2:1099890/141‑1 (MQ=37)
                                                             tCCTTTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAc                                                                               <  2:430111/141‑1 (MQ=255)
                                                                    attaATTCGTTTATGGGCTTAATCAAAAATCCCAAAAAGCGCACATTCGGTCCCCGCCCCCCTTTGGGTTGTGGTTTTGGGGTGGGGGAACCGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        <  1:1089267/141‑1 (MQ=255)
                                                                    attaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        >  2:2703492/1‑141 (MQ=255)
                                                                    attaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        >  2:2500034/1‑141 (MQ=255)
                                                                    attaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        >  2:2703503/1‑141 (MQ=255)
                                                                     ttaATAAGGTTAAGGGCTTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAg                                                                       <  2:761756/141‑1 (MQ=255)
                                                                          atggTTAAGGGCTTAATTACAAATCCTACAACGCGAACAACCCGACCCCCCGCCCCTTGGGGTTGTGGTTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCTCCCTGGAAGGgaga                                                                  <  1:1637591/139‑1 (MQ=19)
                                                                           agtttACGGGCTTAATCAAAAATTCTAAAAAGCGCACAATCGGTCCCCTCGCCCGTTTGGGTTGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 <  1:1311332/138‑1 (MQ=255)
                                                                           agtttAAGGGCTTAAGTAAAATTTCCAAAAAGCGCACAATCCGCCCCCTCCCCCTTTGGGGGTGTGGTTTTGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 <  1:2270057/138‑1 (MQ=255)
                                                                           aGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 >  2:2101760/1‑141 (MQ=255)
                                                                           aGGGTAAGGGCGTAATGACCAATACGACAAAGCGCACAAACCGTCCCCTCGCCCCCTTGGGGAGAGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCCCACCCCAGCCCTCACCCTGGAAGGGagag                                                                 <  2:824413/141‑1 (MQ=25)
                                                                            ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgg                                                                >  2:2702331/1‑141 (MQ=255)
                                                                                      gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:1158244/1‑141 (MQ=255)
                                                                                      gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:2331595/1‑141 (MQ=255)
                                                                                          tGACAAATTCGACAACGCGCACAATCCGTCCCCACCCCCCTTTGGGGAGAGGGTTTGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                  <  1:1910031/141‑1 (MQ=255)
                                                                                          tGAAAATTTCTACAACGCCCACAATCCGTCCCCCCCCCCATTGGGGTCGCGGTTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                  <  1:840813/141‑1 (MQ=255)
                                                                                          tGAAAAATTCGACAACGCGCAAAATCCGTCCCCTCGCCCCTTTGGGGAGTGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                  <  1:1959953/141‑1 (MQ=25)
                                                                                          tGAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTCGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                  <  1:1959925/141‑1 (MQ=255)
                                                                                              aatttCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGa                                              <  1:2660403/138‑1 (MQ=25)
                                                                                              aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGa                                              >  2:1565078/1‑141 (MQ=255)
                                                                                                 ttctaaaaaGCGAACAATCCGTCCCCTCGCCCTTTTGGGTTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                           <  1:331407/135‑1 (MQ=255)
                                                                                                 ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACCTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                           <  2:604233/141‑1 (MQ=255)
                                                                                                 ttCGACAAAGCGCACAATCCGGCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                           <  2:2230115/141‑1 (MQ=255)
                                                                                                      cAAAGCACACAATCCGTCCCCTCGCCCCTTGGGGTAGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                      <  1:181616/141‑1 (MQ=255)
                                                                                                      cAAAGCACACAATCCGTCCCCCCGCCCCTTTGGGGTGAGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                      <  1:181561/141‑1 (MQ=255)
                                                                                                                       ccccTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     >  1:808413/1‑141 (MQ=255)
                                                                                                                       ccccTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     >  1:808216/1‑141 (MQ=255)
                                                                                                                        cccTCGCCCCTTTGGGGAGTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCa                    <  2:2178350/141‑1 (MQ=255)
                                                                                                                         ccTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                   >  1:486718/1‑141 (MQ=255)
                                                                                                                                  tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCGTATGCGAATCCATCAATAATGc          >  1:2412926/1‑141 (MQ=255)
                                                                                                                                     ggggAGTGGTTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtc       <  1:1403695/141‑1 (MQ=255)
                                                                                                                                     ggggAGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGATTCCATCAATAATGCTtc       <  1:1277681/141‑1 (MQ=255)
                                                                                                                                     ggggAGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtc       <  1:1403313/141‑1 (MQ=255)
                                                                                                                                          gagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGACCGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAAAAATGCTTCTCAtt  >  1:451862/1‑141 (MQ=255)
                                                                                                                                          gagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:867578/1‑141 (MQ=255)
                                                                                                                                          gagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  2:30868/1‑141 (MQ=255)
                                                                                                                                          gagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:1917611/1‑141 (MQ=255)
                                                                                                                                          gagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  2:1959689/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624459‑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: