BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,954,504	T→C	intergenic (+8/+44)	ilvA → / ← ilvY	threonine deaminase/IlvY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,954,504	0	T	C	99.5%	591.6 / NA	193	intergenic (+8/+44)	ilvA/ilvY	threonine deaminase/IlvY
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (91/101); minor base A (1/0); total (92/101)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

TCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGC  >  NC_000913/3954414‑3954593

tCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCTCGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGc                                                                                         <  2:2931973/93‑1 (MQ=255)

 cGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGcc                                                                                        <  1:935985/93‑1 (MQ=255)

  gAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCt                                                                                       >  1:4906858/1‑93 (MQ=255)

  gAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCt                                                                                       >  1:4077891/1‑93 (MQ=255)

   aaaCCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTg                                                                                      >  1:3671332/1‑93 (MQ=255)

    accccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGa                                                                                     <  2:2801082/91‑1 (MQ=255)

    aaCCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAAAGCCTGa                                                                                     >  1:5298418/1‑93 (MQ=255)

     aCCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCCGGTTAGGGGAAAACGACTGAt                                                                                    >  1:5431812/1‑93 (MQ=255)

      cccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATa                                                                                   <  2:36705/93‑1 (MQ=255)

      cccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATa                                                                                   <  2:5274483/93‑1 (MQ=255)

      cccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATa                                                                                   >  1:3743402/1‑93 (MQ=255)

      cccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATa                                                                                   >  1:6121125/1‑93 (MQ=255)

       ccGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAAACACTGATAg                                                                                  >  2:3289617/1‑93 (MQ=255)

          gCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAgcgc                                                                               <  1:227982/93‑1 (MQ=255)

          gCTGAATGAGCCGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAACGCCTGATAgcgc                                                                               >  1:1394323/1‑93 (MQ=255)

           cTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGcgct                                                                              <  2:1210218/93‑1 (MQ=255)

            tGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGcgctt                                                                             >  2:1681079/1‑93 (MQ=255)

               aTGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGcgcttcgc                                                                          <  2:1736378/93‑1 (MQ=255)

                tGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGcgcttcgct                                                                         <  1:3312051/93‑1 (MQ=255)

                 gAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGcgcttcgctt                                                                        >  1:1840812/1‑93 (MQ=255)

                  aGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTa                                                                       <  1:4380063/93‑1 (MQ=255)

                  aGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTa                                                                       >  2:3014459/1‑93 (MQ=255)

                    cTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATc                                                                     >  2:2220965/1‑93 (MQ=255)

                    cTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATc                                                                     <  1:1507779/93‑1 (MQ=255)

                    cTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATc                                                                     <  2:3160372/93‑1 (MQ=255)

                       ggCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAgg                                                                  <  1:4931513/93‑1 (MQ=255)

                       ggCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAgg                                                                  <  1:2570614/93‑1 (MQ=255)

                       ggCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAgg                                                                  >  1:1962116/1‑93 (MQ=255)

                        gCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGc                                                                 <  1:3987799/93‑1 (MQ=255)

                        gCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGc                                                                 <  1:2599166/93‑1 (MQ=255)

                         cTACGATTTCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATATCGCTTCGCTTCTAAGGcc                                                                >  2:6136883/1‑93 (MQ=255)

                           aCGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTa                                                              >  2:3803639/1‑93 (MQ=255)

                            cGATTGCCACGACGAAACCAATATCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAACCGCCTGATAGCGCTTCGCTTATCAGGCCTAc                                                             <  2:2395978/93‑1 (MQ=255)

                             gATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTAcc                                                            <  1:2053642/93‑1 (MQ=255)

                             gATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTAcc                                                            <  1:356443/93‑1 (MQ=255)

                             gATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTAcc                                                            <  2:5797672/93‑1 (MQ=255)

                             gATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTAcc                                                            <  1:975136/93‑1 (MQ=255)

                              aTTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTAccc                                                           <  2:620155/93‑1 (MQ=255)

                               ttGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCcg                                                          <  2:5604572/93‑1 (MQ=255)

                               ttGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCcg                                                          <  2:2370078/93‑1 (MQ=255)

                               ttGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCcg                                                          >  1:3087225/1‑93 (MQ=255)

                               ttGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCCTCGCTTATCAAGCCTACCcg                                                          >  1:1752387/1‑93 (MQ=255)

                                tGCCACGACGAAACCAATACCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCcgc                                                         <  2:4714825/93‑1 (MQ=255)

                                    acgacgAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGa                                                     >  1:1110949/1‑93 (MQ=255)

                                      gacgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACGCGCGCGACa                                                   >  2:2543111/1‑93 (MQ=255)

                                        cgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAAc                                                 <  2:4549218/93‑1 (MQ=255)

                                             ccAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAAAGCCTGATAGCGCTTCGCTTATCAGGCCTAACCGCGCGACAACGTCAt                                            >  2:2715955/1‑93 (MQ=255)

                                               aaTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCAttt                                          >  2:2907365/1‑93 (MQ=255)

                                                aTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTtg                                         <  1:2177093/93‑1 (MQ=255)

                                                aTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCGTATCAGGCCTACCCGCGCGACATCGTCAtgtg                                         >  1:6302342/1‑93 (MQ=255)

                                                  aaCCCGGCGTTCAGGTTTTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTtgtg                                       <  2:1357776/93‑1 (MQ=255)

                                                   aCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTgg                                      >  1:1420201/1‑93 (MQ=255)

                                                   aCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTgg                                      <  2:4933669/93‑1 (MQ=255)

                                                      cGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTc                                   <  2:6057233/93‑1 (MQ=255)

                                                      cGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTc                                   >  2:5945627/1‑93 (MQ=255)

                                                       ggCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCg                                  <  1:5158328/93‑1 (MQ=255)

                                                           ttCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaa                              >  2:720904/1‑93 (MQ=255)

                                                           ttCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaa                              >  2:4387106/1‑93 (MQ=255)

                                                            tCAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaaa                             >  2:690901/1‑93 (MQ=255)

                                                             cAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaaaa                            >  2:5477552/1‑93 (MQ=255)

                                                             cAGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCaaaa                            <  1:2947470/93‑1 (MQ=255)

                                                              aGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAAt                           <  2:4079743/93‑1 (MQ=255)

                                                              aGGTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAAt                           <  2:6304/93‑1 (MQ=255)

                                                                gTTCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCt                         <  2:771637/93‑1 (MQ=255)

                                                                 ttCTTTTTGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCtt                        <  2:2743473/93‑1 (MQ=255)

                                                                    tttttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCa                     <  2:4193437/93‑1 (MQ=255)

                                                                     ttttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAg                    >  1:4705489/1‑93 (MQ=255)

                                                                     ttttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAg                    <  2:5283967/93‑1 (MQ=255)

                                                                     ttttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAg                    <  2:2266586/93‑1 (MQ=255)

                                                                      tttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGa                   >  1:4846168/1‑93 (MQ=255)

                                                                      tttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGa                   <  1:2880743/93‑1 (MQ=255)

                                                                       ttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGaa                  <  1:426758/93‑1 (MQ=255)

                                                                       ttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGaa                  <  2:1919070/93‑1 (MQ=255)

                                                                       ttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGaa                  <  2:5386817/93‑1 (MQ=255)

                                                                       ttGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGaa                  <  1:5589679/93‑1 (MQ=255)

                                                                       ttGGCGGGTTAGGGAAAAACGCCTGATAGAGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGGTAGGCAAAATCTTCCAGaa                  >  2:5515984/1‑93 (MQ=39)

                                                                        tGGCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAAt                 >  1:103109/1‑93 (MQ=255)

                                                                         ggCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATg                <  1:143757/93‑1 (MQ=255)

                                                                          gCGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGc               >  1:5794639/1‑93 (MQ=255)

                                                                           cGGGTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGcc              >  1:1563327/1‑93 (MQ=255)

                                                                            gggTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCt             <  2:1821209/93‑1 (MQ=255)

                                                                             ggTTAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTc            >  2:4871096/1‑93 (MQ=255)

                                                                               ttAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCaa          >  1:2457763/1‑93 (MQ=255)

                                                                               ttAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCaa          >  2:2981008/1‑93 (MQ=255)

                                                                                tAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAt         <  1:5652946/93‑1 (MQ=255)

                                                                                tAGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAt         >  1:3015072/1‑93 (MQ=255)

                                                                                 aGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt        <  2:4631219/93‑1 (MQ=255)

                                                                                 aGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt        >  2:4252483/1‑93 (MQ=255)

                                                                                 aGGGAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAAtt        <  1:5220960/93‑1 (MQ=255)

                                                                                  gggAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa       <  2:5089420/93‑1 (MQ=255)

                                                                                  gggAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa       <  2:5154577/93‑1 (MQ=255)

                                                                                  gggAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTa       >  1:3405112/1‑93 (MQ=255)

                                                                                   ggAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAg      <  2:1240787/93‑1 (MQ=255)

                                                                                    gAAAAACGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGc     >  1:1228264/1‑93 (MQ=255)

                                                                                    gAAAAACGACTGATAGCGCTTCGCTTATCAGGCTTACCCGCGCGACCACGCCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGc     >  2:882567/1‑93 (MQ=255)

                                                                                     aaaaaCGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCg    <  2:1060144/93‑1 (MQ=255)

                                                                                     aaaaaCGCCAGAGAGCGCTTCGCTTATCAGGCCGACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCg    <  1:2694836/93‑1 (MQ=255)

                                                                                      aaaaCGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCgg   <  1:4365327/93‑1 (MQ=255)

                                                                                      aaaaCGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCgg   >  2:2152884/1‑93 (MQ=255)

                                                                                       aaaCGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGc  >  2:5264274/1‑93 (MQ=255)

TCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGC  >  NC_000913/3954414‑3954593

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: `‑`