Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,133,8570AG41.3% 106.2 / 22.2 75Q19R (CAG→CGG) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (40/4);  new base G (0/31);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 4.61e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGC  >  NC_000913/1133707‑1133992
                                                                                                                                                        |                                                                                                                                       
gCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACa                                                                                                                                         <  1:1156995/151‑1 (MQ=255)
gCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACa                                                                                                                                         <  1:808589/151‑1 (MQ=255)
  gTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAgg                                                                                                                                       >  1:818906/1‑151 (MQ=255)
   tAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGc                                                                                                                                      >  2:453630/1‑151 (MQ=255)
   tAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGc                                                                                                                                      >  2:1057764/1‑151 (MQ=255)
      cTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGt                                                                                                                                   >  1:84496/1‑151 (MQ=255)
       tATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGt                                                                                                                                   >  2:1363347/1‑150 (MQ=255)
                     ccAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCt                                                                                                                    >  2:1006155/1‑151 (MQ=255)
                       aGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTgg                                                                                                                  >  2:897701/1‑151 (MQ=255)
                               aaaaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGcc                                                                                                          >  1:60480/1‑151 (MQ=255)
                                 aaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                        >  1:509798/1‑151 (MQ=255)
                                 aaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                        >  2:1262430/1‑151 (MQ=255)
                                  aCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                        >  1:536581/1‑150 (MQ=255)
                                   ccccGGACCAGATCCTAAACTCGTTGGTAAACTTCCCCTAATGGCTGAGGGGAAAGCCCATGGAATCACGAATTTTATCCGGCGAGGGGGGCGGCCGGCCAGACTCTGAACAAACAGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCaa                                                                                                      <  2:106538/151‑1 (MQ=255)
                                             gATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGttt                                                                                            >  2:561651/1‑151 (MQ=255)
                                             gATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGttt                                                                                            >  1:473471/1‑151 (MQ=255)
                                                cccaaaCCCGCTGTTTACCTTCCGCTAACCGCTGAGGGAATAGCTAAATGGACCCCTCATTTTATCCCGGGGGGGGAGCACCCAGCCAGACCCTAAAACAAGGGGCGATAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTc                                                                                          <  1:910345/146‑1 (MQ=255)
                                                 ctaaaCGCGTTGTTTAACTTGCCTTAACCGCGGACGGGAGATCTAAGGGATTCACGCATTTTATCCCGGGAGGGGAACACCCACCCAGACACGAAAACAACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcg                                                                                         <  2:774332/147‑1 (MQ=255)
                                                 cccAACCCGGTGTTTACCTTACGATAACCGCGGAGGGGAAAGCTCAAGGGATCACGAATTTTACCCCGGGAGGGGAGCAGCCACCCAGACCCTAAAAAAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgc                                                                                        <  2:1205054/149‑1 (MQ=255)
                                                 cTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgc                                                                                        >  1:1312061/1‑151 (MQ=255)
                                                  tCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcg                                                                                       >  1:1418781/1‑151 (MQ=255)
                                                   cAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgc                                                                                      >  1:1430273/1‑151 (MQ=255)
                                                    aaCACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgcg                                                                                     >  2:319518/1‑151 (MQ=255)
                                                      cacGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc                                                                                    >  2:231359/1‑150 (MQ=255)
                                                        cGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGt                                                                                 >  2:58233/1‑151 (MQ=255)
                                                         gCTGTTTAACTTACGCTAACGGCGGAGGGGATAGCTCAATGGATCCCGCAATTTATCCCGCGATGGGAGCACCCAGCCAGACACTGAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGtt                                                                                <  1:809330/151‑1 (MQ=255)
                                                           tGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTg                                                                               >  2:864201/1‑150 (MQ=255)
                                                            gtttACCTTACGCAAATCGCTGAGGGGATAGCTAAAGGGATCACGAATTTTATCCCGGGAGGGGAGCGCCCAGCCAGACCCAAAATCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc                                                                                    <  1:517637/142‑1 (MQ=255)
                                                            ggTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGa                                                                              >  2:672859/1‑150 (MQ=255)
                                                            ggTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGa                                                                              >  2:648337/1‑150 (MQ=255)
                                                              taactttaccatAACCGCTGGCGGGATACCTAAAGGGACCCCGAATTTTATCCCGGGGGGGGAGCAGCCAGCCAGACACTAAAACAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGtt                                                                                <  2:1065518/135‑1 (MQ=255)
                                                                aaCTTCCGCAAACCGCTGGGGGGAAAGCTCAAGGAATCACGAATTTTACCCCGGGGGGGGGCCGGCCAGCCAGACACTGAAACAAGGGGGGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGc                                                                          <  2:1384485/150‑1 (MQ=255)
                                                                aaCTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCt                                                                         >  2:124496/1‑151 (MQ=255)
                                                                   ttACCCAAATCGCGGAGGGAATAGCTAAATGAATCACCAATTTTACCCGGGGGGGGGAGCACCCACCCAGACACTAAACCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt                                                                       <  2:1436413/150‑1 (MQ=255)
                                                                    tACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAc                                                                     >  1:62853/1‑151 (MQ=255)
                                                                      cgtaaaTCGGTGAGGGGATACCTAAAGATATCACGAATTTTATCCCGGGAGGGGACCGCCCACCGAGACACTAAATCAACGGGCGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAc                                                                     <  1:1355756/145‑1 (MQ=255)
                                                                       gCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                  >  1:689188/1‑151 (MQ=255)
                                                                         tAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgcg                                                                >  2:675309/1‑151 (MQ=255)
                                                                          aaccgcgggggggATAGCTAAATGGATCACAAATTTTATCCCGCGAGGGGAACACCCAGCCAGACACTAAAACAACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                  <  1:1363312/138‑1 (MQ=255)
                                                                          aaTCGCTGAGGGGATAGCTAAATGGATCACGCAATTTATCCCGCGAGGGGAGCAGCCAGCCAGACACTAAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                  <  2:36965/148‑1 (MQ=255)
                                                                          aaTCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCgg                                                               >  2:1375752/1‑151 (MQ=255)
                                                                             cGCTGGGGGGTTGGCTCAATGGATCACGCATTTTACCCCGCGAGGGGAGCAGCAACCCAAACACTAAACCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt                                                                       <  2:1224298/140‑1 (MQ=255)
                                                                                tGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGaa                                                                             >  1:413326/1‑131 (MQ=255)
                                                                                  aCGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTg                                                       >  2:113174/1‑151 (MQ=255)
                                                                                    gggATAGCTCAATGACCCACACAATTTATCCCGGGGGGGGGGCACCCACCCAGACACTGAACCAACGGGGGTAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCa                                                          <  1:424975/146‑1 (MQ=255)
                                                                                    gggATAGCTAAATGGATCACGAATTTTATCCCGCGATGGGAGCAGCCAGCCAGACACTAAAACAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTg                                                       <  2:79550/149‑1 (MQ=255)
                                                                                             cAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtc                                            >  2:403131/1‑151 (MQ=255)
                                                                                             aaaTGTACCACGAATTTTACACCGCGAGGGGAGCAGCCACCAAGACACAGAACCAGCGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtc                                            <  1:114427/150‑1 (MQ=255)
                                                                                              aaTGAATCCCCCATTTTATCCCGGGAGGGGACCACCCACCCACACCCTGAATCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtc                                            <  2:1135415/150‑1 (MQ=255)
                                                                                                tGGATCACGCACTTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtctc                                          >  2:866491/1‑150 (MQ=255)
                                                                                                gggACCACGAAATTTATCCCGCGAGGGGGGCAGCCACCCAGACACTAAAACAGCGGGCGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtctc                                          <  1:374257/149‑1 (MQ=255)
                                                                                                 ggACCACGAATTTTATCCGGCGAGGGGAGCAGCCAGCCAGCCACTGAATCAAGGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTctct                                         <  2:954215/150‑1 (MQ=255)
                                                                                                  gATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGc                                       >  1:703795/1‑151 (MQ=255)
                                                                                                  gATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGc                                       >  2:529341/1‑151 (MQ=255)
                                                                                                        cgaaaTTTATCCGGCGAGGGGAGCAGCAAGCCAGCCACTGAACCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGccc                                     <  1:791225/145‑1 (MQ=255)
                                                                                                         gCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcgc                                 >  2:1051481/1‑151 (MQ=255)
                                                                                                             tttATCCGGCGGGGGGACCAGCCAGCCAGACACTAAACAAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcgc                                 <  1:911656/147‑1 (MQ=255)
                                                                                                             tttATCCCGGGGGGGGAGGACCCAGCAAGCAACTAAACAAACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACgcgc                               <  1:251620/149‑1 (MQ=255)
                                                                                                              ttATCCCGCGAGGGGAGCAGCCAGCCAGACCCGGAATCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCAc                             <  2:214409/150‑1 (MQ=255)
                                                                                                              ttATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACg                            >  1:1307886/1‑151 (MQ=255)
                                                                                                               tataCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                           >  2:1099015/1‑151 (MQ=255)
                                                                                                                  ccggggagggggGCACCCAGCCAGACATGAAACAAACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                          <  2:1111972/141‑1 (MQ=255)
                                                                                                                  aCCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                        >  2:1454268/1‑151 (MQ=255)
                                                                                                                   cggggaggggAGCCGCCACCAAGACACTAAACCAGCGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGt                       <  1:1049549/144‑1 (MQ=255)
                                                                                                                   ccGCGAGGGGAGCAGCCAGCCAGACACTGAATCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                          <  1:176703/148‑1 (MQ=255)
                                                                                                                   ccGCGAGGGGAGCACCAAGCCAAATACAAAAACAACGGGCGAAACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                        <  2:192013/150‑1 (MQ=255)
                                                                                                                          gggAGCACCCAGCAAGACACAAAATCAACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTc                      <  2:1430313/145‑1 (MQ=255)
                                                                                                                             agcagcCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTc              >  2:522810/1‑150 (MQ=255)
                                                                                                                                  ccagccagACCCGAAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCa             <  1:561635/146‑1 (MQ=255)
                                                                                                                                  ccagccagACACTGAATAAACGGGCGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAc           <  2:1270855/148‑1 (MQ=255)
                                                                                                                                  ccagccagACACTAAACCAACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCa             <  2:619568/146‑1 (MQ=255)
                                                                                                                                  ccacccacACACGAAACCAACGGGCGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAc           <  1:1215558/148‑1 (MQ=255)
                                                                                                                                     gccagACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCGGTGGAAAGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCgg      >  2:156985/1‑150 (MQ=255)
                                                                                                                                      ccagCCACTAAACCAACGGGCGGTACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCgg      <  2:1363378/149‑1 (MQ=255)
                                                                                                                                      ccagACACTAAATCAAGGGGCGGTACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCggg     <  1:897678/150‑1 (MQ=255)
                                                                                                                                            acTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcgc  >  1:1166539/1‑147 (MQ=255)
                                                                                                                                             ctaaaccAACGGGCGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgc    <  2:820729/138‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                       
GCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGATCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATCAACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGC  >  NC_000913/1133707‑1133992

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