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.8%	1544.8 / ‑4.3	414	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/0); new base A (219/194); total (220/194)
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.87e-01

Only 100 of 463 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AA‑ATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTA  >  NC_000913/4474731‑4474909

ttAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat                                                                                          >  1:2760938/1‑93 (MQ=255)

ttAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat                                                                                          <  2:7555092/93‑1 (MQ=255)

ttAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat                                                                                          <  1:23720/93‑1 (MQ=255)

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

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

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

     gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa‑aat                                                                                   >  1:2579668/1‑91 (MQ=255)

     gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa‑aat                                                                                   >  2:5028601/1‑91 (MQ=255)

     gAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aa‑aat                                                                                   <  1:3976093/93‑3 (MQ=255)

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

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

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

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

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

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

              ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATg                                                                           >  2:3245223/1‑93 (MQ=255)

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

               taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGc                                                                          <  2:1637558/93‑1 (MQ=255)

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

                   taaaCTCAATCCACCGGCATATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCggg                                                                      <  2:8961008/93‑1 (MQ=255)

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

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

                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGGAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATg                                                                   <  1:5987299/93‑1 (MQ=255)

                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATg                                                                   >  1:9783892/1‑93 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AACATAATGTGATGCCGGGATGGt                                                                 >  1:4049635/1‑93 (MQ=255)

                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGtt                                                                <  1:6789312/93‑1 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgt                                                             <  2:7337311‑M1/93‑3 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAGTAGAACACTTTGGGTGGAAGCATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgta                                                            <  1:8978427‑M1/93‑4 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgta                                                            >  1:4479501‑M1/1‑90 (MQ=255)

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

                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtat                                                           <  1:7513662‑M1/93‑5 (MQ=255)

                                ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtattt                                                         >  2:4346405‑M1/1‑87 (MQ=255)

                                  ggAACATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcc                                                       <  1:9057456‑M1/93‑9 (MQ=255)

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

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

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttccc                                                      <  1:1142572‑M1/93‑10 (MQ=255)

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

                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                    >  1:8051673‑M1/1‑82 (MQ=255)

                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccgga                                                   >  2:8069086‑M1/1‑81 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                 >  2:2033503‑M1/1‑79 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                 >  1:9457536‑M1/1‑79 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                 >  1:3839208‑M1/1‑79 (MQ=255)

                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                 <  2:8745849‑M1/93‑15 (MQ=255)

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

                                            cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                             >  1:2157047‑M1/1‑75 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                            <  1:9402299‑M1/93‑20 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                            <  2:1163198‑M1/93‑20 (MQ=255)

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

                                              tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                           <  2:3973458‑M1/93‑21 (MQ=255)

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

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                          <  2:766633‑M1/93‑22 (MQ=255)

                                                taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                         <  1:77465‑M1/93‑23 (MQ=255)

                                                taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                         <  1:8836557‑M1/93‑23 (MQ=255)

                                                taGCAATAGAACACTTTGGGAGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                         >  1:4941827‑M1/1‑71 (MQ=255)

                                                 aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                        <  2:6490850‑M1/93‑24 (MQ=255)

                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                       >  1:6424387‑M1/1‑69 (MQ=255)

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

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                    <  2:9606127‑M1/93‑28 (MQ=255)

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

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

                                                        gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGGTTgtatttcccggcatctttatagcgatagca                                 >  1:5667833‑M1/1‑63 (MQ=255)

                                                         aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                <  2:7810595‑M1/93‑32 (MQ=255)

                                                         aaCACTTGTGGTGCAAGAATAGACCTATCACTGCAT‑AAGATAATGTGGTGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                <  1:3255337‑M1/93‑32 (MQ=255)

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

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

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

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

                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                           <  2:3785425‑M1/93‑37 (MQ=255)

                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                           <  2:9486235‑M1/93‑37 (MQ=255)

                                                              tttGGGTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatt                           <  2:9111012‑M1/93‑37 (MQ=255)

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

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                        >  2:5654501‑M1/1‑54 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattacttac                        <  2:7178180‑M1/93‑40 (MQ=255)

                                                                  ggTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                       >  1:5148942‑M1/1‑53 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                      <  1:4157042‑M1/93‑42 (MQ=255)

                                                                   gTGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                      <  1:3392581‑M1/93‑42 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                     <  2:1886747‑M1/93‑43 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                    <  2:2803457‑M1/93‑44 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                    <  2:5191655‑M1/93‑44 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                    <  2:3491908‑M1/93‑44 (MQ=255)

                                                                      gaagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagt                   >  2:1775149‑M1/1‑49 (MQ=255)

                                                                       aagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGCTGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtg                  >  1:1677114‑M1/1‑48 (MQ=255)

                                                                       aagaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGATTgtatttcccggcatctttatagcgatagcaattatttactgagtg                  >  2:1028312‑M1/1‑48 (MQ=255)

                                                                         gaaTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtc                >  2:6099098‑M1/1‑46 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgc              >  1:10016482‑M1/1‑44 (MQ=255)

                                                                             aGACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcga            <  2:7113964‑M1/93‑52 (MQ=255)

                                                                              gACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac           <  2:8001053‑M1/93‑53 (MQ=255)

                                                                              gACCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgac           >  2:9901155‑M1/1‑41 (MQ=255)

                                                                               aCCTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgacc          >  2:5864190‑M1/1‑40 (MQ=255)

                                                                                 cTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat        >  1:2921102‑M1/1‑38 (MQ=255)

                                                                                 cTATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccat        <  1:9601473‑M1/93‑56 (MQ=255)

                                                                                  tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc       >  1:2006507‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc       >  2:4727896‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc       >  1:1460614‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc       >  1:8229113‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc       >  1:8493310‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatc       >  1:8395772‑M1/1‑37 (MQ=255)

                                                                                    tCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac     <  2:1326482‑M1/93‑59 (MQ=255)

                                                                                    tCACTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcac     <  1:7906277‑M1/93‑59 (MQ=255)

                                                                                       cTGCATAAA‑ATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgc  >  2:9188724‑M1/1‑32 (MQ=255)

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AA‑ATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTA  >  NC_000913/4474731‑4474909

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