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.4%	539.6 / ‑3.7	153	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (87/65); total (88/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 162 total aligned reads displayed.

TAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474732‑4474910

  aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaa                                                                                      >  1:2323806/1‑93 (MQ=255)

    gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaat                                                                                    >  1:277316/1‑91 (MQ=255)

     aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaata                                                                                   <  2:2397605/93‑4 (MQ=255)

        aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtg                                                                                 <  2:1481846/93‑1 (MQ=255)

        aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtg                                                                                 >  2:690532/1‑93 (MQ=255)

         gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgt                                                                                >  1:2136436/1‑93 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:112372/1‑93 (MQ=255)

           gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:572443/1‑93 (MQ=255)

            tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  1:1136394/1‑93 (MQ=255)

            tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:3705834/1‑93 (MQ=255)

             ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAGCACTGTGGGTGGAAGAATAGACCTATCTCTGCATAAAATAATGTGATg                                                                            <  1:116862/93‑1 (MQ=255)

             ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATg                                                                            >  1:2534454/1‑93 (MQ=255)

             ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATg                                                                            >  1:4516183/1‑93 (MQ=255)

              taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGc                                                                           <  1:2760126/93‑1 (MQ=255)

               aaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGcc                                                                          >  1:4650205/1‑93 (MQ=255)

                 ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCgg                                                                        >  2:1261497/1‑93 (MQ=255)

                 ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCgg                                                                        <  2:531158/93‑1 (MQ=255)

                   aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGa                                                                      <  1:4573536/93‑1 (MQ=255)

                    aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  2:602433/1‑93 (MQ=255)

                       tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  >  1:2049984/1‑93 (MQ=255)

                       tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  <  2:1682674/93‑1 (MQ=255)

                       tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  >  2:723743/1‑93 (MQ=255)

                        aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 >  2:4511908/1‑93 (MQ=255)

                         aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGttt                                                                >  1:789358/1‑93 (MQ=255)

                          aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtttt                                                               >  2:3953431/1‑92 (MQ=255)

                            ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                             <  2:2087834‑M1/93‑4 (MQ=255)

                            ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                             >  1:2386256‑M1/1‑90 (MQ=255)

                             cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            >  1:1413915‑M1/1‑89 (MQ=255)

                              aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                           <  2:3374637‑M1/93‑6 (MQ=255)

                              aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                           >  1:2780745‑M1/1‑88 (MQ=255)

                               ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtattt                                                          >  1:2924690‑M1/1‑87 (MQ=255)

                               ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtattt                                                          >  1:1680397‑M1/1‑87 (MQ=255)

                                cGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                         <  1:4664615‑M1/93‑8 (MQ=255)

                                 ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:2232631‑M1/1‑85 (MQ=255)

                                 ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  2:3723324‑M1/1‑85 (MQ=255)

                                 ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  2:1324248‑M1/93‑9 (MQ=255)

                                 ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  1:2232630‑M1/1‑85 (MQ=255)

                                  gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  2:2415231‑M1/93‑10 (MQ=255)

                                   aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      <  1:2673390‑M1/93‑11 (MQ=255)

                                   aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      <  1:3278480‑M1/93‑11 (MQ=255)

                                     aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    >  2:1368023‑M1/1‑81 (MQ=255)

                                      aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGTTGGTTTgtatttcccggct                                                   >  1:1597348‑M1/1‑80 (MQ=255)

                                      aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                   <  1:4407574‑M1/93‑14 (MQ=255)

                                      aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                   >  2:4159878‑M1/1‑80 (MQ=255)

                                          aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  2:2528334‑M1/1‑76 (MQ=255)

                                          aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               <  2:583285‑M1/93‑18 (MQ=255)

                                          aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  2:598407‑M1/1‑76 (MQ=255)

                                           cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              <  1:3691054‑M1/93‑19 (MQ=255)

                                           cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              >  2:3073767‑M1/1‑75 (MQ=255)

                                           cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              >  1:2423045‑M1/1‑75 (MQ=255)

                                            gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  1:1257783‑M1/1‑74 (MQ=255)

                                            gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  2:1276963‑M1/1‑74 (MQ=255)

                                            gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  1:2412716‑M1/93‑20 (MQ=255)

                                            gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  1:604422‑M1/93‑20 (MQ=255)

                                             tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            <  2:1129306‑M1/93‑21 (MQ=255)

                                              ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  1:3018678‑M1/93‑22 (MQ=255)

                                                aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                         >  2:2605331‑M1/1‑70 (MQ=255)

                                                 gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        <  1:2748881‑M1/93‑25 (MQ=255)

                                                  cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                       <  1:3784377‑M1/93‑26 (MQ=255)

                                                   aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      >  2:3751318‑M1/1‑67 (MQ=255)

                                                    aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  1:3299056‑M1/1‑66 (MQ=255)

                                                       gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  <  1:4065657‑M1/93‑31 (MQ=255)

                                                        aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 <  1:1733351‑M1/93‑32 (MQ=255)

                                                        aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 <  2:2089422‑M1/93‑32 (MQ=255)

                                                         acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                >  2:2944466‑M1/1‑61 (MQ=255)

                                                           acTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              <  2:3114116‑M1/93‑35 (MQ=255)

                                                            cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             <  1:4363307‑M1/93‑36 (MQ=255)

                                                             tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            <  1:4039905‑M1/93‑37 (MQ=255)

                                                             tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            >  1:1236178‑M1/1‑57 (MQ=255)

                                                             tttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                            >  1:4273279‑M1/1‑57 (MQ=255)

                                                               tGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                          >  2:4630499‑M1/1‑55 (MQ=255)

                                                                gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  1:197838‑M1/93‑40 (MQ=255)

                                                                gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         >  2:1178017‑M1/1‑54 (MQ=255)

                                                                gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:1125997‑M1/93‑40 (MQ=255)

                                                                 ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        <  2:2322350‑M1/93‑41 (MQ=255)

                                                                 ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        >  1:3140169‑M1/1‑53 (MQ=255)

                                                                    ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  2:2510398‑M1/1‑50 (MQ=255)

                                                                    ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  2:2510399‑M1/1‑50 (MQ=255)

                                                                     gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                    >  2:1274176‑M1/1‑49 (MQ=255)

                                                                     gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTggatttcccggcatctttatagcgatagcaattatttactgagt                    >  1:3195543‑M1/1‑49 (MQ=255)

                                                                      aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                   >  2:1884684‑M1/1‑48 (MQ=255)

                                                                       agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgt                  >  1:3517189‑M1/1‑47 (MQ=255)

                                                                        gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 <  1:1377050‑M1/93‑48 (MQ=255)

                                                                        gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 >  1:1378857‑M1/1‑46 (MQ=255)

                                                                        gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                 <  2:3299752‑M1/93‑48 (MQ=255)

                                                                          aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               >  2:863189‑M1/1‑44 (MQ=255)

                                                                          aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc               <  2:3800639‑M1/93‑50 (MQ=255)

                                                                            aGACCTATCACTGCATAAAATAATGTGATGCCGGGATTGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga             <  1:1621536‑M1/93‑52 (MQ=255)

                                                                             gACCTATCACAGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac            <  1:175000‑M1/93‑53 (MQ=255)

                                                                               ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca          <  1:3627434‑M1/93‑55 (MQ=255)

                                                                               ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca          <  1:3239409‑M1/93‑55 (MQ=255)

                                                                               ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacca          >  1:1082037‑M1/1‑39 (MQ=255)

                                                                                cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat         <  2:3321878‑M1/93‑56 (MQ=255)

                                                                                  aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttctcggcatctttatagcgatagcaattatttactgagtgtcgcgaccatca       >  1:3453603‑M1/1‑36 (MQ=255)

                                                                                   tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac      <  2:4050593‑M1/93‑59 (MQ=255)

                                                                                   tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac      <  1:3821052‑M1/93‑59 (MQ=255)

                                                                                    cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcacc     >  1:1011401‑M1/1‑34 (MQ=255)

                                                                                      cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgc   <  1:755513‑M1/93‑62 (MQ=255)

                                                                                       tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  >  2:4344897‑M1/1‑31 (MQ=255)

                                                                                       tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct  <  1:2395477‑M1/93‑63 (MQ=255)

TAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474732‑4474910

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