BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,011,752	T→C	S168P (TCC→CCC)	fliF →	FliF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,011,752	0	T	C	99.2%	357.2 / ‑5.1	117	S168P (TCC→CCC)	fliF	FliF
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (56/60); total (56/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 121 total aligned reads displayed.

TCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAG  >  NC_000913/2011662‑2011840

tcncGTACCATCGAAACTATCNNCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAccc                                                                                        <  1:2302048/93‑1 (MQ=255)

   cGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCt                                                                                     <  2:3424677/93‑1 (MQ=255)

   cGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCt                                                                                     <  1:3170664/93‑1 (MQ=255)

     tACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTc                                                                                   >  1:1106437/1‑93 (MQ=255)

      aCCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCt                                                                                  <  1:406808/93‑1 (MQ=255)

       ccATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCtc                                                                                 >  1:1828666/1‑92 (MQ=255)

         aTCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCa                                                                               >  1:818780/1‑93 (MQ=255)

          tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt                                                                              >  1:831938/1‑93 (MQ=255)

          tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt                                                                              <  2:4310177/93‑1 (MQ=255)

          tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt                                                                              >  2:1963283/1‑93 (MQ=255)

           cGAATCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATc                                                                             <  2:347246/93‑1 (MQ=255)

            gAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCg                                                                            >  1:4300811/1‑93 (MQ=255)

             aaaCTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCgg                                                                           >  1:4207220/1‑93 (MQ=255)

              aaCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGt                                                                          >  2:3230393/1‑93 (MQ=255)

              aaCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAGGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAGCCCCCTTCTGCATCGGt                                                                          <  1:2202901/93‑1 (MQ=255)

               aCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTg                                                                         <  1:744145/93‑1 (MQ=255)

                cTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGa                                                                        <  1:2455307/93‑1 (MQ=255)

                  aTCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACg                                                                      <  1:135750/93‑1 (MQ=255)

                  aTCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACACCCTTCTGCATCGGTGACg                                                                      >  2:2871216/1‑93 (MQ=255)

                   tCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACgg                                                                     >  2:2447307/1‑93 (MQ=255)

                      gCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTaa                                                                  >  1:1592289/1‑93 (MQ=255)

                       cccAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTaaa                                                                 <  2:1955411/93‑1 (MQ=255)

                        ccagTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAAt                                                                <  1:2591108/90‑1 (MQ=255)

                        ccagTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAAt                                                                >  1:3368157/4‑93 (MQ=255)

                        ccagTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAAt                                                                >  1:3021438/4‑93 (MQ=255)

                           gTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTg                                                             <  2:87572/93‑1 (MQ=255)

                            tAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGt                                                            <  1:373281/93‑1 (MQ=255)

                            tAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGt                                                            >  1:3248606/1‑93 (MQ=255)

                              aaaGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTa                                                          <  2:536821/93‑1 (MQ=255)

                              aaaGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTa                                                          >  1:1708546/1‑93 (MQ=255)

                               aaGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTAc                                                         >  1:3157439/1‑93 (MQ=255)

                               aaGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTAc                                                         >  2:4422734/1‑93 (MQ=255)

                                aGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTAcc                                                        <  2:3676272/93‑1 (MQ=255)

                                aGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTAcc                                                        >  1:2162593/1‑93 (MQ=255)

                                  gggCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCg                                                      <  1:1064309/93‑1 (MQ=255)

                                  gggCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCg                                                      >  1:1676216/1‑93 (MQ=255)

                                   ggCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCgg                                                     <  2:2729477/93‑1 (MQ=255)

                                    gcgcgcGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGc                                                    <  2:3085748/93‑1 (MQ=255)

                                      gcgcgTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCcg                                                  <  2:3715254/93‑1 (MQ=255)

                                           tACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCAc                                             <  2:693077/93‑1 (MQ=255)

                                           tACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCAc                                             <  1:2618554/93‑1 (MQ=255)

                                           tACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCAc                                             <  2:4434860/93‑1 (MQ=255)

                                           tACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCAc                                             <  1:689352/93‑1 (MQ=255)

                                            aCATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACt                                            <  1:1173328/93‑1 (MQ=255)

                                              aTCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCg                                          <  2:1442504/93‑1 (MQ=255)

                                               tCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGa                                         >  2:1658635/1‑93 (MQ=255)

                                               tCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGa                                         >  2:3019800/1‑93 (MQ=255)

                                               tCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGa                                         >  2:3535499/1‑93 (MQ=255)

                                                 tGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATg                                       >  1:1526103/1‑93 (MQ=255)

                                                 tGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATg                                       >  2:148759/1‑93 (MQ=255)

                                                  ggCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGTTAAATCTGTTACCCGGCCGCGCACTCGATGa                                      >  2:2507529/1‑93 (MQ=255)

                                                  ggCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGa                                      <  1:4668781/93‑1 (MQ=255)

                                                  ggCAATCCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGa                                      <  2:190821/93‑1 (MQ=255)

                                                   gCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGaa                                     <  1:3225810/93‑1 (MQ=255)

                                                    cAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAg                                    <  2:3364112/93‑1 (MQ=255)

                                                    cAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAg                                    <  2:4231842/93‑1 (MQ=255)

                                                     aaTGCCGAAACCGTCTTTATTCGTCCGTGAACAAAACCCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAgg                                   <  2:4437921/93‑1 (MQ=255)

                                                     aaTGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTGCCCGGCCGCGCACTCGATGAAgg                                   <  1:379124/93‑1 (MQ=255)

                                                     aaTGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAgg                                   >  1:1262834/1‑93 (MQ=255)

                                                     aaTGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAgg                                   >  2:1482099/1‑93 (MQ=255)

                                                       tGCCGAAGCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCTGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGc                                 <  1:2513079/93‑1 (MQ=255)

                                                          cGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCaaa                              <  1:3443268/93‑1 (MQ=255)

                                                          cGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCaaa                              <  2:2176230/93‑1 (MQ=255)

                                                          cGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCaaa                              >  2:3730152/1‑93 (MQ=255)

                                                          cGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCaaa                              >  1:2753486/1‑93 (MQ=255)

                                                           gAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAAt                             >  1:2197309/1‑93 (MQ=255)

                                                            aaaCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAAtt                            >  2:3721174/1‑93 (MQ=255)

                                                            aaaCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAAtt                            >  1:1541531/1‑93 (MQ=255)

                                                            aaaCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAAtt                            <  1:2561496/93‑1 (MQ=255)

                                                             aaCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTa                           <  2:236336/93‑1 (MQ=255)

                                                              aCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAg                          <  1:3334359/93‑1 (MQ=255)

                                                              aCCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAg                          <  1:3736465/93‑1 (MQ=255)

                                                                 gTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAgcgc                       >  1:2529889/1‑93 (MQ=255)

                                                                 gTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAgcgc                       <  1:1830898/93‑1 (MQ=255)

                                                                 gTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGCGCAAATTAgcgc                       >  1:560138/1‑93 (MQ=255)

                                                                  tCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGcc                      >  2:3188082/1‑93 (MQ=255)

                                                                  tCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGcc                      >  1:304887/1‑93 (MQ=255)

                                                                  tCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGGGACGGTAAATCTGTTACACGGCCGCGCACTCGATGAAGGGCAAATTAGCGcc                      >  2:3310130/1‑93 (MQ=255)

                                                                   ttttATTCGTCCGTGAACAAAAGCCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCa                     <  2:973492/92‑1 (MQ=255)

                                                                    tttATTCGTCTGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCAt                    <  2:752312/93‑1 (MQ=255)

                                                                     ttattaGTCCGTGAACAAAAACCCCCTTCCGCATCGGTGCCGGCAAATCTGTTACCCGGCCGCGCACTCGACGAACGGCAAATTCGCGCCagg                   >  2:1749609/1‑91 (MQ=255)

                                                                     ttATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCAtt                   <  2:3029262/93‑1 (MQ=255)

                                                                     ttATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCAtt                   >  1:3122080/1‑93 (MQ=255)

                                                                      tATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATtg                  <  2:383657/93‑1 (MQ=255)

                                                                       aTTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATtgt                 >  2:3620934/1‑93 (MQ=255)

                                                                       aTTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATtgt                 >  1:3463702/1‑93 (MQ=255)

                                                                         tCGTCCGTGAACAAAACCCCCCTCCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGc               <  2:2411494/93‑1 (MQ=255)

                                                                          cGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCa              >  1:2701977/1‑93 (MQ=255)

                                                                           gTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCAt             >  1:603588/1‑93 (MQ=255)

                                                                           gTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCAt             >  1:4230153/1‑93 (MQ=255)

                                                                            tCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATc            <  1:864166/93‑1 (MQ=255)

                                                                            tCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATc            >  2:1158962/1‑93 (MQ=255)

                                                                              cGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTg          >  2:826480/1‑93 (MQ=255)

                                                                              cGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTg          >  2:83758/1‑93 (MQ=255)

                                                                                tGAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGt        <  2:4505157/93‑1 (MQ=255)

                                                                                 gAACAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGtt       <  2:2416385/93‑1 (MQ=255)

                                                                                  aaCAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGttt      >  1:3855714/1‑93 (MQ=255)

                                                                                  aaCAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGttt      >  2:711272/1‑93 (MQ=255)

                                                                                    cAAAAACCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTcc    <  1:96711/93‑1 (MQ=255)

                                                                                     aaaaaCCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCa   <  2:4015593/93‑1 (MQ=255)

TCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAATCTGTTACCCGGCCGCGCACTCGATGAAGGGCAAATTAGCGCCATTGTGCATCTGGTTTCCAG  >  NC_000913/2011662‑2011840

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