BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	781,326	A→T	6.9%	intergenic (+104/‑43)	lysY → / → lysZ	tRNA‑Lys/tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,326	0	A	T	6.9%	453.8 / 18.8	131	intergenic (+104/‑43)	lysY/lysZ	tRNA‑Lys/tRNA‑Lys
Reads supporting (aligned to +/- strand): ref base A (66/56); new base T (6/3); total (72/59)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

Only 100 of 141 total aligned reads displayed.

TTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGAACG  >  NC_000913/781180‑781472

tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                                                                                                   <  4:290681/149‑1 (MQ=255)

tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                                                                                                   <  3:724799/149‑1 (MQ=255)

tttAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaa                                                                                                                                                    <  2:422461/148‑1 (MQ=255)

    aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAtt                                                                                                                                               <  3:551199/149‑1 (MQ=255)

     tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                                                                                              <  4:289078/149‑1 (MQ=255)

      cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGAt                                                                                                                                                <  1:445780/146‑1 (MQ=255)

      cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                                                                                              >  2:327963/1‑148 (MQ=255)

      cAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTg                                                                                                                                              >  1:598126/1‑148 (MQ=255)

       aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                                                                                                   <  4:5828/142‑1 (MQ=255)

         ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATGATACTCAGATGATAAGCGGGAAGGATAACGTTGCGCCAGCAACGGCCCGGAGGGAGAGCGAAGCGAGTCATCCAGCACGACCCACCAATCTTAAAGATTGGccc                                                                                                                                          >  4:238867/1‑149 (MQ=17)

         ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                                                                                                   <  4:644658/140‑1 (MQ=255)

         ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTaaa                                                                                                                                                   >  3:644658/1‑140 (MQ=255)

         ttGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGcc                                                                                                                                           <  2:583773/148‑1 (MQ=255)

            gTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCg                                                                                                                                        <  2:620460/148‑1 (MQ=255)

             tCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGa                                                                                                                                       >  2:239531/1‑148 (MQ=255)

               gCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGt                                                                                                                                     >  2:740559/1‑148 (MQ=255)

                cAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggt                                                                                                                                            <  1:136052/140‑2 (MQ=255)

                   gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTTAGATTggtgg                                                                                                                                          >  3:852457/1‑136 (MQ=255)

                   gTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTTAGATTggtgg                                                                                                                                          <  4:852457/139‑4 (MQ=255)

                     tCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAAt                                                                                                                               >  1:91265/1‑148 (MQ=255)

                        aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                                                                                           >  4:377052/1‑149 (MQ=255)

                        aaTCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCtt                                                                                                                            >  1:461457/1‑148 (MQ=255)

                          tCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                                                                                          <  1:546619/148‑1 (MQ=255)

                           ccTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAg                                                                                                                        <  4:828775/149‑1 (MQ=255)

                            cTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAgg                                                                                                                       >  4:706185/1‑149 (MQ=255)

                             tGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCATAAAGATTGGCCCCGAGTAAAAATCTTACAgg                                                                                                                       <  1:259017/148‑1 (MQ=255)

                              gCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGt                                                                                                                      >  2:199156/1‑148 (MQ=255)

                               cACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTaa                                                                                                                    <  4:70807/149‑1 (MQ=255)

                                   aCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                                                                                          <  2:35018/139‑1 (MQ=255)

                                   aCCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTc                                                                                                                          >  1:35018/1‑139 (MQ=255)

                                     ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCg                                                                                                             >  4:53694/1‑149 (MQ=255)

                                     ccaccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCg                                                                                                             >  3:185060/1‑149 (MQ=255)

                                       accaGTTTTAGCATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGGGAAGCGAGTCATCCTGCACGACCCACCAATCTGAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTa                                                                                                           >  4:586071/1‑149 (MQ=255)

                                       accaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTTAAGATTggtg                                                                                                                                           >  1:251809/1‑116 (MQ=255)

                                        ccaGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTAt                                                                                                          >  4:392046/1‑149 (MQ=255)

                                             tttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGt                                                                                                      >  2:841947/1‑148 (MQ=255)

                                             tttAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGt                                                                                                      >  1:526370/1‑148 (MQ=255)

                                                   aTCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAac‑ac                                                                                                                <  4:220033/132‑1 (MQ=255)

                                                    tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAACGACCCATACGGGtgttacc                                                                                                <  2:543690/148‑7 (MQ=21)

                                                    tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTc                                                                                                     >  4:285245/1‑142 (MQ=255)

                                                    tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTc                                                                                                     <  3:285245/142‑1 (MQ=255)

                                                    tCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCt                                                                                              >  3:249463/1‑149 (MQ=255)

                                                      aaaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTc                                                                                             <  1:146771/148‑1 (MQ=255)

                                                       aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTa                                                                                                                                                     >  2:607074/1‑92 (MQ=255)

                                                       aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTa                                                                                                                                                     <  1:607074/92‑1 (MQ=255)

                                                       aaCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATggg                                                                                                       >  1:521913/1‑137 (MQ=255)

                                                        aCTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAg                                                                                           <  1:417895/148‑1 (MQ=255)

                                                         cTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGt                                                                                          <  1:740559/148‑1 (MQ=255)

                                                          tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                                                                                           >  3:749117/1‑115 (MQ=255)

                                                          tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCttt                                                                                                                           <  4:749117/115‑1 (MQ=255)

                                                          tCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTg                                                                                        >  4:666399/1‑149 (MQ=255)

                                                           cAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTgg                                                                                       <  3:579570/149‑1 (MQ=255)

                                                            aGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCGTAAAGATTGGCCCCGAGTAAAACTCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGt                                                                                      <  3:588837/149‑1 (MQ=255)

                                                             gATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGt                                                                                      <  1:327963/148‑1 (MQ=255)

                                                               tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGt                                                                                                   <  3:615079/133‑1 (MQ=255)

                                                               tGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGACGt                                                                                                   >  4:615079/1‑133 (MQ=255)

                                                                 ttAAGCGTGAATGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGc                                                                                 <  3:161322/149‑1 (MQ=255)

                                                                 ttAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTagag                                                                                  >  1:86261/1‑148 (MQ=255)

                                                                  tAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                <  3:32455/149‑1 (MQ=255)

                                                                    aGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                                                                              >  4:444935/1‑149 (MQ=255)

                                                                     gCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGt                                                                              <  1:199156/148‑1 (MQ=255)

                                                                       gTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                           <  4:317541/149‑1 (MQ=255)

                                                                         gAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt                                                                         <  3:392046/149‑1 (MQ=255)

                                                                           aGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                        <  2:862990/148‑1 (MQ=255)

                                                                           aGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                        <  1:841947/148‑1 (MQ=255)

                                                                              aTAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                                                                     >  2:337526/1‑148 (MQ=255)

                                                                                    ttGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGt                                                                                                   <  3:260839/112‑1 (MQ=255)

                                                                                    ttGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGt                                                                                                   >  4:260839/1‑112 (MQ=255)

                                                                                                cGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGtt                                                   <  2:72380/148‑1 (MQ=255)

                                                                                                 ggCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAt                                                                   >  3:309367/1‑131 (MQ=255)

                                                                                                 ggCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAt                                                                   <  4:309367/131‑1 (MQ=255)

                                                                                                       tAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                            <  2:526370/148‑1 (MQ=255)

                                                                                                       tAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                            <  2:461457/148‑1 (MQ=255)

                                                                                                               gcgaagcgaGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTc                                                                                                     >  4:769137/1‑83 (MQ=255)

                                                                                                                    gcgaGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacc                               >  1:662090/1‑148 (MQ=255)

                                                                                                                    gcgaGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacc                               >  2:770734/1‑148 (MQ=255)

                                                                                                                          cATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTa                         >  2:769621/1‑148 (MQ=255)

                                                                                                                            tCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                       <  3:732347/100‑1 (MQ=255)

                                                                                                                            tCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGAGTCGTTAGCTCAGGTGGTAGAGCAGTTGACttt                                                                       >  4:732347/1‑100 (MQ=255)

                                                                                                                                gCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAggtgg                  <  3:706185/149‑1 (MQ=255)

                                                                                                                                    gACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTa               <  1:762723/148‑1 (MQ=255)

                                                                                                                                      cccACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACt             <  2:152520/148‑1 (MQ=255)

                                                                                                                                         accaATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTa         >  3:585430/1‑149 (MQ=255)

                                                                                                                                          ccaATCTTTAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTa         >  1:158363/1‑148 (MQ=255)

                                                                                                                                          ccaATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTa         >  1:811710/1‑148 (MQ=255)

                                                                                                                                          ccaATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTa         >  1:912903/1‑148 (MQ=255)

                                                                                                                                          ccaATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTa         >  1:120131/1‑148 (MQ=255)

                                                                                                                                            aaTCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTagag      >  3:816973/1‑149 (MQ=255)

                                                                                                                                            aaTCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTaga       >  1:66291/1‑148 (MQ=255)

                                                                                                                                              tCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGaa    >  4:462523/1‑149 (MQ=255)

                                                                                                                                              tCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGaa    <  4:638272/149‑1 (MQ=255)

                                                                                                                                              tCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGaa    >  4:262527/1‑149 (MQ=255)

                                                                                                                                                ttAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTgg           <  4:118335/140‑1 (MQ=255)

                                                                                                                                                ttAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTgg           >  3:118335/1‑140 (MQ=255)

                                                                                                                                                ttAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGAACg  >  3:377110/1‑149 (MQ=255)

                                                                                                                                                ttAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGAACg  >  4:800307/1‑149 (MQ=255)

                                                                                                                                                 tAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTaa                        >  2:483316/1‑126 (MQ=255)

                                                                                                                                                 tAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTaa                        <  1:483316/126‑1 (MQ=255)

                                                                                                                                                 tAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGAACg  <  2:442758/148‑1 (MQ=255)

                                                                                                                                                 tAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGAACg  >  2:788406/1‑148 (MQ=255)

TTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCAAACTCAGATGTTAAGCGTGAAGGATAACGTTGCGCCAGCAACGGCCCGTAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATCTTAAAGATTGGCCCCGAGTAAAAATCTTTCAGGTAAC‑ACCCGTATGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAGGTGGTTACTGGTAGAGAACG  >  NC_000913/781180‑781472

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