BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,091,676	A→G	Q53R (CAG→CGG)	yqgF →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,091,676	0	A	G	99.7%	918.6 / ‑5.0	287	Q53R (CAG→CGG)	yqgF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (101/184); total (101/186)
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.60e-01

Only 100 of 320 total aligned reads displayed.

TTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTA  >  NC_000913/3091586‑3091765

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                         >  1:9439372/1‑93 (MQ=255)

ttACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcggc                                                                                         >  2:8082923/1‑93 (MQ=255)

  aCCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCg                                                                                       >  1:3847620/1‑93 (MQ=255)

    cGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGa                                                                                     >  2:6174309/1‑93 (MQ=255)

     ggCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGAc                                                                                    >  2:1390743/1‑93 (MQ=255)

     ggCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGAc                                                                                    <  2:9771770/93‑1 (MQ=255)

      gCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTggcgggcgggcg                                                                                   >  1:7374543/1‑93 (MQ=255)

      gCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACg                                                                                   <  2:5537782/93‑1 (MQ=255)

        aCCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGACGaa                                                                                 >  2:492087/1‑93 (MQ=255)

         ccGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGaaa                                                                                <  1:8468637/93‑1 (MQ=255)

         ccGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGaaa                                                                                <  1:6980802/93‑1 (MQ=255)

          cGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGgaaa                                                                               >  1:8740904/1‑92 (MQ=255)

          cGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAat                                                                               >  1:655409/1‑93 (MQ=255)

            cTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAatca                                                                             >  2:9439729/1‑93 (MQ=255)

            cTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatca                                                                             <  2:9862690/93‑1 (MQ=255)

            cTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCCGGCGAAatca                                                                             >  2:1736389/1‑93 (MQ=255)

             tCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatcat                                                                            <  1:3290957/93‑1 (MQ=255)

             ccGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatcat                                                                            <  1:2182206/92‑1 (MQ=255)

              cGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatcatc                                                                           <  2:6874224/93‑1 (MQ=255)

              cGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAatcatc                                                                           <  1:6304141/93‑1 (MQ=255)

              cGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCCGACCAAatcatc                                                                           >  2:8247412/1‑93 (MQ=255)

               gCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGGGGGAAACAtcg                                                                          >  1:1887570/1‑93 (MQ=255)

                cccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgt                                                                         <  2:3283567/93‑1 (MQ=255)

                cccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgt                                                                         <  2:4037258/93‑1 (MQ=255)

                 ccTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgtc                                                                        <  1:8675896/93‑1 (MQ=255)

                  cTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGGCGCAATCAtcctcg                                                                       >  1:6284533/1‑93 (MQ=255)

                  cTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAtcgtcg                                                                       <  1:111623/93‑1 (MQ=255)

                   tttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGGAAAAATCATCgg                                                                      >  1:7794215/1‑93 (MQ=255)

                   tttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGCAATCATCGTCgg                                                                      >  2:8593620/1‑93 (MQ=255)

                    ttGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGt                                                                     <  1:297150/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGtt                                                                    <  1:1432659/93‑1 (MQ=255)

                     tGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGtt                                                                    <  2:6590083/93‑1 (MQ=255)

                      gCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGGGGCGGCCGG‑CGAAATCATCGTCCGtttt                                                                  >  2:3840669/1‑92 (MQ=255)

                        cTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTCGGGTTGc                                                                 >  1:1507459/1‑93 (MQ=255)

                         tGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGCCGAAATCATCGTAGGTTTGcc                                                                >  2:6742322/1‑93 (MQ=255)

                         tGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAAACATCGTCGGTTTGcc                                                                >  1:427286/1‑93 (MQ=255)

                           cAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGc                                                              <  2:1105340/93‑1 (MQ=255)

                             aTTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCTTCGGTTTGCCGCTg                                                            >  2:2723178/1‑93 (MQ=255)

                             aTTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAAACATCGTCGGTTTGCCGCTg                                                            >  1:317554/1‑93 (MQ=255)

                              atAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGa                                                           <  1:7569681/92‑1 (MQ=255)

                                 aaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAata                                                        <  1:8845401/93‑1 (MQ=255)

                                 aaGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAata                                                        <  2:6426216/93‑1 (MQ=255)

                                   gCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATg                                                      <  2:4255207/93‑1 (MQ=255)

                                   gCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATg                                                      >  2:3181410/1‑93 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGa                                                    <  1:8347496/93‑1 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGa                                                    <  2:6049892/93‑1 (MQ=255)

                                     acaGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGa                                                    <  1:6157598/93‑1 (MQ=255)

                                      caGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAAACATCGTCGGGTTTCCCCTGAATATGGag                                                   >  1:8266843/1‑92 (MQ=255)

                                      caGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGAc                                                   <  2:7759067/93‑1 (MQ=255)

                                      caGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGAc                                                   <  1:5793496/93‑1 (MQ=255)

                                      caGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGAc                                                   <  2:2903244/93‑1 (MQ=255)

                                       aGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGTCCCTGAATAAGGGCg                                                  >  1:5297378/1‑93 (MQ=255)

                                          aCGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGCCGGCAGCTCCACTCCTTCGCCGGTTTGCGGCTGAACATGGTCGGct                                               <  2:5789358/93‑2 (MQ=255)

                                           cGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCAc                                              <  2:6107016/93‑1 (MQ=255)

                                           cGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCAc                                              <  1:7722634/93‑1 (MQ=255)

                                           cGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCAc                                              <  2:2787938/93‑1 (MQ=255)

                                            ggTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGGCGGGCGAAATCATCGTCGGGTTGCCGCTGAATATGGAAGGCAcc                                             >  2:3627263/1‑93 (MQ=255)

                                              tACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGa                                           <  1:10001352/93‑1 (MQ=255)

                                                    gtaCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATAATCCCCGGCTTGGCGCCGAAACTCGGCGGCCCCCAGCAccc                                     >  1:9841266/3‑93 (MQ=255)

                                                     gACTGGACCANTATCGAGCGTTTNATTGAGGAGTGGCGGCCAGACGAAATCAGCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCa                                    <  1:9546677/93‑1 (MQ=255)

                                                     gACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCa                                    <  1:248292/93‑1 (MQ=255)

                                                     gACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCa                                    <  2:5256606/93‑1 (MQ=255)

                                                      aCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCAt                                   <  1:5923149/93‑1 (MQ=255)

                                                      aCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCAt                                   <  2:8850961/93‑1 (MQ=255)

                                                      aCTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGTCGCTGGATATTGACGGCACCGCGCAACCAt                                   >  1:9143484/1‑93 (MQ=255)

                                                       cTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAAACATCGTCGGTTTGCCGCCGAAAATGGAAGGCACCGCGCAACCact                                  >  1:2041423/1‑91 (MQ=255)

                                                        tGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAATCATCGTCGGTTTTCCCCTGGATATTGACCGCACCGCGCAACCAttt                                 >  2:4805906/1‑92 (MQ=255)

                                                        tGGAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATATGGACGGGACCGAGCAACCactt                                 >  1:3716772/1‑90 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGa                                <  2:2151329/93‑1 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGa                                <  2:4373652/93‑1 (MQ=255)

                                                         ggAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGCCGCCGAATATGGACGGCACCCAGCAAACAttta                                >  2:347642/1‑91 (MQ=255)

                                                          gAACATTATCGAGCGTTTACTGAAAGAGTGGCGGCCCGGCGAAATCATCGTCGGTTTTCCGCTGTATATGGACGGCACCCAGCAACCATTTAc                               >  2:8315609/1‑93 (MQ=255)

                                                          gAACATCATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCAGCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGAc                               <  2:4456979/93‑1 (MQ=255)

                                                           aaCATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACt                              <  1:8888224/93‑1 (MQ=255)

                                                           aaCATTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACt                              <  1:5872906/93‑1 (MQ=255)

                                                            aCATTGTCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCCGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTg                             <  1:2934258/93‑1 (MQ=255)

                                                             cATTATCGAGCGTGTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGc                            <  1:976225/93‑1 (MQ=255)

                                                              aTTATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGcc                           <  2:5563786/93‑1 (MQ=255)

                                                               ttATCGAGCGTTTACTGAAAGAGTGGCGGCCGGGCGAAAACATCGTCGGTTTGCCGCTGAATATGGACGGGACCGAGCAACCATTGACTGCCa                          >  2:1653590/1‑93 (MQ=255)

                                                               ttATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCa                          <  1:7315336/93‑1 (MQ=255)

                                                               ttATCGAGCGTTGACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCa                          <  2:1362805/93‑1 (MQ=255)

                                                                tATCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGGTTGCCGCTGAATATGGACGGCACCGAGCAAACATTGACCGCCag                         >  1:1214952/1‑93 (MQ=255)

                                                                 ttCGAGCGTTGACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGCCTGCCaga                        <  1:5385397/92‑1 (MQ=255)

                                                                  tCGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCagag                       <  1:2730987/93‑1 (MQ=255)

                                                                   cGAGCGTTTACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAgc                      <  1:4674916/93‑1 (MQ=255)

                                                                     aGCGTTGAGTGAAAGAGTGGCCGAGGAACGAAATCATCGTCAGTTTGCCGCTGAATATGGTCGGCACCGAGAAACCATTGACTGCCAGAgcgc                    <  2:8623838/93‑1 (MQ=255)

                                                                        gTTTACTGAAAGAGTGGCGGCCGGACGAAAATATCGTCGGGTTGCCTCTGAATATGGACGGCACCCAGCCACCACTGACTGCCACAGCGCGTa                 >  2:8055770/1‑93 (MQ=255)

                                                                         tttACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGGACCGAGCAAACACTGCATGCCCGAGCGCGTaa                >  2:5475282/1‑93 (MQ=255)

                                                                          ttACTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTaaa               <  2:5889961/93‑1 (MQ=255)

                                                                            aCTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAAtt             <  2:2693258/93‑1 (MQ=255)

                                                                            aCTGAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGGAAAtt             >  2:929780/1‑93 (MQ=255)

                                                                               gAAAGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGGAACCATTGACTCCCAGAGCGCGTAAAttttc          >  2:2099135/1‑91 (MQ=255)

                                                                                aaaGAGTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGcc         <  2:9185332/93‑1 (MQ=255)

                                                                                  agagTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCaa       <  1:3479951/93‑1 (MQ=255)

                                                                                  agagTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCaa       <  2:7360216/93‑1 (MQ=255)

                                                                                    agTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAAcc     <  1:3916493/93‑1 (MQ=255)

                                                                                     gTGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAAACg    >  2:8882314/1‑93 (MQ=255)

                                                                                      tGGCGGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGt   <  1:7505053/93‑1 (MQ=255)

                                                                                       ggcggcCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTa  <  1:5533088/93‑1 (MQ=255)

                                                                                       ggcggcCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTa  <  2:7891735/93‑1 (MQ=255)

TTACCGGCACCGCTCGCCCTTTGCCTGCAATTAAAGCACAGGACGGTACGCCGGACTGGAACATTATCGAGCGTTTACTGAAAGAGTGGCAGCCGGACGAAATCATCGTCGGTTTGCCGCTGAATATGGACGGCACCGAGCAACCATTGACTGCCAGAGCGCGTAAATTTGCCAACCGTA  >  NC_000913/3091586‑3091765

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