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

GCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGC  >  NC_000913/1133707‑1133992
                                                                                                                                                        |                                                                                                                                       
gCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACa                                                                                                                                         <  1:1156995/151‑1 (MQ=255)
gCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACa                                                                                                                                         <  1:808589/151‑1 (MQ=255)
  gTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAgg                                                                                                                                       >  1:818906/1‑151 (MQ=255)
   tAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGc                                                                                                                                      >  2:453630/1‑151 (MQ=255)
   tAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGc                                                                                                                                      >  2:1057764/1‑151 (MQ=255)
      cTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGt                                                                                                                                   >  1:84496/1‑151 (MQ=255)
       tATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGt                                                                                                                                   >  2:1363347/1‑150 (MQ=255)
                     ccAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCt                                                                                                                    >  2:1006155/1‑151 (MQ=255)
                       aGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTgg                                                                                                                  >  2:897701/1‑151 (MQ=255)
                               aaaaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGcc                                                                                                          >  1:60480/1‑151 (MQ=255)
                                 aaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                        >  1:509798/1‑151 (MQ=255)
                                 aaCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                        >  2:1262430/1‑151 (MQ=255)
                                  aCCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTc                                                                                                        >  1:536581/1‑150 (MQ=255)
                                   ccccGGACCAGATCCTAAACTCGTTGGTAAACTTCCCCTAATGGCTGAGGGGAAAGCCCATGGAA‑TCACGAATTTTATCCGGCGAGGGGGGCGGCCGGCCAGACTCTGAACA‑AACAGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCaa                                                                                                      <  2:106538/151‑1 (MQ=255)
                                             gATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGttt                                                                                            >  2:561651/1‑151 (MQ=255)
                                             gATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGttt                                                                                            >  1:473471/1‑151 (MQ=255)
                                                cccaaaCCCGCTGTTTACCTTCCGCTAACCGCTGAGGGAATAGCTAAATGGA‑CCCCTCATTTTATCCCGGGGGGGGAGCACCCAGCCAGACCCTAAAAC‑AAGGGGCGATAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTc                                                                                          <  1:910345/146‑1 (MQ=255)
                                                 ctaaaCGCGTTGTTTAACTTGCCTTAACCGCGGACGGGAGATCT‑AAGGGATTCACGCATTTTATCCCGGGAGGGGAACACCCACCCAGACACGAAAAC‑AACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcg                                                                                         <  2:774332/147‑1 (MQ=255)
                                                 cccAACCCGGTGTTTACCTTACGATAACCGCGGAGGGGAAAGCTCAAGGGA‑TCACGAATTTTACCCCGGGAGGGGAGCAGCCACCCAGACCCTAAAAA‑AACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgc                                                                                        <  2:1205054/149‑1 (MQ=255)
                                                 cTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgc                                                                                        >  1:1312061/1‑151 (MQ=255)
                                                  tCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcg                                                                                       >  1:1418781/1‑151 (MQ=255)
                                                   cAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgc                                                                                      >  1:1430273/1‑151 (MQ=255)
                                                    aaCACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTcgcgcg                                                                                     >  2:319518/1‑151 (MQ=255)
                                                      cacGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc                                                                                    >  2:231359/1‑150 (MQ=255)
                                                        cGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGt                                                                                 >  2:58233/1‑151 (MQ=255)
                                                         gCTGTTTAACTTACGCTAACGGCGGAGGGGATAGCTCAATGGA‑TCCCGCAATTTATCCCGCGATGGGAGCACCCAGCCAGACACTGAATC‑AACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGtt                                                                                <  1:809330/151‑1 (MQ=255)
                                                           tGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTg                                                                               >  2:864201/1‑150 (MQ=255)
                                                            gtttACCTTACGCAAATCGCTGAGGGGATAGCTAAAGGGA‑TCACGAATTTTATCCCGGGAGGGGAGCGCCCAGCCAGACCCAAAATC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCgcgcgc                                                                                    <  1:517637/142‑1 (MQ=255)
                                                            ggTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGa                                                                              >  2:672859/1‑150 (MQ=255)
                                                            ggTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGa                                                                              >  2:648337/1‑150 (MQ=255)
                                                              taactttaccatAACCGCTGGCGGGATACCTAAAGGGA‑CCCCGAATTTTATCCCGGGGGGGGAGCAGCCAGCCAGACACTAAAAC‑AACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGtt                                                                                <  2:1065518/135‑1 (MQ=255)
                                                                aaCTTCCGCAAACCGCTGGGGGGAAAGCTCAAGGAA‑TCACGAATTTTACCCCGGGGGGGGGCCGGCCAGCCAGACACTGAAAC‑AAGGGGGGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGc                                                                          <  2:1384485/150‑1 (MQ=255)
                                                                aaCTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCt                                                                         >  2:124496/1‑151 (MQ=255)
                                                                   ttACCCAAATCGCGGAGGGAATAGCTAAATGAA‑TCACCAATTTTACCCGGGGGGGGGAGCACCCACCCAGACACTAAACC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt                                                                       <  2:1436413/150‑1 (MQ=255)
                                                                    tACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAc                                                                     >  1:62853/1‑151 (MQ=255)
                                                                      cgtaaaTCGGTGAGGGGATACCTAAAGATA‑TCACGAATTTTATCCCGGGAGGGGACCGCCCACCGAGACACTAAATC‑AACGGGCGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAc                                                                     <  1:1355756/145‑1 (MQ=255)
                                                                       gCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                  >  1:689188/1‑151 (MQ=255)
                                                                         tAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcgcg                                                                >  2:675309/1‑151 (MQ=255)
                                                                          aaccgcgggggggATAGCTAAATGGA‑TCACAAATTTTATCCCGCGAGGGGAACACCCAGCCAGACACTAAAAC‑AACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                  <  1:1363312/138‑1 (MQ=255)
                                                                          aaTCGCTGAGGGGATAGCTAAATGGA‑TCACGCAATTTATCCCGCGAGGGGAGCAGCCAGCCAGACACTAAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTAcgcg                                                                  <  2:36965/148‑1 (MQ=255)
                                                                          aaTCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCgg                                                               >  2:1375752/1‑151 (MQ=255)
                                                                             cGCTGGGGGGTTGGCTCAATGGA‑TCACGCATTTTACCCCGCGAGGGGAGCAGCAACCCAAACACTAAACC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCttt                                                                       <  2:1224298/140‑1 (MQ=255)
                                                                                tGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGaa                                                                             >  1:413326/1‑131 (MQ=255)
                                                                                  aCGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTg                                                       >  2:113174/1‑151 (MQ=255)
                                                                                    gggATAGCTCAATGAC‑CCACACAATTTATCCCGGGGGGGGGGCACCCACCCAGACACTGAACC‑AACGGGGGTAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCa                                                          <  1:424975/146‑1 (MQ=255)
                                                                                    gggATAGCTAAATGGA‑TCACGAATTTTATCCCGCGATGGGAGCAGCCAGCCAGACACTAAAAC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTg                                                       <  2:79550/149‑1 (MQ=255)
                                                                                             cAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtc                                            >  2:403131/1‑151 (MQ=255)
                                                                                             aaaTGTA‑CCACGAATTTTACACCGCGAGGGGAGCAGCCACCAAGACACAGAACC‑AGCGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtc                                            <  1:114427/150‑1 (MQ=255)
                                                                                              aaTGAA‑TCCCCCATTTTATCCCGGGAGGGGACCACCCACCCACACCCTGAATC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtc                                            <  2:1135415/150‑1 (MQ=255)
                                                                                                tGGA‑TCACGCACTTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtctc                                          >  2:866491/1‑150 (MQ=255)
                                                                                                gggA‑CCACGAAATTTATCCCGCGAGGGGGGCAGCCACCCAGACACTAAAAC‑AGCGGGCGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTtctc                                          <  1:374257/149‑1 (MQ=255)
                                                                                                 ggA‑CCACGAATTTTATCCGGCGAGGGGAGCAGCCAGCCAGCCACTGAATC‑AAGGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTctct                                         <  2:954215/150‑1 (MQ=255)
                                                                                                  gA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGc                                       >  1:703795/1‑151 (MQ=255)
                                                                                                  gA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGc                                       >  2:529341/1‑151 (MQ=255)
                                                                                                        cgaaaTTTATCCGGCGAGGGGAGCAGCAAGCCAGCCACTGAACC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGccc                                     <  1:791225/145‑1 (MQ=255)
                                                                                                         gCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcgc                                 >  2:1051481/1‑151 (MQ=255)
                                                                                                             tttATCCGGCGGGGGGACCAGCCAGCCAGACACTAAACA‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCAcgc                                 <  1:911656/147‑1 (MQ=255)
                                                                                                             tttATCCCGGGGGGGGAGGACCCAGCAAGCAACTAAACA‑AACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACgcgc                               <  1:251620/149‑1 (MQ=255)
                                                                                                              ttATCCCGCGAGGGGAGCAGCCAGCCAGACCCGGAATC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCAc                             <  2:214409/150‑1 (MQ=255)
                                                                                                              ttATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACg                            >  1:1307886/1‑151 (MQ=255)
                                                                                                               tataCCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGt                           >  2:1099015/1‑151 (MQ=255)
                                                                                                                  ccggggagggggGCACCCAGCCAGACA‑TGAAACAAACGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                          <  2:1111972/141‑1 (MQ=255)
                                                                                                                  aCCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                        >  2:1454268/1‑151 (MQ=255)
                                                                                                                   cggggaggggAGCCGCCACCAAGACACTAAACC‑AGCGGGGGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGt                       <  1:1049549/144‑1 (MQ=255)
                                                                                                                   ccGCGAGGGGAGCAGCCAGCCAGACACTGAATC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGtt                          <  1:176703/148‑1 (MQ=255)
                                                                                                                   ccGCGAGGGGAGCACCAAGCCAAATACAAAAAC‑AACGGGCGAAACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTgg                        <  2:192013/150‑1 (MQ=255)
                                                                                                                          gggAGCACCCAGCAAGACACAAAATC‑AACGGGGGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTc                      <  2:1430313/145‑1 (MQ=255)
                                                                                                                             agcagcCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTc              >  2:522810/1‑150 (MQ=255)
                                                                                                                                  ccagccagACCCGAAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCa             <  1:561635/146‑1 (MQ=255)
                                                                                                                                  ccagccagACACTGAATA‑AACGGGCGGAAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAc           <  2:1270855/148‑1 (MQ=255)
                                                                                                                                  ccagccagACACTAAACC‑AACGGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCa             <  2:619568/146‑1 (MQ=255)
                                                                                                                                  ccacccacACACGAAACC‑AACGGGCGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAAc           <  1:1215558/148‑1 (MQ=255)
                                                                                                                                     gccagACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCGGTGGAAAGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCgg      >  2:156985/1‑150 (MQ=255)
                                                                                                                                      ccagCCACTAAACC‑AACGGGCGGTACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCgg      <  2:1363378/149‑1 (MQ=255)
                                                                                                                                      ccagACACTAAATC‑AAGGGGCGGTACCCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCggg     <  1:897678/150‑1 (MQ=255)
                                                                                                                                            acTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgcgc  >  1:1166539/1‑147 (MQ=255)
                                                                                                                                             ctaaacc‑AACGGGCGTTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGgc    <  2:820729/138‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                       
GCGTAACTATCAGTCTAACGCCCAGACCATCAAAACCCAGGACCAGATCCTCAACACGCTGGTTAACTTACGCTAATCGCTGACGGGATAGCTCAATGGA‑TCACGCAATTTATACCGCGATGGGAGCAGCCAGCCAGACACTGAATC‑AACAGGCGGTAACCGCCAGTAATCTGGCCAATGCCTCAACGCCCGGTTTTCGCGCGCAGTTGAATGCTTTACGCGCGGTGCCAGTGGAAGGGCTTTCTCTGCCCACGCGCACGTTGGTCACGGCGTCAACGCCGGGCGC  >  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: ‑