BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	362,884	A→G	21.8%	L99S (TTA→TCA)	lacY ←	lactose permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	362,884	0	A	G	21.8%	314.4 / 88.8	165	L99S (TTA→TCA)	lacY	lactose permease
Reads supporting (aligned to +/- strand): ref base A (65/64); new base G (19/17); total (84/81)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 188 total aligned reads displayed.

ACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGAA  >  NC_000913/362785‑362982

 cTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTaa                                                                                                   >  1:980387/1‑100 (MQ=255)

  tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAAc                                                                                                  >  2:2761327/1‑100 (MQ=255)

  tttCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAAc                                                                                                  <  1:797670/100‑1 (MQ=255)

     ctcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAAAAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGt                                                                                               >  2:938428/1‑100 (MQ=255)

      tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                              >  1:701871/1‑100 (MQ=255)

      tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                              >  2:374513/1‑100 (MQ=255)

      tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                              >  1:3909472/1‑100 (MQ=255)

      tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                              <  1:2576137/100‑1 (MQ=255)

      tcAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                              >  1:1546543/1‑100 (MQ=255)

       cAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTgg                                                                                             >  1:3433036/1‑100 (MQ=255)

           aaaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCg                                                                                         <  2:3110406/100‑1 (MQ=255)

           aaaTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCg                                                                                         <  1:1207173/100‑1 (MQ=255)

            aaTGCCTCTACTGCTAGCGCACCGACGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGccc                                                                                          <  2:1632186/98‑1 (MQ=255)

             aTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGaa                                                                                       <  2:2758205/100‑1 (MQ=255)

             aTGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGaa                                                                                       <  2:3161251/100‑1 (MQ=255)

                ccTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGAt                                                                                    >  1:3404336/1‑100 (MQ=255)

                  tctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaa                                                                                  >  2:603983/1‑100 (MQ=255)

                   ctACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaa                                                                                 <  2:3154216/100‑1 (MQ=255)

                    tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaa                                                                                <  1:3695680/100‑1 (MQ=255)

                    tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaa                                                                                >  1:1169434/1‑100 (MQ=255)

                    tACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaa                                                                                >  1:1024340/1‑100 (MQ=255)

                     aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATaaaaa                                                                               >  2:728601/1‑100 (MQ=255)

                     aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaaa                                                                               <  1:226962/100‑1 (MQ=255)

                     aCTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATaaaaa                                                                               >  2:422087/1‑100 (MQ=255)

                      ctgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAAt                                                                              >  1:299074/1‑100 (MQ=255)

                       tgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATa                                                                             >  2:2349725/1‑100 (MQ=255)

                       tgctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATa                                                                             >  1:3191252/1‑100 (MQ=255)

                        gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATaa                                                                            >  1:710503/1‑100 (MQ=255)

                        gctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATaa                                                                            >  1:3962154/1‑100 (MQ=255)

                         ctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATaaa                                                                           <  2:1728880/100‑1 (MQ=255)

                         ctgGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATaaa                                                                           <  1:3750122/100‑1 (MQ=255)

                           ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGa                                                                         >  1:3177083/1‑100 (MQ=255)

                           ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGa                                                                         >  1:2818750/1‑100 (MQ=255)

                           ggCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGa                                                                         >  1:1988196/1‑100 (MQ=255)

                             cgcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAAc                                                                       >  1:1098340/1‑100 (MQ=255)

                              gcACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACg                                                                      >  1:3090769/1‑100 (MQ=255)

                                  cGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGgcgc                                                                  <  1:1964825/100‑1 (MQ=255)

                                   ggCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGAGGCCCGAAGATAAAAATAAAGAACGGCGCa                                                                 >  2:1511706/1‑100 (MQ=255)

                                   agCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCa                                                                 <  2:3246899/99‑1 (MQ=255)

                                      gTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAAc                                                              >  2:3566117/1‑100 (MQ=255)

                                       ttAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACa                                                             <  1:422090/100‑1 (MQ=255)

                                        tAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACAt                                                            >  1:152675/1‑100 (MQ=255)

                                         aaaaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATc                                                           >  1:2219217/1‑100 (MQ=255)

                                           ctcAAAAGCCTAGTTAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                         <  2:572890/98‑1 (MQ=255)

                                           aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                         <  1:218317/100‑1 (MQ=255)

                                           aaCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAc                                                         <  1:2390110/100‑1 (MQ=255)

                                             cAAAAGCCTAGTTAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCAcag                                                       >  1:572895/1‑98 (MQ=255)

                                                aaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACa                                                    <  1:2463744/100‑1 (MQ=255)

                                                aaGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACa                                                    <  1:2033201/100‑1 (MQ=255)

                                                  gcatAGATAAATACCATCAACAATCGATACTACTAAAATGTTGTATTGTAACAGTGGCCCTAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATg                                                  <  2:3468512/97‑1 (MQ=255)

                                                  gCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATg                                                  >  1:177913/1‑100 (MQ=255)

                                                   ccTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGc                                                 <  1:654695/100‑1 (MQ=255)

                                                      aGATAAATACCACAAACAATCGATCCTACTAGAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCgg                                              >  1:2781414/1‑100 (MQ=255)

                                                       gATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGt                                             <  1:389117/100‑1 (MQ=255)

                                                         tAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaa                                           >  2:2765784/1‑100 (MQ=255)

                                                          aaaTACCACCAACAATCGATCCTACTAAAATGTTGGATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaat                                          <  2:422546/100‑1 (MQ=255)

                                                           aaTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGtaata                                         >  1:1821420/1‑100 (MQ=255)

                                                            aTACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTGACATGCCGGtaataa                                        >  1:4072351/1‑100 (MQ=255)

                                                             tACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTaataat                                       <  2:3464269/100‑1 (MQ=255)

                                                               ccaccaACAATCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATc                                      >  2:2312570/1‑99 (MQ=255)

                                                                caccaACAATCGATCCTGCTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATGca                                    >  2:2544035/1‑100 (MQ=255)

                                                                caccaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCca                                    >  1:3328676/1‑100 (MQ=255)

                                                                 accaACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                   >  1:2544931/1‑100 (MQ=255)

                                                                 accaACAATCGATCATACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                                   <  2:3695360/100‑1 (MQ=255)

                                                                   caacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcag                                 >  1:930905/1‑100 (MQ=255)

                                                                   caacaaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcag                                 >  2:2884316/1‑100 (MQ=255)

                                                                    aacaaTCGATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                                <  2:4037159/100‑1 (MQ=255)

                                                                    aacaaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcagc                                <  2:3316705/100‑1 (MQ=255)

                                                                       aaTCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAgg                             >  2:749446/1‑100 (MQ=255)

                                                                           gATCCTACTAAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAtt                         <  2:2813458/100‑1 (MQ=255)

                                                                            aTCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTAttt                        >  1:3218703/1‑100 (MQ=255)

                                                                             tCCTACTAAATTGTTGTATTGTAACAGTGGCCCGAAGATGCAGGTAAAGAACGGCGCAAACATCACTATCATGCCTGTAATAATCCACAGCAGGTATTTg                       <  2:2903823/100‑1 (MQ=255)

                                                                             tCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTg                       <  2:299071/100‑1 (MQ=255)

                                                                             tCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTg                       <  2:1024334/100‑1 (MQ=255)

                                                                               ctactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcg                     >  2:2429539/1‑100 (MQ=255)

                                                                                tactaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                    <  1:1355909/100‑1 (MQ=255)

                                                                                tactaAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc                    >  1:1467074/1‑100 (MQ=255)

                                                                                  ctaAAATGTTGTATTGTGACAGTGGCCCGAGGATAAAAATAAAGAACGGCGCAAACATCACTAACATCCCGGTAATAATCCACAGAAGGTATTTCCGCa                   >  2:1162895/1‑99 (MQ=255)

                                                                                  ctaAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                  >  1:3795826/1‑100 (MQ=255)

                                                                                  ctaAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                  >  1:2714737/1‑100 (MQ=255)

                                                                                  ctaAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                  <  1:1019049/100‑1 (MQ=255)

                                                                                  ctaAAATGTTGCATTATAACAGAAGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGACGGTAATAATCGGCAGCAGGTATTTtcg                     >  2:816705/1‑97 (MQ=255)

                                                                                   taAAATGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAg                  <  2:1555887/99‑1 (MQ=255)

                                                                                   taAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGc                 <  2:3471758/100‑1 (MQ=255)

                                                                                    aaaaTGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGcc                >  2:2626945/1‑100 (MQ=255)

                                                                                      aaTGTTGTATTGTGACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCg              <  2:1907100/100‑1 (MQ=255)

                                                                                       aTGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGa             >  1:1946283/1‑100 (MQ=255)

                                                                                         gttgtATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGt           >  1:3014212/1‑100 (MQ=255)

                                                                                          ttgtattgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt          <  2:1983634/100‑1 (MQ=255)

                                                                                          ttctattgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt          <  1:1849048/97‑1 (MQ=255)

                                                                                            gtattgtgACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTg        >  1:1628123/1‑100 (MQ=255)

                                                                                            gtattgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTg        <  1:1294255/100‑1 (MQ=255)

                                                                                             tattgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGt       >  2:791405/1‑100 (MQ=255)

                                                                                              attgtaACAGTAGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc      <  2:2040875/100‑1 (MQ=255)

                                                                                               ttgtgACAGTGGCGCGAAGATAAAAATAAAGACCGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTc      <  2:1680499/99‑1 (MQ=255)

                                                                                               ttgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     <  1:1077276/100‑1 (MQ=255)

                                                                                               ttgtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCa     <  2:3806737/100‑1 (MQ=255)

                                                                                                 gtaACAGTGGCCCGAAGATAGAAATACAGCACGGCGCAAACATCACTAACATGCCGGTAATACTCCACAGCAGCTATTTGCGAAGCCCGAGTTTGTCAg    >  2:1240229/1‑99 (MQ=255)

                                                                                                 gtaACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGa   >  1:1760326/1‑100 (MQ=255)

                                                                                                  taACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGaa  >  1:3573436/1‑100 (MQ=255)

ACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTTTGTCAGAA  >  NC_000913/362785‑362982

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