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%	890.4 / ‑1.2	250	coding (1849/1851 nt)	yjgL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base A (147/101); total (149/101)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 270 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAA‑TAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

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

    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaa                                                                                      >  2:5191840/1‑92 (MQ=255)

    tGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaa                                                                                      >  1:561779/1‑92 (MQ=255)

       aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat‑aaaataa                                                                                   >  1:6081954/1‑89 (MQ=255)

         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtg                                                                                 >  1:5028921/1‑93 (MQ=255)

         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtg                                                                                 >  1:7266094/1‑93 (MQ=255)

         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtg                                                                                 >  2:3683632/1‑93 (MQ=255)

         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtg                                                                                 >  1:2723365/1‑93 (MQ=255)

          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtgt                                                                                >  2:7562592/1‑93 (MQ=255)

          gCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAAtgt                                                                                >  2:3804280/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGGCCTATCACTGCATAAAA‑TAATGTGa                                                                              >  1:310965/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGa                                                                              >  1:6982512/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGa                                                                              >  1:8378712/1‑93 (MQ=255)

            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGa                                                                              >  1:2446873/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                                             >  1:3081499/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                                             >  1:3316890/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                                             >  2:5925318/1‑93 (MQ=255)

             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGAt                                                                             >  2:1514259/1‑93 (MQ=255)

              ttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATg                                                                            >  1:3905083/1‑93 (MQ=255)

                 aataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCg                                                                         <  1:8509184/93‑1 (MQ=255)

                   taaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCggg                                                                       >  1:5726922/1‑93 (MQ=255)

                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGAt                                                                     >  2:7304982/1‑93 (MQ=255)

                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATg                                                                    <  1:7623971/93‑1 (MQ=255)

                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATg                                                                    <  2:912495/93‑1 (MQ=255)

                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATg                                                                    >  2:320808/1‑93 (MQ=255)

                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATg                                                                    >  2:7902470/1‑93 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGt                                                                  <  2:1168458/93‑1 (MQ=255)

                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGt                                                                  >  2:5456873/1‑93 (MQ=255)

                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATggct                                                                 >  2:4244897/1‑91 (MQ=255)

                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGAGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGtt                                                                 <  1:7415048/93‑1 (MQ=255)

                          aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGttt                                                                <  2:3641214/93‑1 (MQ=255)

                          aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGttt                                                                >  1:7049201/1‑93 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgt                                                              >  1:463225‑M1/1‑91 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgt                                                              <  1:1872316‑M1/93‑3 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgt                                                              >  1:5937931‑M1/1‑91 (MQ=255)

                            tCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgt                                                              <  1:6213437‑M1/93‑3 (MQ=255)

                             ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgta                                                             <  1:1081891‑M1/93‑4 (MQ=255)

                               aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatt                                                           <  1:7615878‑M1/93‑6 (MQ=255)

                                ccGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtattt                                                          >  1:128046‑M1/1‑87 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  2:8537552‑M1/93‑10 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  1:2134528‑M1/1‑84 (MQ=255)

                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  2:6573747‑M1/93‑10 (MQ=255)

                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccgg                                                     >  1:3052342‑M1/1‑82 (MQ=255)

                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    <  1:1043195‑M1/93‑13 (MQ=255)

                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  1:5914756‑M1/1‑76 (MQ=255)

                                            cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              <  2:484080‑M1/93‑19 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  2:612705‑M1/1‑74 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  1:7697874‑M1/1‑74 (MQ=255)

                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  2:6673565‑M1/93‑20 (MQ=255)

                                              tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            <  1:1474681‑M1/93‑21 (MQ=255)

                                              tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            >  2:4790586‑M1/1‑73 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  1:6151364‑M1/93‑22 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  1:412526‑M1/93‑22 (MQ=255)

                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  2:189742‑M1/93‑22 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:4541887‑M1/1‑66 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:3975209‑M1/1‑66 (MQ=255)

                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  1:3464661‑M1/1‑66 (MQ=255)

                                                       aGAACACTTTGGGTTGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttctcggcatctttatagcgataga                                   >  1:2353954‑M1/1‑64 (MQ=255)

                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                   >  1:1440304‑M1/1‑64 (MQ=255)

                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                   <  1:8695382‑M1/93‑30 (MQ=255)

                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgataga                                   >  2:3066383‑M1/1‑64 (MQ=255)

                                                          acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                <  1:4899588‑M1/93‑33 (MQ=255)

                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               >  1:812277‑M1/1‑60 (MQ=255)

                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               <  1:5986603‑M1/93‑34 (MQ=255)

                                                            acTTTGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              >  2:161425‑M1/1‑59 (MQ=255)

                                                                tGGGTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                          >  2:2302490‑M1/1‑55 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:549303‑M1/93‑40 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         >  2:560905‑M1/1‑54 (MQ=255)

                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         >  2:3465692‑M1/1‑54 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      >  2:4395920‑M1/1‑51 (MQ=255)

                                                                    tGGAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                      <  1:6891861‑M1/93‑43 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  1:357281‑M1/1‑50 (MQ=255)

                                                                     ggAAGAATAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     <  1:6381783‑M1/93‑44 (MQ=255)

                                                                      gaagaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                    >  1:2923026‑M1/1‑49 (MQ=255)

                                                                       aagaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg                   <  1:6255840‑M1/93‑46 (MQ=255)

                                                                       aagaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg                   >  2:6886954‑M1/1‑48 (MQ=255)

                                                                        agaaTAGCCCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtattttccggcatctttatagcgatagcaattatttactgagcgt                  >  2:1358233‑M1/1‑47 (MQ=255)

                                                                        agaaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgt                  <  1:7191696‑M1/93‑47 (MQ=255)

                                                                          aaTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcg                <  2:8732241‑M1/93‑49 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               >  1:1638666‑M1/1‑44 (MQ=255)

                                                                           aTAGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc               <  1:2217612‑M1/93‑50 (MQ=255)

                                                                             aGACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcga             >  2:1037691‑M1/1‑42 (MQ=255)

                                                                              gACCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            >  2:2829178‑M1/1‑41 (MQ=255)

                                                                               aCCTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc           <  1:5299035‑M1/93‑54 (MQ=255)

                                                                                ccTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          >  1:3902928‑M1/1‑39 (MQ=255)

                                                                                ccTATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          <  1:8300795‑M1/93‑55 (MQ=255)

                                                                                 cTATCACTGCAT‑AAAGTAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccat         >  1:8587707‑M1/1‑38 (MQ=255)

                                                                                  tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc        >  1:4437231‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc        >  2:4088806‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc        >  2:1089133‑M1/1‑37 (MQ=255)

                                                                                  tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc        <  2:2022014‑M1/93‑57 (MQ=255)

                                                                                  tATCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtagcgaccatc        >  2:828600‑M1/1‑37 (MQ=255)

                                                                                    tCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac      <  1:4389366‑M1/93‑59 (MQ=255)

                                                                                    tCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac      <  2:390581‑M1/93‑59 (MQ=255)

                                                                                    tCACTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac      >  2:2572746‑M1/1‑35 (MQ=255)

                                                                                      aCTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccg    >  1:7970798‑M1/1‑33 (MQ=255)

                                                                                       cTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   <  2:6762892‑M1/93‑62 (MQ=255)

                                                                                       cTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   <  2:8165142‑M1/93‑62 (MQ=255)

                                                                                       cTGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   >  1:6517463‑M1/1‑32 (MQ=255)

                                                                                        tGCATAAAA‑TAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  <  1:3760236‑M1/93‑63 (MQ=255)

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCAT‑AAA‑TAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

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