Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,355,0490AC68.8% 77.8 / 79.0 77L141W (TTG→TGG) sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (15/9);  new base C (29/24);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTT  >  NC_000913/1354813‑1355280
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
tCGCCGCCGCTTTTACATAGTTCTGTTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCTGTTGTTGGCGCAACCGACAGGTTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGcaa                                                                                                                                                                                                                                        <  1:1043282/238‑1 (MQ=255)
 cgccgccgcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGccacc                                                                                                                                                                                                                                      <  1:691224/239‑1 (MQ=255)
        gcTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAg                                                                                                                                                                                                                               <  1:164621/239‑1 (MQ=255)
                         gTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACa                                                                                                                                                                                                              >  1:717744/1‑239 (MQ=255)
                              acacTTCGATAGTGCTGATACGCATCAGGCGGATCCCTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGccac                                                                                                                                                                                                                                       <  1:1189901/209‑1 (MQ=255)
                              acacTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGccac                                                                                                                                                                                                                                       >  2:1189901/1‑209 (MQ=255)
                              acacTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTc                                                                                                                                                                                                         >  2:242971/1‑239 (MQ=255)
                                    cGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGTAAACAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGa                                                                                                                                                                                                   >  1:1119360/1‑239 (MQ=255)
                                        aGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAgcg                                                                                                                                                                                               <  1:1178606/239‑1 (MQ=255)
                                            cTGATACGCATCAGGCGGATCACTTCGGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCgcca                                                                                                                                                                                            >  2:99168/1‑238 (MQ=255)
                                            cTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGAGAACAAAACAGCGTCAAAAGAAGAgcca                                                                                                                                                                                            >  2:1110108/1‑238 (MQ=255)
                                                      tCAGTCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACAcc                                                                                                                                                                                 >  1:828491/1‑239 (MQ=255)
                                                         ggCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGCTTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCCGCCACAACACCcgg                                                                                                                                                                              >  1:644629/1‑239 (MQ=255)
                                                           cgtaTCCCTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAAAATAATATGGCGGATTGCGCTAATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGa                                                                                                                                                                             <  1:30626/235‑1 (MQ=255)
                                                               tCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAgcg                                                                                                                                                                                               <  1:778927/216‑1 (MQ=255)
                                                               tCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAgcg                                                                                                                                                                                               >  2:778927/1‑216 (MQ=255)
                                                               tCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCCGAACACCGGGATTGaa                                                                                                                                                                        >  1:494614/1‑239 (MQ=255)
                                                                 aCTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGaaaa                                                                                                                                                                      <  1:242971/239‑1 (MQ=255)
                                                                  cTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAg                                                                                                                                                                     <  2:1130295/239‑1 (MQ=255)
                                                                        ttgttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCgtga                                                                                                                                                                                                                          >  1:547117/1‑180 (MQ=255)
                                                                        ttgttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCgtga                                                                                                                                                                                                                          <  2:547117/180‑1 (MQ=255)
                                                                           ttgGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGAAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAg                                                                                                                                                            >  1:56009/1‑239 (MQ=255)
                                                                             ggCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca                                                                                                                                                                                   <  2:589463/214‑1 (MQ=255)
                                                                             ggCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAaca                                                                                                                                                                                   >  1:589463/1‑214 (MQ=255)
                                                                             ggCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGc                                                                                                                                                          >  1:204548/1‑239 (MQ=255)
                                                                             ggCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGc                                                                                                                                                          >  1:140288/1‑239 (MQ=255)
                                                                                   aCCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGc                                                                                                                                                    >  2:722015/1‑239 (MQ=255)
                                                                                          gggTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGc                                                                                                                                                                                                            <  2:255703/176‑1 (MQ=255)
                                                                                          gggTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGc                                                                                                                                                                                                            >  1:255703/1‑176 (MQ=255)
                                                                                                  aCGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTAt                                                                                                                                      >  1:409334/1‑238 (MQ=255)
                                                                                                                  ggCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTa                                                                                                                                                              <  2:227836/198‑1 (MQ=255)
                                                                                                                  ggCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTa                                                                                                                                                              >  1:227836/1‑198 (MQ=255)
                                                                                                                    cGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCg                                                                                                                                                         >  2:833845/1‑201 (MQ=255)
                                                                                                                    cGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCg                                                                                                                                                         <  1:833845/201‑1 (MQ=255)
                                                                                                                    cGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCt                                                                                                                                   <  2:51972/223‑1 (MQ=255)
                                                                                                                    cGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCt                                                                                                                                   >  1:51972/1‑223 (MQ=255)
                                                                                                                    cGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGt                                                                                                                   >  2:57985/1‑239 (MQ=255)
                                                                                                                    cGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGt                                                                                                                   >  2:721184/1‑239 (MQ=255)
                                                                                                                      gATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCt                                                                                                                                   >  1:954351/1‑221 (MQ=255)
                                                                                                                      gATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCt                                                                                                                                   <  2:954351/221‑1 (MQ=255)
                                                                                                                      gATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTaa                                                                                                                 >  2:217623/1‑239 (MQ=255)
                                                                                                                            gcTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCgg                                                                                                           <  1:99168/239‑1 (MQ=255)
                                                                                                                                          cATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACg                                                                                             >  1:165059/1‑239 (MQ=255)
                                                                                                                                            tCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACggg                                                                                           >  1:364069/1‑239 (MQ=255)
                                                                                                                                             cccGCCATGGAGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGACCAGCAAGGCGATGGCGTTGATCAAGTTAACCCGCCACTTATGGCGCGAAATCCCGGCAATAATCCCCACgggg                                                                                          >  2:1086856/1‑239 (MQ=255)
                                                                                                                                              ccGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGa                                                                                         <  2:828491/239‑1 (MQ=255)
                                                                                                                                                  cATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGcc                                                                                     <  1:721184/239‑1 (MQ=255)
                                                                                                                                                     ggTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGccac                                                                                                                                                                                                                                       <  1:812043/90‑1 (MQ=255)
                                                                                                                                                     ggTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGccac                                                                                                                                                                                                                                       >  2:812043/1‑90 (MQ=255)
                                                                                                                                                        gAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAAt                                                                               <  2:1163033/239‑1 (MQ=255)
                                                                                                                                                           tCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCa                                                                             >  1:109408/1‑238 (MQ=255)
                                                                                                                                                                         atcaatcaACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAgg                                                              <  1:722015/239‑1 (MQ=255)
                                                                                                                                                                         atcaatcaACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAgg                                                              <  2:409334/239‑1 (MQ=255)
                                                                                                                                                                           caatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCg                                                            <  1:57985/239‑1 (MQ=255)
                                                                                                                                                                            aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTa                                                                                                                                       >  1:860688/1‑163 (MQ=255)
                                                                                                                                                                            aatcaaCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTa                                                                                                                                       <  2:860688/163‑1 (MQ=255)
                                                                                                                                                                                 aCGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGAt                                                       >  2:27073/1‑238 (MQ=255)
                                                                                                                                                                                   gCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCaa                                                                       >  1:1008984/1‑220 (MQ=255)
                                                                                                                                                                                   gCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCaa                                                                       <  2:1008984/220‑1 (MQ=255)
                                                                                                                                                                                   gCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGca                                                    <  1:217623/239‑1 (MQ=255)
                                                                                                                                                                                       aaCCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGccagccag                                                                                                                                                                                       >  2:261541/1‑104 (MQ=255)
                                                                                                                                                                                       aaCCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGccagccag                                                                                                                                                                                       <  1:261541/104‑1 (MQ=255)
                                                                                                                                                                                                 cGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAtcgtc                                       >  1:1108564/1‑238 (MQ=255)
                                                                                                                                                                                                   gTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAt                                           <  2:47945/232‑1 (MQ=255)
                                                                                                                                                                                                   gTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAt                                           >  1:47945/1‑232 (MQ=255)
                                                                                                                                                                                                   gTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAAACAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCAgcgttacc                                    >  2:660216/1‑235 (MQ=255)
                                                                                                                                                                                                    tttCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGcc                                    <  1:383999/238‑1 (MQ=255)
                                                                                                                                                                                                          ttCGTAGGGCAGATCGAAACGCCCTGAACCGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                             <  1:844374/239‑1 (MQ=255)
                                                                                                                                                                                                          ttCGTAGAGCAGATCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGa                             <  2:717744/239‑1 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAacac                                                                                                                                                                                  <  2:650595/78‑1 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAacac                                                                                                                                                                                  >  1:650595/1‑78 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAat                                                                                                       <  2:565403/153‑1 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAat                                                                                                       >  1:565403/1‑153 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTa                    >  1:472390/1‑236 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTa                    >  1:138931/1‑236 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTa                    >  2:850970/1‑236 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTa                    >  1:761362/1‑236 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTa                    >  2:985426/1‑236 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGccagccag                                                                                                                                                                                       <  2:641418/73‑1 (MQ=255)
                                                                                                                                                                                                                      aTCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGccagccag                                                                                                                                                                                       >  1:641418/1‑73 (MQ=255)
                                                                                                                                                                                                                                     aaCGGGCCACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGtt  <  2:1108564/239‑1 (MQ=255)
                                                                                                                                                                                                                                         ggCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAgcg                                                                                                                                                                                               >  1:524065/1‑46 (MQ=255)
                                                                                                                                                                                                                                         ggCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAgcg                                                                                                                                                                                               <  2:524065/46‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                       
TCGCCGCCGCTTTTACATAGTTCTGGTCATACACTTCGATAGTGCTGATACGCATCAGGCGGATCACTTCAGTTGTTGGCGCAACCGACAGGGTGATCACGGGCAATATCATATGGCGGATTGCGCTCATGATCATTTCATCCCGCCATGGTGAGTCCGAAAGCCAGGCATCAATCAACGCAAAACCGGTAATCGGTTTCACTTCGTAGAGCAGATCGAAACGCCCTGAAACGGGCAACCAACCTAGCGTGAGTGAACAAAACAGCGTCAACAGAAGCGCCAGCCAGAACACCGGGATTGAAAAGCCCAGTAAGGCGATGGCGTTGATCAAGTTATCCTGCCACTTATGGCGCGTAATCCCGGCAATCATCCCCACGGGGATGCCAACAATCAGCGCAAAACCGAAGGCGAGGATGCACAGTTCCATCGTCGCCGGGAAGACCTCTTTTAACTGTTCGGCAATCGGTT  >  NC_000913/1354813‑1355280

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: