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	98.9%	678.1 / ‑0.9	185	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base A (135/47); total (137/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 206 total aligned reads displayed.

TAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AA‑ATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474732‑4474910

tAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑a                                                                                         >  2:342496/1‑93 (MQ=255)

tAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑a                                                                                         >  1:3419042/1‑93 (MQ=255)

 aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa                                                                                        >  2:5951169/1‑93 (MQ=255)

   tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa‑aa                                                                                     >  2:1565955/1‑92 (MQ=255)

    gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGCGTGGAGGAATAGACCTATCACTACat‑aa‑aat                                                                                    >  2:5420820/1‑91 (MQ=255)

     aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGCAGAATAGACCTATCACTGCat‑aa‑aata                                                                                   >  2:6033515/1‑90 (MQ=255)

     aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa‑aata                                                                                   >  2:5015168/1‑90 (MQ=255)

      aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa‑aataa                                                                                  >  2:4675667/1‑89 (MQ=255)

        aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAAtg                                                                                 >  1:45386/1‑93 (MQ=255)

         gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAAtgt                                                                                >  1:2355491/1‑93 (MQ=255)

          cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAAtgtg                                                                               >  1:2825353/1‑93 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGa                                                                              >  1:594960/1‑93 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGa                                                                              <  1:4755129/93‑1 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGa                                                                              <  2:2938951/93‑1 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGa                                                                              >  1:3234041/1‑93 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGa                                                                              >  1:1722121/1‑93 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGa                                                                              >  2:5898847/1‑93 (MQ=255)

            tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGAt                                                                             >  2:4765403/1‑93 (MQ=255)

            tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGAt                                                                             >  1:2250184/1‑93 (MQ=255)

              taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGc                                                                           >  1:100725/1‑93 (MQ=255)

               aaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGcc                                                                          <  1:6134878/93‑1 (MQ=255)

                aataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCg                                                                         >  1:2982064/1‑93 (MQ=255)

                 ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCgg                                                                        >  1:5563144/1‑93 (MQ=255)

                 ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGAGGAAGAATAGAACTATCACTGCATAAA‑ATAATGTGATGCCgg                                                                        >  1:3165079/1‑93 (MQ=255)

                  taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCggg                                                                       >  2:3995772/1‑93 (MQ=255)

                  taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGCGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCggg                                                                       >  2:4121345/1‑93 (MQ=255)

                   aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGa                                                                      <  2:4341809/93‑1 (MQ=255)

                   aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGa                                                                      >  2:4198923/1‑93 (MQ=255)

                    aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGAt                                                                     >  1:3932139/1‑93 (MQ=255)

                       tAAATCCACCGGAGAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGt                                                                  <  1:3123003/93‑1 (MQ=255)

                       tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGt                                                                  >  2:5066301/1‑93 (MQ=255)

                        aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGtt                                                                 >  2:2493677/1‑93 (MQ=255)

                         aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGttt                                                                >  1:2338699/1‑93 (MQ=255)

                         aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGttt                                                                <  2:5853518/93‑1 (MQ=255)

                           tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgt                                                              >  1:5504509‑M1/1‑91 (MQ=255)

                            ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgta                                                             <  1:330918‑M1/93‑4 (MQ=255)

                             cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtat                                                            >  2:1787278‑M1/1‑89 (MQ=255)

                                 ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGGTTgtatttcc                                                        >  1:5642346‑M1/1‑85 (MQ=255)

                                  gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  2:5327009‑M1/1‑84 (MQ=255)

                                   aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      >  2:898866‑M1/1‑83 (MQ=255)

                                          aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               <  1:356423‑M1/93‑18 (MQ=255)

                                          aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AACATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  2:4304981‑M1/1‑76 (MQ=255)

                                             tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            >  2:1728596‑M1/1‑73 (MQ=255)

                                              ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  1:5422106‑M1/1‑72 (MQ=255)

                                              ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           >  2:1809059‑M1/1‑72 (MQ=255)

                                              ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttccccgcatctttata                                           <  2:6049234‑M1/93‑22 (MQ=255)

                                               taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                          >  2:1169834‑M1/1‑71 (MQ=255)

                                               taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                          >  2:3944705‑M1/1‑71 (MQ=255)

                                                  cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       <  2:3812306‑M1/93‑26 (MQ=255)

                                                   aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      >  1:2342827‑M1/1‑67 (MQ=255)

                                                   aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      >  2:3648206‑M1/1‑67 (MQ=255)

                                                   aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      >  2:5531806‑M1/1‑67 (MQ=255)

                                                    aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:4169251‑M1/1‑66 (MQ=255)

                                                    aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:5083436‑M1/1‑66 (MQ=255)

                                                     tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  2:4252231‑M1/1‑65 (MQ=255)

                                                     tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  1:2996833‑M1/1‑65 (MQ=255)

                                                     tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  2:4898934‑M1/1‑65 (MQ=255)

                                                     tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  1:3624455‑M1/1‑65 (MQ=255)

                                                        aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 >  1:858954‑M1/1‑62 (MQ=255)

                                                         acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                >  2:1997794‑M1/1‑61 (MQ=255)

                                                         acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                >  2:858132‑M1/1‑61 (MQ=255)

                                                           acTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              <  2:4513392‑M1/93‑35 (MQ=255)

                                                           acTTTGGGGGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              >  1:750361‑M1/1‑59 (MQ=255)

                                                            cTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             <  1:4638897‑M1/93‑36 (MQ=255)

                                                            cTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:976056‑M1/1‑58 (MQ=255)

                                                            cTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:1019782‑M1/1‑58 (MQ=255)

                                                            cTTTGGGAGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:3787086‑M1/1‑58 (MQ=255)

                                                             tttGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            >  1:5695039‑M1/1‑57 (MQ=255)

                                                              ttGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt                           <  2:3095212‑M1/93‑38 (MQ=255)

                                                                gggTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:3044809‑M1/93‑40 (MQ=255)

                                                                   tGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      >  1:711958‑M1/1‑51 (MQ=255)

                                                                   tGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      >  1:4952552‑M1/1‑51 (MQ=255)

                                                                    ggAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  1:4435566‑M1/1‑50 (MQ=255)

                                                                    ggAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttagagcgatagcaattatttactgag                     >  1:4301212‑M1/1‑50 (MQ=255)

                                                                     gaagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                    <  2:2799261‑M1/93‑45 (MQ=255)

                                                                     gaagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                    <  2:3855829‑M1/93‑45 (MQ=255)

                                                                     gaagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                    >  2:3338939‑M1/1‑49 (MQ=255)

                                                                      aagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg                   >  1:4488460‑M1/1‑48 (MQ=255)

                                                                         aaTAGACCTATCACTGCATGAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcg                <  2:3738197‑M1/93‑49 (MQ=255)

                                                                          aTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  2:4216875‑M1/1‑44 (MQ=255)

                                                                          aTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  2:1299029‑M1/1‑44 (MQ=255)

                                                                           tAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttaatgagcgtcgcg              >  1:1089295‑M1/1‑43 (MQ=255)

                                                                            aGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcga             >  2:432283‑M1/1‑42 (MQ=255)

                                                                             gACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            >  1:3883577‑M1/1‑41 (MQ=255)

                                                                              aCCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc           >  2:3684646‑M1/1‑40 (MQ=255)

                                                                              aCCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc           >  1:1890799‑M1/1‑40 (MQ=255)

                                                                               ccTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          <  2:3178786‑M1/93‑55 (MQ=255)

                                                                                cTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccat         >  1:411691‑M1/1‑38 (MQ=255)

                                                                                cTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccat         <  2:3562327‑M1/93‑56 (MQ=255)

                                                                                 tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc        >  1:3846635‑M1/1‑37 (MQ=255)

                                                                                 tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagagtcgcgcccatc        >  2:700911‑M1/1‑37 (MQ=255)

                                                                                  aTCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatca       <  2:219833‑M1/93‑58 (MQ=255)

                                                                                     aCTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccg    >  2:3462560‑M1/1‑33 (MQ=255)

                                                                                     aCTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccg    >  1:5918916‑M1/1‑33 (MQ=255)

                                                                                      cTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   >  1:2966180‑M1/1‑32 (MQ=255)

                                                                                      cTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   <  2:5683864‑M1/93‑62 (MQ=255)

                                                                                       tGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  1:5051003‑M1/1‑31 (MQ=255)

                                                                                       tGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  >  2:2015308‑M1/1‑31 (MQ=255)

                                                                                       tGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  >  1:1605777‑M1/1‑31 (MQ=255)

                                                                                       tGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacagct  >  2:2280588‑M1/1‑31 (MQ=255)

TAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AA‑ATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474732‑4474910

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