BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,474,826	+A	intergenic (+1/+44)	yjgL → / ← argI	hypothetical protein/ornithine carbamoyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,474,823	1	.	A	99.5%	708.8 / ‑3.7	189	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (113/75); total (114/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

Only 100 of 202 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑A‑AATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

  aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑a‑a                                                                                       >  2:5681232/1‑93 (MQ=255)

  aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑a‑a                                                                                       >  2:4759330/1‑93 (MQ=255)

  aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATAACTGCat‑a‑a                                                                                       >  2:1827695/1‑93 (MQ=255)

       aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑a‑aaataa                                                                                  >  2:801012/1‑89 (MQ=255)

         aGCGTTTAAATAAACTAAATCCTCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAAtg                                                                                 <  2:5772811/93‑1 (MQ=255)

         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAAtg                                                                                 <  1:4236265/93‑1 (MQ=255)

          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAAtgt                                                                                >  2:1424382/1‑93 (MQ=255)

           cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAAtgtg                                                                               >  1:6131420/1‑93 (MQ=255)

           cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAAtgtg                                                                               >  1:1626061/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGa                                                                              >  1:1426005/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGa                                                                              >  2:253061/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGa                                                                              >  2:2914406/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGa                                                                              >  1:3044535/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGAt                                                                             <  1:819455/93‑1 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGAt                                                                             >  1:3858957/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGAt                                                                             >  2:447948/1‑93 (MQ=255)

              ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATg                                                                            >  2:5022455/1‑93 (MQ=255)

               taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGc                                                                           >  2:6253102/1‑93 (MQ=255)

                   taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCggg                                                                       >  1:1390378/1‑93 (MQ=255)

                       cTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATgg                                                                   <  2:4276029/93‑1 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGt                                                                  >  1:6250025/1‑93 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGt                                                                  <  1:4354615/93‑1 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgt                                                              >  1:4654216‑M1/1‑91 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgta                                                             >  1:2387039‑M1/1‑90 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgta                                                             >  1:5368782‑M1/1‑90 (MQ=255)

                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtat                                                            <  1:4817115‑M1/93‑5 (MQ=255)

                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtat                                                            >  2:3749439‑M1/1‑89 (MQ=255)

                               aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatt                                                           <  1:5121459‑M1/93‑6 (MQ=255)

                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  1:6188448‑M1/93‑9 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  1:1874019‑M1/1‑84 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  1:3451470‑M1/93‑10 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  1:1908403‑M1/1‑84 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  1:2939848‑M1/93‑10 (MQ=255)

                                    aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      <  2:1877398‑M1/93‑11 (MQ=255)

                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    <  1:6267432‑M1/93‑13 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  1:4429171‑M1/1‑79 (MQ=255)

                                         ttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatc                                                 >  2:3085510‑M1/1‑78 (MQ=255)

                                          tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatct                                                >  1:2042671‑M1/1‑77 (MQ=255)

                                          tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatct                                                <  1:5864078‑M1/93‑17 (MQ=255)

                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  2:4417342‑M1/1‑76 (MQ=255)

                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  1:5151736‑M1/1‑76 (MQ=255)

                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               <  1:4759330‑M1/93‑18 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  2:5268792‑M1/1‑74 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  1:5730937‑M1/93‑20 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  1:352271‑M1/1‑72 (MQ=255)

                                                 aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                         >  2:1291615‑M1/1‑70 (MQ=255)

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        <  1:5476236‑M1/93‑25 (MQ=255)

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        >  2:3992435‑M1/1‑69 (MQ=255)

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        >  2:4247583‑M1/1‑69 (MQ=255)

                                                   cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       <  2:5049414‑M1/93‑26 (MQ=255)

                                                   cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       <  2:1093536‑M1/93‑26 (MQ=255)

                                                   cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       >  1:3298950‑M1/1‑68 (MQ=255)

                                                    aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      <  1:4944815‑M1/93‑27 (MQ=255)

                                                    aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagagat                                      >  2:3108042‑M1/1‑67 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:3019429‑M1/1‑66 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     <  1:6072187‑M1/93‑28 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  1:12057‑M1/1‑66 (MQ=255)

                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                   >  1:2584624‑M1/1‑64 (MQ=255)

                                                        gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  1:2812936‑M1/1‑63 (MQ=255)

                                                         aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 <  2:3291175‑M1/93‑32 (MQ=255)

                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               >  1:453326‑M1/1‑60 (MQ=255)

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             <  2:3254223‑M1/93‑36 (MQ=255)

                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:1566235‑M1/1‑58 (MQ=255)

                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            >  2:2209688‑M1/1‑57 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         >  2:1678804‑M1/1‑54 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  1:567675‑M1/93‑40 (MQ=255)

                                                                  ggTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        <  1:2367125‑M1/93‑41 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                       <  1:5217927‑M1/93‑42 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      >  1:5769616‑M1/1‑51 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      >  1:4140364‑M1/1‑51 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      <  1:3985816‑M1/93‑43 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     <  1:253125‑M1/93‑44 (MQ=255)

                                                                      gaagaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    >  1:4633476‑M1/1‑49 (MQ=255)

                                                                       aagaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                   <  1:1099085‑M1/93‑46 (MQ=255)

                                                                        agaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  2:3480821‑M1/1‑47 (MQ=255)

                                                                        agaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  <  1:2656261‑M1/93‑47 (MQ=255)

                                                                        agaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  <  1:2546745‑M1/93‑47 (MQ=255)

                                                                        agaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  <  2:5479762‑M1/93‑47 (MQ=255)

                                                                         gaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 <  2:4453000‑M1/93‑48 (MQ=255)

                                                                         gaaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  2:1774525‑M1/1‑46 (MQ=255)

                                                                          aaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcg                >  2:2350620‑M1/1‑45 (MQ=255)

                                                                          aaTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcg                >  2:4361207‑M1/1‑45 (MQ=255)

                                                                           aTAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               <  2:5368782‑M1/93‑50 (MQ=255)

                                                                            tAGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcg              >  1:5005995‑M1/1‑43 (MQ=255)

                                                                             aGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga             <  1:333980‑M1/93‑52 (MQ=255)

                                                                             aGACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga             <  2:646593‑M1/93‑52 (MQ=255)

                                                                              gACCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac            <  2:2117719‑M1/93‑53 (MQ=255)

                                                                               aCCTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacc           <  1:3790670‑M1/93‑54 (MQ=255)

                                                                                 cTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         <  1:1087848‑M1/93‑56 (MQ=255)

                                                                                 cTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         >  2:1722333‑M1/1‑38 (MQ=255)

                                                                                 cTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         >  1:5598209‑M1/1‑38 (MQ=255)

                                                                                 cTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         <  2:4487617‑M1/93‑56 (MQ=255)

                                                                                 cTATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         >  1:3168479‑M1/1‑38 (MQ=255)

                                                                                  tATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        >  2:210086‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc        >  2:5918926‑M1/1‑37 (MQ=255)

                                                                                     cACTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc     >  2:3187548‑M1/1‑34 (MQ=255)

                                                                                      aCTGCAT‑AGAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccg    >  2:6136611‑M1/1‑33 (MQ=255)

                                                                                       cTGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgc   >  1:1924148‑M1/1‑32 (MQ=255)

                                                                                        tGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  2:4985706‑M1/1‑31 (MQ=255)

                                                                                        tGCATAA‑AATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  1:1190270‑M1/1‑31 (MQ=255)

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑A‑AATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

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