BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,472,602	A→G	Y56C (TAC→TGC)	b4252 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,472,602	0	A	G	99.4%	1673.6 / ‑2.6	528	Y56C (TAC→TGC)	b4252	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (271/254); total (272/256)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 553 total aligned reads displayed.

CACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCAAC  >  NC_000913/4472512‑4472691

  cacGTTACGGCAGCACCGCCAACGGGCAGGCACGATATCGACGGCGATCGACTGTTTTATCTTCTCTCGGAAGATATGACCGAGCCGTGCGag                                                                                       <  1:8984450/93‑2 (MQ=255)

  cacGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGGGCGaa                                                                                       >  1:3357944/1‑93 (MQ=255)

    cGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGc                                                                                     <  2:3351254/93‑1 (MQ=255)

     gtgaCGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCt                                                                                    <  2:360652/90‑1 (MQ=255)

     gTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCt                                                                                    >  2:7272971/1‑93 (MQ=255)

     gTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCt                                                                                    >  1:5286582/1‑93 (MQ=255)

        aCGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGc                                                                                 <  1:1351100/93‑1 (MQ=255)

           gCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGt                                                                              <  2:7844366/93‑1 (MQ=255)

             aGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGc                                                                            <  2:9961342/93‑1 (MQ=255)

             aGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGc                                                                            <  1:8297787/93‑1 (MQ=255)

              gAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCGGTGCGAAGCTCGCGGTGCg                                                                           >  1:3467641/1‑93 (MQ=255)

              gAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCg                                                                           >  2:5088111/1‑93 (MQ=255)

              gAAACGCCAAACGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCg                                                                           >  2:8267428/1‑93 (MQ=255)

                 aCGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAg                                                                        <  2:5764038/93‑1 (MQ=255)

                 aCGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAg                                                                        <  1:654453/93‑1 (MQ=255)

                  cGCCAAAGGGCAAGCACGATCTCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGt                                                                       <  2:5397428/93‑1 (MQ=255)

                   gCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGATCGCCGTGCGGAGTa                                                                      >  2:1739910/1‑93 (MQ=255)

                      aaaGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCa                                                                   <  2:3106323/93‑1 (MQ=255)

                      aaaGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGCAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCa                                                                   <  1:6929520/93‑1 (MQ=255)

                         gggCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGc                                                                <  2:8316648/93‑1 (MQ=255)

                          ggCAAGCACGATATCGAAGGCAATGGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGcc                                                               <  2:3984987/93‑1 (MQ=255)

                          ggCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGcc                                                               >  2:4743921/1‑93 (MQ=255)

                          ggCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGcc                                                               <  2:9647282/93‑1 (MQ=255)

                             aaGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGc                                                            <  2:4561909/93‑1 (MQ=255)

                             aaGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGc                                                            <  2:1849922/93‑1 (MQ=255)

                              aGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCt                                                           >  1:5511967/1‑93 (MQ=255)

                              aGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCt                                                           >  2:6153210/1‑93 (MQ=255)

                              aGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCt                                                           <  2:2555890/93‑1 (MQ=255)

                                 aCGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATc                                                        >  1:3826066/1‑93 (MQ=255)

                                 aCGATATCGAAGGCAATCGACTGTTCCATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCTGCTATc                                                        >  2:8968560/1‑93 (MQ=255)

                                  cGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGACCGCTATCt                                                       >  1:4732114/1‑93 (MQ=255)

                                  cGATATCGAAGGCAATCGACTGTTTTATCATATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCt                                                       >  1:2694613/1‑93 (MQ=255)

                                   gATATCGAAGGTAATCGACTGTTTTATCTTATCACGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCtt                                                      <  2:6282238/93‑1 (MQ=255)

                                   gATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCtt                                                      >  1:9966023/1‑93 (MQ=255)

                                    atatCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTg                                                     >  2:9065283/1‑93 (MQ=255)

                                    atatCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTg                                                     >  2:9175881/1‑93 (MQ=255)

                                     tatCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGa                                                    <  1:7165628/93‑1 (MQ=255)

                                     tatCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGa                                                    >  2:2281426/1‑93 (MQ=255)

                                       tCGAAGCCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACa                                                  <  1:8089820/93‑1 (MQ=255)

                                         gAAGGCAATCGACTGTTTTATCTGATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACAtt                                                <  2:5532111/93‑1 (MQ=255)

                                          aaGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTc                                               <  2:2742930/93‑1 (MQ=255)

                                           aGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCa                                              >  1:9162020/1‑93 (MQ=255)

                                           aGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCa                                              <  2:3812970/93‑1 (MQ=255)

                                            ggCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAg                                             >  1:10029434/1‑93 (MQ=255)

                                              agaTCGGCTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGAt                                           <  2:5828976/91‑1 (MQ=255)

                                               aaTCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGAtt                                          <  2:7607014/93‑1 (MQ=255)

                                                 tCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATtgt                                        >  1:3703857/1‑93 (MQ=255)

                                                 tCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATtgt                                        >  1:4546818/1‑93 (MQ=255)

                                                 tCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATtgt                                        <  2:6638710/93‑1 (MQ=255)

                                                  cGACTGTTTTATCTTAACTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATtgtg                                       <  1:3047018/93‑1 (MQ=255)

                                                    aCTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGtt                                     <  1:7671662/93‑1 (MQ=255)

                                                      tGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCTTGCCCGCTATCTTGACATTCAGATTGTGTTaa                                   >  1:5962657/1‑93 (MQ=255)

                                                      tGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTaa                                   >  2:8756897/1‑93 (MQ=255)

                                                      tGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTaa                                   >  1:1870241/1‑93 (MQ=255)

                                                       gTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTaaa                                  <  2:3267820/93‑1 (MQ=255)

                                                       gTTTTATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTaaa                                  >  1:5852854/1‑93 (MQ=255)

                                                        ttttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTaaaa                                 <  1:6398581/93‑1 (MQ=255)

                                                        ttttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTaaaa                                 >  1:3961578/1‑93 (MQ=255)

                                                         tttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAg                                >  2:6778709/1‑93 (MQ=255)

                                                         tttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAg                                <  2:9401587/93‑1 (MQ=255)

                                                         tttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAg                                >  1:2159228/1‑93 (MQ=255)

                                                         tttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAg                                <  1:1518887/93‑1 (MQ=255)

                                                         tttATCTTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAg                                >  1:62353/1‑93 (MQ=255)

                                                         tttATCCTATCTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTTTCTTGACATTCAGATTGTGTTAAAAg                                <  1:317118/93‑1 (MQ=255)

                                                          ttatcttatcTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAgg                               <  1:4444156/93‑1 (MQ=255)

                                                          ttatcttatcTCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAgg                               <  1:4184075/93‑1 (MQ=255)

                                                           tatcttatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGt                              >  2:4479501/1‑93 (MQ=255)

                                                            atcttatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTc                             <  2:1168493/93‑1 (MQ=255)

                                                            atcttatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTc                             <  2:8406824/93‑1 (MQ=255)

                                                             tcttatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCa                            <  1:2024978/93‑1 (MQ=255)

                                                             tcttatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCa                            <  2:5319643/93‑1 (MQ=255)

                                                             tcttatctCGGAAGATATGACCGAGCCGCGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCa                            >  1:4995142/1‑93 (MQ=255)

                                                              cttatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAg                           >  1:9451258/1‑93 (MQ=255)

                                                                tatctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGa                         <  2:6502148/93‑1 (MQ=255)

                                                                 atctCGGAAGATATGCCCGACCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGaa                        <  2:5271702/93‑1 (MQ=255)

                                                                 atctCGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGaa                        <  1:6823490/93‑1 (MQ=255)

                                                                    tcGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGc                     <  1:2842939/93‑1 (MQ=255)

                                                                    tcGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGc                     >  2:4840506/1‑93 (MQ=255)

                                                                    tcGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGc                     >  1:7962724/1‑93 (MQ=255)

                                                                     cGGAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCa                    <  2:9775986/93‑1 (MQ=255)

                                                                       gAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGGCATTCAGATTGTGTTAAAAGGACAGGAAGGCGTg                  >  1:869514/1‑93 (MQ=255)

                                                                       gAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATg                  <  1:961622/93‑1 (MQ=255)

                                                                       gAAGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATg                  >  1:8851249/1‑93 (MQ=255)

                                                                         aGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGAc                >  2:1556116/1‑93 (MQ=255)

                                                                         aGATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGAc                >  2:7496302/1‑93 (MQ=255)

                                                                          gATATGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGAcc               >  1:6116391/1‑93 (MQ=255)

                                                                           atatGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCt              >  2:4300198/1‑93 (MQ=255)

                                                                           atatGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCt              <  1:2433089/93‑1 (MQ=255)

                                                                             atGACCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTc            >  1:1604081/1‑93 (MQ=255)

                                                                               ggccGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAg          <  1:7445/91‑1 (MQ=255)

                                                                                aCCGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGc         >  2:7694301/1‑93 (MQ=255)

                                                                                 ccGAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCa        <  2:2438990/93‑1 (MQ=255)

                                                                                 ccGAGCCGTGCGAAGCTCGCCATGCGGAGGACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGTGCAGGAAGGCATGACCTTCAGcc        >  2:6954966/1‑92 (MQ=255)

                                                                                   gAGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACg      <  2:7997771/93‑1 (MQ=255)

                                                                                    aGCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGc     <  1:8461121/93‑1 (MQ=255)

                                                                                     gCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCa    <  2:4232860/93‑1 (MQ=255)

                                                                                     gCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCa    <  2:8500946/93‑1 (MQ=255)

                                                                                     gCCGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCa    >  1:1170644/1‑93 (MQ=255)

                                                                                      ccGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCaa   >  1:5412857/1‑93 (MQ=255)

                                                                                       cGTGCGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCAAc  <  1:2246605/93‑1 (MQ=255)

CACACGTTACGGCAGAAACGCCAAAGGGCAAGCACGATATCGAAGGCAATCGACTGTTTTATCTTATCTCGGAAGATATGACCGAGCCGTACGAAGCTCGCCGTGCGGAGTACCATGCCCGCTATCTTGACATTCAGATTGTGTTAAAAGGTCAGGAAGGCATGACCTTCAGCACGCAAC  >  NC_000913/4472512‑4472691

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