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.

CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAA  >  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)
                gAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctg                                                                                                                                      >  4:669304/1‑150 (MQ=255)
                 aaGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgt                                                                                                                                     <  2:1135219/150‑1 (MQ=255)
                      gtAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGtt                                                                                                                                >  2:293935/1‑150 (MQ=255)
                           tCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGttattta                                                                                                                           >  2:269767/1‑150 (MQ=255)
                              gCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTc                                                                                                                        >  4:633550/1‑150 (MQ=255)
                                        acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTc                                                                                                                        >  4:624793/1‑140 (MQ=255)
                                        acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTc                                                                                                                        <  3:624793/140‑1 (MQ=255)
                                        acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  1:586603/150‑1 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAGCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  2:854136/147‑1 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              >  1:854136/1‑147 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  2:1005533/147‑1 (MQ=255)
                                           ctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              >  1:1005533/1‑147 (MQ=255)
                                              acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgt                                                                                                                                     <  2:846228/121‑1 (MQ=255)
                                              acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgt                                                                                                                                     >  1:846228/1‑121 (MQ=255)
                                              acCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTc                                                                                                        >  1:1134915/1‑150 (MQ=255)
                                                cTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttt                                                                                                                                            <  2:907878/113‑1 (MQ=255)
                                                cTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttt                                                                                                                                            >  1:907878/1‑113 (MQ=255)
                                                        gTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTg                                                                                                    >  4:440467/1‑144 (MQ=255)
                                                        gTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTg                                                                                                    <  3:440467/144‑1 (MQ=255)
                                                                  gtTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    <  3:28173/150‑1 (MQ=255)
                                                                     gACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              <  1:518878/121‑1 (MQ=255)
                                                                     gACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAttttt                                                                                                              >  2:518878/1‑121 (MQ=255)
                                                                       ccAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATc                                                                               <  3:669304/150‑1 (MQ=255)
                                                                            atTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCaa                                                                          <  2:30596/150‑1 (MQ=255)
                                                                             ttGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGaa                                                                                            <  3:522412/131‑1 (MQ=255)
                                                                             ttGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGaa                                                                                            >  4:522412/1‑131 (MQ=255)
                                                                                aTAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                                          >  2:773486/1‑83 (MQ=255)
                                                                                aTAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                                          <  1:773486/83‑1 (MQ=255)
                                                                                 tAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATgaga                                                                     <  1:269767/150‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgt                                                                                                                                     >  1:342116/1‑85 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgt                                                                                                                                     <  2:342116/85‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTc                                                                                                             >  2:277914/1‑109 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTc                                                                                                             <  1:277914/109‑1 (MQ=255)
                                                                                  aaCAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGagag                                                                    >  2:546573/1‑150 (MQ=255)
                                                                                   actaaCCTGCGCTGATTTGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCTCCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                                   <  2:272188/147‑1 (MQ=255)
                                                                                    cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttt                                                                                                                                            >  1:895254/1‑77 (MQ=255)
                                                                                    cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCtttt                                                                                                                                            <  2:895254/77‑1 (MQ=255)
                                                                                    cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAggg                                                                  <  1:450989/150‑1 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGc                                                                                                                                                   <  2:674603/69‑1 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGc                                                                                                                                                   >  1:674603/1‑69 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCt                                                                                                            >  2:981948/1‑107 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCt                                                                                                            <  1:981948/107‑1 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGc                                                                 >  2:1082098/1‑150 (MQ=255)
                                                                                     aaaCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGc                                                                 <  3:633550/150‑1 (MQ=255)
                                                                                       aCCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGa                                                               >  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)
                                                                                           gcgcTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATaaa                                                           >  1:1022203/1‑150 (MQ=255)
                                                                                               tGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAAtctc                                                       <  1:293935/150‑1 (MQ=255)
                                                                                                     gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCgg                                                                             <  4:58365/122‑1 (MQ=255)
                                                                                                     gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCgg                                                                             >  3:58365/1‑122 (MQ=255)
                                                                                                       aaCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGAt                                               >  1:620211/1‑150 (MQ=255)
                                                                                                           aataataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTCCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGACAGGGCGATAAATCCCCACCCGCTCGcc                                           >  2:857340/1‑150 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttt                                                                                                                                           <  4:689455/52‑1 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttt                                                                                                                                           >  3:689455/1‑52 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                                 <  1:820344/109‑1 (MQ=255)
                                                                                                              aataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                                 >  2:820344/1‑109 (MQ=255)
                                                                                                                     ttGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                                   >  2:1050487/1‑116 (MQ=255)
                                                                                                                     ttGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                                   <  1:1050487/116‑1 (MQ=255)
                                                                                                                         tctcAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGttat                                                                                                                              <  2:679668/53‑1 (MQ=255)
                                                                                                                         tctcAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGttat                                                                                                                              >  1:679668/1‑53 (MQ=255)
                                                                                                                           tcAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTTCTCGTTGCGTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCg                           >  2:928341/1‑150 (MQ=255)
                                                                                                                            cAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGc                          <  3:92505/150‑1 (MQ=255)
                                                                                                                             aaaGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGttat                                                                                                                              >  1:256862/1‑49 (MQ=255)
                                                                                                                             aaaGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGttat                                                                                                                              <  2:256862/49‑1 (MQ=255)
                                                                                                                                   gTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGttt                   >  4:46174/1‑150 (MQ=255)
                                                                                                                                    ttAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGtttt                  >  3:261039/1‑150 (MQ=255)
                                                                                                                                    ttAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGtttt                  >  4:366704/1‑150 (MQ=255)
                                                                                                                                    agaaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGtttt                  <  3:333167/148‑1 (MQ=255)
                                                                                                                                      aaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    <  3:660365/82‑1 (MQ=255)
                                                                                                                                      aaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    >  4:660365/1‑82 (MQ=255)
                                                                                                                                        aaaaGGTGCTCAATGAGCACCTTTTTTCGGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTGGCCTGCGAAGCTCTTTCGGATCGGCAAGGAGAGGGCGATAAAGCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCg              <  1:857340/150‑1 (MQ=255)
                                                                                                                                         aaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCgg             >  1:460633/1‑150 (MQ=255)
                                                                                                                                         aaaGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTACTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCgg             >  3:622396/1‑150 (MQ=255)
                                                                                                                                         aaaGGTGATCATTGAGCACCTTTTTTAAAGTCTGTAATTAATTCGCTGATGTTTCTGCTCGTTGCCTGCGAATCTCTTTCGGATCGGCAATGAAGGGGCGATAAATCTCCACCCGATCGCCATCATACACACTATCGCTTAGTTTTGCCg              >  4:644534/1‑150 (MQ=255)
                                                                                                                                          aaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGa            >  3:282224/1‑150 (MQ=255)
                                                                                                                                          aaGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    >  4:231758/1‑78 (MQ=255)
                                                                                                                                          aaGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                                    <  3:231758/78‑1 (MQ=255)
                                                                                                                                           aGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCttt                                                                                     >  2:904137/1‑76 (MQ=255)
                                                                                                                                           aGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTGTTCTGCTCGTTGCCTGCGAAGCTCttt                                                                                     <  1:904137/76‑1 (MQ=255)
                                                                                                                                             gTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACgg         >  2:71329/1‑150 (MQ=255)
                                                                                                                                             gTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACgg         <  4:282224/150‑1 (MQ=255)
                                                                                                                                               gCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCt       <  3:79475/150‑1 (MQ=255)
                                                                                                                                                cTCATTGAGCACCTTTTTTAACGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCt       <  2:1212250/150‑1 (MQ=255)
                                                                                                                                                 tCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTTCTCTTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGt     >  2:1253612/1‑150 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTTCTTTCTTTTATTTATTCGCTGATTCTTCTGCTCGTTGCCTGCGAAGCTCGTTCGGATCGGCAATGAGAGGCCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTaa   >  2:365927/1‑150 (MQ=255)
                                                                                                                                                   aaTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTaa   >  3:709665/1‑150 (MQ=255)
                                                                                                                                                    aTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTaaa  <  2:1022203/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAA  >  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: