Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi2,753,9860AT23.7% 205.5 / 71.0 97intergenic (+40/+22)bamE/ratBlipoprotein component of BamABCDE OM biogenesis complex/UPF0125 family protein
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (11/12);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAA  >  NC_000913_3_bme_pgi/2753838‑2754135
                                                                                                                                                    |                                                                                                                                                      
ccAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                       <  1:597129/150‑1 (MQ=255)
    cAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGc                                                                                                                                                   <  1:880981/150‑1 (MQ=255)
      aCCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCAc                                                                                                                                                 <  2:1075198/150‑1 (MQ=255)
         aGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtt                                                                                                                                              <  2:137702/150‑1 (MQ=255)
          ggTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttt                                                                                                                                             <  1:327036/150‑1 (MQ=255)
          ggTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCttt                                                                                                                                             <  3:442261/150‑1 (MQ=255)
                gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttt                                                                                                                                             <  1:342180/144‑1 (MQ=255)
                gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttt                                                                                                                                             >  2:342180/1‑144 (MQ=255)
                gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                                          <  2:783231/147‑1 (MQ=255)
                gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                                          >  1:783231/1‑147 (MQ=255)
                gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTt‑ctg                                                                                                                                      >  4:669304/1‑150 (MQ=255)
                 aaGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTt‑ctgt                                                                                                                                     <  2:1135219/150‑1 (MQ=255)
                      gtAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGtt                                                                                                                                >  2:293935/1‑150 (MQ=255)
                           tCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGttattta                                                                                                                           >  2:269767/1‑150 (MQ=255)
                              gCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTc                                                                                                                        >  4:633550/1‑150 (MQ=255)
                                        acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTc                                                                                                                        >  4:624793/1‑140 (MQ=255)
                                        acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTc                                                                                                                        <  3:624793/140‑1 (MQ=255)
                                        acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  1:586603/150‑1 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAGCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  2:854136/147‑1 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              >  1:854136/1‑147 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  2:1005533/147‑1 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              >  1:1005533/1‑147 (MQ=255)
                                              acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTt‑ctgt                                                                                                                                     <  2:846228/121‑1 (MQ=255)
                                              acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTt‑ctgt                                                                                                                                     >  1:846228/1‑121 (MQ=255)
                                              acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTc                                                                                                        >  1:1134915/1‑150 (MQ=255)
                                                cTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttt                                                                                                                                            <  2:907878/113‑1 (MQ=255)
                                                cTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttt                                                                                                                                            >  1:907878/1‑113 (MQ=255)
                                                        gTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTg                                                                                                    >  4:440467/1‑144 (MQ=255)
                                                        gTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTg                                                                                                    <  3:440467/144‑1 (MQ=255)
                                                                  gtTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    <  3:28173/150‑1 (MQ=255)
                                                                     gACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  1:518878/121‑1 (MQ=255)
                                                                     gACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGAttttt                                                                                                              >  2:518878/1‑121 (MQ=255)
                                                                       ccAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATc                                                                               <  3:669304/150‑1 (MQ=255)
                                                                            atTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCaa                                                                          <  2:30596/150‑1 (MQ=255)
                                                                             ttGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGaa                                                                                            <  3:522412/131‑1 (MQ=255)
                                                                             ttGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGaa                                                                                            >  4:522412/1‑131 (MQ=255)
                                                                                aTAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                                          >  2:773486/1‑83 (MQ=255)
                                                                                aTAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                                          <  1:773486/83‑1 (MQ=255)
                                                                                 tAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATgaga                                                                     <  1:269767/150‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTt‑ctgt                                                                                                                                     >  1:342116/1‑85 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTt‑ctgt                                                                                                                                     <  2:342116/85‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTc                                                                                                             >  2:277914/1‑109 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTc                                                                                                             <  1:277914/109‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGagag                                                                    >  2:546573/1‑150 (MQ=255)
                                                                                   actaaCCTGCGCTGATTTGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCTCCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                                   <  2:272188/147‑1 (MQ=255)
                                                                                    cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttt                                                                                                                                            >  1:895254/1‑77 (MQ=255)
                                                                                    cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttt                                                                                                                                            <  2:895254/77‑1 (MQ=255)
                                                                                    cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAggg                                                                  <  1:450989/150‑1 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGc                                                                                                                                                   <  2:674603/69‑1 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGc                                                                                                                                                   >  1:674603/1‑69 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCt                                                                                                            >  2:981948/1‑107 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCt                                                                                                            <  1:981948/107‑1 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGc                                                                 >  2:1082098/1‑150 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGc                                                                 <  3:633550/150‑1 (MQ=255)
                                                                                       aCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGa                                                               >  2:566117/1‑150 (MQ=255)
                                                                                          tGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttttaac                                                                                                                                       <  2:914870/76‑4 (MQ=255)
                                                                                          tGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttttaac                                                                                                                                       >  1:914870/1‑73 (MQ=255)
                                                                                          tGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                       <  2:621641/60‑1 (MQ=255)
                                                                                          tGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                                       >  1:621641/1‑60 (MQ=255)
                                                                                           gcgcTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATaaa                                                           >  1:1022203/1‑150 (MQ=255)
                                                                                               tGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAAtctc                                                       <  1:293935/150‑1 (MQ=255)
                                                                                                     gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCgg                                                                             <  4:58365/122‑1 (MQ=255)
                                                                                                     gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCgg                                                                             >  3:58365/1‑122 (MQ=255)
                                                                                                       aaCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGAt                                               >  1:620211/1‑150 (MQ=255)
                                                                                                           aataataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTCCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGACAGGGCGATAAATCCCCACCCGCTCGcc                                           >  2:857340/1‑150 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttt                                                                                                                                           <  4:689455/52‑1 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttt                                                                                                                                           >  3:689455/1‑52 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                                 <  1:820344/109‑1 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                                 >  2:820344/1‑109 (MQ=255)
                                                                                                                     ttGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                                   >  2:1050487/1‑116 (MQ=255)
                                                                                                                     ttGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                                   <  1:1050487/116‑1 (MQ=255)
                                                                                                                         tctcAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGttat                                                                                                                              <  2:679668/53‑1 (MQ=255)
                                                                                                                         tctcAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGttat                                                                                                                              >  1:679668/1‑53 (MQ=255)
                                                                                                                           tcAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTTCTCGTTGCGTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCg                           >  2:928341/1‑150 (MQ=255)
                                                                                                                            cAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGc                          <  3:92505/150‑1 (MQ=255)
                                                                                                                             aaaGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGttat                                                                                                                              >  1:256862/1‑49 (MQ=255)
                                                                                                                             aaaGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGttat                                                                                                                              <  2:256862/49‑1 (MQ=255)
                                                                                                                                   gTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGttt                   >  4:46174/1‑150 (MQ=255)
                                                                                                                                    ttAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGtttt                  >  3:261039/1‑150 (MQ=255)
                                                                                                                                    ttAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGtttt                  >  4:366704/1‑150 (MQ=255)
                                                                                                                                    agaaaaaaGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGtttt                  <  3:333167/148‑1 (MQ=255)
                                                                                                                                      aaaaaaGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    <  3:660365/82‑1 (MQ=255)
                                                                                                                                      aaaaaaGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    >  4:660365/1‑82 (MQ=255)
                                                                                                                                        aaaaGGTGCTCAATGAGCACCTTTTTT‑CGGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTGGCCTGCGAAGCTCTTTCGGATCGGCAAGGAGAGGGCGATAAAGCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCg              <  1:857340/150‑1 (MQ=255)
                                                                                                                                         aaaGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCgg             >  1:460633/1‑150 (MQ=255)
                                                                                                                                         aaaGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTACTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCgg             >  3:622396/1‑150 (MQ=255)
                                                                                                                                         aaaGGTGATCATTGAGCACCTTTTTTAAAGTCTGTAATTAATTCGCTGATGTTTCTGCTCGTTGCCTGCGAATCTCTTTCGGATCGGCAATGAAGGGGCGATAAATCTCCACCCGATCGCCATCATACACACTATCGCTTAGTTTTGCCg              >  4:644534/1‑150 (MQ=255)
                                                                                                                                          aaGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGa            >  3:282224/1‑150 (MQ=255)
                                                                                                                                          aaGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    >  4:231758/1‑78 (MQ=255)
                                                                                                                                          aaGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    <  3:231758/78‑1 (MQ=255)
                                                                                                                                           aGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCttt                                                                                     >  2:904137/1‑76 (MQ=255)
                                                                                                                                           aGGTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTGTTCTGCTCGTTGCCTGCGAAGCTCttt                                                                                     <  1:904137/76‑1 (MQ=255)
                                                                                                                                             gTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACgg         >  2:71329/1‑150 (MQ=255)
                                                                                                                                             gTGCTCATTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACgg         <  4:282224/150‑1 (MQ=255)
                                                                                                                                               gCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCt       <  3:79475/150‑1 (MQ=255)
                                                                                                                                                cTCATTGAGCACCTTTTTTAACGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCt       <  2:1212250/150‑1 (MQ=255)
                                                                                                                                                 tCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTTCTCTTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGt     >  2:1253612/1‑150 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTT‑CTTTCTTTTATTTATTCGCTGATTCTTCTGCTCGTTGCCTGCGAAGCTCGTTCGGATCGGCAATGAGAGGCCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTaa   >  2:365927/1‑150 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTaa   >  3:709665/1‑150 (MQ=255)
                                                                                                                                                    aTGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTaaa  <  2:1022203/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTT‑CTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAA  >  NC_000913_3_bme_pgi/2753838‑2754135

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