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.2%	2011.2 / 4.5	496	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/1); new base A (260/231); total (263/233)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 542 total aligned reads displayed.

ATATTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑T‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGC  >  NC_000913/4474728‑4474912

      aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCa‑t‑aaaata                                                                                     <  1:7710276/96‑4 (MQ=255)

      aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCa‑t‑aaaata                                                                                     <  1:600756/96‑4 (MQ=255)

      aTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCa‑t‑aaaata                                                                                     >  1:3069984/1‑93 (MQ=255)

       tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCa‑t‑aaaataa                                                                                    >  1:3315274/1‑92 (MQ=255)

         aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAAtg                                                                                   <  1:9084158/96‑1 (MQ=255)

         aaCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAAtg                                                                                   >  2:13423262/1‑96 (MQ=255)

          aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAAtgt                                                                                  <  1:5205182/96‑1 (MQ=255)

          aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACGGCA‑TAAAATAAtgt                                                                                  >  1:3763860/1‑96 (MQ=255)

           cAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAAtgtg                                                                                 <  2:11110585/96‑1 (MQ=255)

           cAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAAtgtg                                                                                 >  2:12062003/1‑96 (MQ=255)

             gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGAt                                                                               <  1:2119172/96‑1 (MQ=255)

              cctttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATg                                                                              >  1:11860761/3‑96 (MQ=255)

              cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATg                                                                              <  2:3115586/96‑1 (MQ=255)

               gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGATTAGACCTATCACTGCA‑TAAAATAATGTGATGc                                                                             >  1:4188375/1‑96 (MQ=255)

               gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGc                                                                             <  1:11081107/96‑1 (MQ=255)

               gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGc                                                                             <  2:8131512/96‑1 (MQ=255)

                 ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCg                                                                           >  2:4437667/1‑96 (MQ=255)

                 ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCg                                                                           >  2:6261640/1‑96 (MQ=255)

                  taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCgg                                                                          <  1:4959266/96‑1 (MQ=255)

                   aaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCggg                                                                         <  2:9258384/96‑1 (MQ=255)

                   aaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCggg                                                                         >  2:1490179/1‑96 (MQ=255)

                    aataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGa                                                                        >  2:13966438/1‑96 (MQ=255)

                     ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGAt                                                                       <  2:13084948/96‑1 (MQ=255)

                       aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATgg                                                                     >  1:14135744/1‑96 (MQ=255)

                       aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATgg                                                                     >  2:3352381/1‑96 (MQ=255)

                       aaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATgg                                                                     >  1:2251567/1‑96 (MQ=255)

                        aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGt                                                                    >  2:955853/1‑96 (MQ=255)

                         aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGtt                                                                   <  2:1231882/96‑1 (MQ=255)

                          cTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGttt                                                                  >  2:1252299/1‑96 (MQ=255)

                          cTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGttt                                                                  >  1:9885987/1‑96 (MQ=255)

                            aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGCTGGTTTgt                                                                <  2:12963152‑M1/96‑3 (MQ=255)

                            aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgt                                                                <  2:11740822‑M1/96‑3 (MQ=255)

                            aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgt                                                                <  1:8181975‑M1/96‑3 (MQ=255)

                             aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAGA‑AAATAACGTGATGCCGGGATGGTTTgta                                                               >  2:10680733‑M1/1‑93 (MQ=255)

                             aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgta                                                               <  1:6216890‑M1/96‑4 (MQ=255)

                             aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgta                                                               <  2:2207506‑M1/96‑4 (MQ=255)

                               tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatt                                                             >  1:8947811‑M1/1‑91 (MQ=255)

                                ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatat                                                            <  1:10174646‑M1/96‑7 (MQ=255)

                                 cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                           >  1:1873834‑M1/1‑89 (MQ=255)

                                 cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                           >  2:8685809‑M1/1‑89 (MQ=255)

                                 cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                           <  2:7803281‑M1/96‑8 (MQ=255)

                                  aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                          >  2:5694359‑M1/1‑88 (MQ=255)

                                  aCCGGAAAATTTACGGATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TCAAATAATGTGATGCCGGGATGGTTTgcatttcc                                                          <  2:13714868‑M1/96‑9 (MQ=255)

                                   ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAAACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                         >  1:2408500‑M1/1‑87 (MQ=255)

                                     ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                       <  1:4295259‑M1/96‑12 (MQ=255)

                                     ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccgg                                                       <  2:5383280‑M1/96‑12 (MQ=255)

                                        aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcgt                                                    <  1:14281097‑M1/96‑15 (MQ=255)

                                         aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGCTGGTTTgtattttccggcatc                                                   >  1:4833622‑M1/1‑81 (MQ=255)

                                         aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatc                                                   <  1:12825591‑M1/96‑16 (MQ=255)

                                          aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatct                                                  <  2:4730644‑M1/96‑17 (MQ=255)

                                            ttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                                >  2:859397‑M1/1‑78 (MQ=255)

                                             tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttca                                               <  1:9359390‑M1/96‑20 (MQ=255)

                                              aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                              <  1:10451678‑M1/96‑21 (MQ=255)

                                              aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                              <  1:13539082‑M1/96‑21 (MQ=255)

                                                 tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                           >  2:7220341‑M1/1‑73 (MQ=255)

                                                  ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                          >  1:11618278‑M1/1‑72 (MQ=255)

                                                  ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                          >  1:4004071‑M1/1‑72 (MQ=255)

                                                  ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                          >  2:2671758‑M1/1‑72 (MQ=255)

                                                   taGCAATAGAACACTTTGGGTGGTAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                         <  1:12542362‑M1/96‑26 (MQ=255)

                                                   taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcga                                         >  2:1382529‑M1/1‑71 (MQ=255)

                                                     gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                       >  2:10812995‑M1/1‑69 (MQ=255)

                                                     gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                       >  2:9733024‑M1/1‑69 (MQ=255)

                                                      cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                      >  1:7271394‑M1/1‑68 (MQ=255)

                                                      cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                      >  1:8677251‑M1/1‑68 (MQ=255)

                                                      cAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                      >  1:7932159‑M1/1‑68 (MQ=255)

                                                       aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                     <  1:3553187‑M1/96‑30 (MQ=255)

                                                       aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                     >  2:1910800‑M1/1‑67 (MQ=255)

                                                       aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                     <  2:6925284‑M1/96‑30 (MQ=255)

                                                        aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                    >  2:640091‑M1/1‑66 (MQ=255)

                                                         tggaaCACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                   >  1:10623145‑M1/3‑65 (MQ=255)

                                                          aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                  <  1:6936376‑M1/96‑33 (MQ=255)

                                                          aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                  >  1:3232342‑M1/1‑64 (MQ=255)

                                                          aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                  >  2:3928954‑M1/1‑64 (MQ=255)

                                                           gAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                                 <  2:9009241‑M1/96‑34 (MQ=255)

                                                                ccttggGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                            >  1:10823589‑M1/3‑58 (MQ=255)

                                                                cTTTGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                            >  1:6497748‑M1/1‑58 (MQ=255)

                                                                   tGGGTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactg                         >  2:170260‑M1/1‑55 (MQ=255)

                                                                    gggTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                        <  2:8301278‑M1/96‑43 (MQ=255)

                                                                    gggTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                        >  1:2584507‑M1/1‑54 (MQ=255)

                                                                    gggTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                        <  1:4560617‑M1/96‑43 (MQ=255)

                                                                     ggTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                       <  1:5500895‑M1/96‑44 (MQ=255)

                                                                      gTGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                      >  1:7417168‑M1/1‑52 (MQ=255)

                                                                       tGGAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg                     >  2:2304261‑M1/1‑51 (MQ=255)

                                                                        ggAAGAATAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgt                    <  1:955853‑M1/96‑47 (MQ=255)

                                                                        ggAAGAATAGACCGATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgt                    >  1:13884599‑M1/1‑50 (MQ=255)

                                                                         gaagaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                   >  2:7526248‑M1/1‑49 (MQ=255)

                                                                         gaagaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                   <  2:11337860‑M1/96‑48 (MQ=255)

                                                                         gaagaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                   >  2:9381164‑M1/1‑49 (MQ=255)

                                                                         gaagaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                   <  2:3743224‑M1/96‑48 (MQ=255)

                                                                         gaagaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                   <  2:6761172‑M1/96‑48 (MQ=255)

                                                                          aagaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcg                  <  1:1208545‑M1/96‑49 (MQ=255)

                                                                            gaaTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcg                <  2:5143458‑M1/96‑51 (MQ=255)

                                                                              aTAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac              >  2:10180185‑M1/1‑44 (MQ=255)

                                                                               tAGACCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc             >  2:4197656‑M1/1‑43 (MQ=255)

                                                                                  aCCTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc          <  1:8504047‑M1/96‑57 (MQ=255)

                                                                                   ccTATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatca         <  2:5917895‑M1/96‑58 (MQ=255)

                                                                                     tATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacc       <  2:4095044‑M1/96‑60 (MQ=255)

                                                                                     tATCACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacc       <  2:14161323‑M1/96‑60 (MQ=255)

                                                                                        cACTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct    <  1:1577048‑M1/96‑63 (MQ=255)

                                                                                          cTGCA‑TAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttaaagcgatagcaattatttactgagcgtcgcgaccatcaccgcttt  <  2:4938345‑M1/96‑65 (MQ=255)

ATATTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCA‑T‑AAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGC  >  NC_000913/4474728‑4474912

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