BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,095,111	G→A	intergenic (‑82/‑206)	rhaB ← / → rhaS	rhamnulokinase/RhaS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,095,111	0	G	A	99.8%	2359.6 / ‑6.2	643	intergenic (‑82/‑206)	rhaB/rhaS	rhamnulokinase/RhaS
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (420/221); total (421/222)
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.39e-01

Only 100 of 648 total aligned reads displayed.

AAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095022‑4095200

gtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGcc                                                                                        >  1:7905781‑M2/12‑93 (MQ=255)

 tcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCa                                                                                       <  2:2803354‑M2/83‑1 (MQ=255)

  cgacatatgTTGATCCTGCTGAATGTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAt                                                                                      <  2:6592434‑M2/84‑1 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     >  2:311639‑M2/9‑93 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     <  1:9687701‑M2/85‑1 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     >  2:4234073‑M2/9‑93 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCAtt                                                                                     <  2:206882‑M2/85‑1 (MQ=255)

   gacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGACAtt                                                                                     >  1:5521758‑M2/9‑93 (MQ=255)

     catatgGTGATCCTGCTGAAGTTCATTACGACCAGTCTAAAAATCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTgg                                                                                   <  2:4989098‑M2/87‑1 (MQ=255)

       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  1:524120‑M2/5‑93 (MQ=255)

       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  1:5124755‑M2/5‑93 (MQ=255)

       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  2:6080200‑M2/5‑93 (MQ=255)

       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 >  2:5317437‑M2/5‑93 (MQ=255)

       tatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCa                                                                                 <  1:1382961‑M2/89‑1 (MQ=255)

        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCaa                                                                                <  1:3158158‑M2/90‑1 (MQ=255)

        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCaa                                                                                >  1:7223903‑M2/4‑93 (MQ=255)

        atgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCaa                                                                                >  2:4448292‑M2/4‑93 (MQ=255)

         tgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAAc                                                                               >  1:4331194‑M2/3‑93 (MQ=255)

              aTCGTGCTGAATTTCATTACGACCACTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAggg                                                                          <  2:7464017/93‑1 (MQ=255)

              aTCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAggg                                                                          >  1:8581942/1‑93 (MQ=255)

               tCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGa                                                                         >  2:9037578/1‑93 (MQ=255)

                ccTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaa                                                                        <  2:6975953/93‑1 (MQ=255)

                ccTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaa                                                                        >  1:6650601/1‑93 (MQ=255)

                ccTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCGCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGaa                                                                        <  2:7988340/93‑1 (MQ=255)

                  tgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAg                                                                      <  1:5857883/93‑1 (MQ=255)

                  tgctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAg                                                                      >  2:1884616/1‑93 (MQ=255)

                    ctgAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGAt                                                                    >  2:2964633/1‑93 (MQ=255)

                      gAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGa                                                                  >  1:3035309/1‑93 (MQ=255)

                      gAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGa                                                                  >  2:8611140/1‑93 (MQ=255)

                         tttCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACg                                                               >  2:5466500/1‑93 (MQ=255)

                           tCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtg                                                             <  1:1073331/93‑1 (MQ=255)

                             aTTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGTAACCAGGGAAAGATGAACGtgat                                                           >  2:7209338/1‑93 (MQ=255)

                             aTTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgat                                                           <  2:1130969/93‑1 (MQ=255)

                              ttACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgatg                                                          >  2:3622469/1‑93 (MQ=255)

                               tACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgatga                                                         <  1:6362461/93‑1 (MQ=255)

                               tACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGtgatga                                                         >  2:6636481/1‑93 (MQ=255)

                                   aCCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGtt                                                     >  1:4497172/1‑93 (MQ=255)

                                   aCCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGtt                                                     >  1:4229719/1‑93 (MQ=255)

                                   aCCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGtt                                                     >  1:9797769/1‑93 (MQ=255)

                                     cAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTca                                                   <  1:1754899/93‑1 (MQ=255)

                                        tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  2:5958832/1‑93 (MQ=255)

                                        tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                <  2:3703065/93‑1 (MQ=255)

                                        tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  1:811279/1‑93 (MQ=255)

                                        tCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACaa                                                >  2:2254202/1‑93 (MQ=255)

                                         cTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAt                                               >  2:63744/1‑93 (MQ=255)

                                         cTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAt                                               >  1:1290940/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              >  1:4088015/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              >  1:1448254/1‑93 (MQ=255)

                                          tAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAtt                                              <  1:6128959/93‑1 (MQ=255)

                                           aaaaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAttt                                             >  1:8223111/1‑93 (MQ=255)

                                            aaaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTg                                            >  2:885190/1‑93 (MQ=255)

                                            aaaaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTg                                            >  1:3251784/1‑93 (MQ=255)

                                              aaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCt                                          <  1:1128542/93‑1 (MQ=255)

                                              aaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCt                                          <  1:6101000/93‑1 (MQ=255)

                                              aaGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCt                                          <  1:9081617/93‑1 (MQ=255)

                                               aGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTg                                         <  2:825365/93‑1 (MQ=255)

                                                gcgcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGa                                        <  2:1000508/93‑1 (MQ=255)

                                                gcgcCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGa                                        >  2:8429673/1‑93 (MQ=255)

                                                    cTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtg                                    <  2:5646482/93‑1 (MQ=255)

                                                     tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   <  2:7795102/93‑1 (MQ=255)

                                                     tGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATtgt                                   >  2:2178000/1‑93 (MQ=255)

                                                       aaTTCGCGACCTTCTCGTTACTGACAGGAAAATTAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtgg                                 <  2:4829630/93‑1 (MQ=255)

                                                       aaTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtgg                                 <  1:249577/93‑1 (MQ=255)

                                                        aTTCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggt                                >  2:1684314/1‑93 (MQ=255)

                                                         ttCGCGACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggtg                               >  1:7399767/1‑93 (MQ=255)

                                                           cgcgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgat                             <  2:943909/93‑1 (MQ=255)

                                                            gcgACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatg                            >  1:6509191/1‑93 (MQ=255)

                                                              gACCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtg                          <  1:4978598/93‑1 (MQ=255)

                                                               aCCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtga                         <  2:2439370/93‑1 (MQ=255)

                                                               aCCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtga                         <  1:4019018/93‑1 (MQ=255)

                                                               aCCTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtga                         >  2:4669436/1‑93 (MQ=255)

                                                                ccTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtgat                        <  2:6217550/93‑1 (MQ=255)

                                                                 cTTCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGtgatgtgatg                       >  1:2313104/1‑93 (MQ=255)

                                                                  ttCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGc                      >  2:9843257/1‑93 (MQ=255)

                                                                  ttCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGc                      >  2:146471/1‑93 (MQ=255)

                                                                  ttCTCGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGc                      >  1:8999902/1‑93 (MQ=255)

                                                                   tctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCt                     >  1:3531836/1‑93 (MQ=255)

                                                                   tctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCt                     >  1:537909/1‑93 (MQ=255)

                                                                    ctcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTc                    >  1:2667372/1‑93 (MQ=255)

                                                                     tcGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCa                   >  2:5173144/1‑93 (MQ=255)

                                                                      cGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAc                  <  1:2902232/93‑1 (MQ=255)

                                                                      cGTTACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCAc                  <  1:9276049/93‑1 (MQ=255)

                                                                        ttACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCg                >  1:4625117/1‑93 (MQ=255)

                                                                        ttACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCg                >  1:950617/1‑93 (MQ=255)

                                                                         tACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGc               >  2:9636438/1‑93 (MQ=255)

                                                                         tACTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGc               >  1:6245357/1‑93 (MQ=255)

                                                                          aCTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCa              >  1:8822719/1‑93 (MQ=255)

                                                                          aCTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCa              >  2:2552788/1‑93 (MQ=255)

                                                                           cTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAt             >  2:3398613/1‑93 (MQ=255)

                                                                           cTGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAt             <  1:2200144/93‑1 (MQ=255)

                                                                            tGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAtt            >  1:7590617/1‑93 (MQ=255)

                                                                            tGACAGGAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAtt            >  2:9281243/1‑93 (MQ=255)

                                                                                 ggAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTg       >  2:1265891/1‑93 (MQ=255)

                                                                                 ggAAAATGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTg       >  1:6977783/1‑93 (MQ=255)

                                                                                   aaaaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaa     >  1:160862/1‑93 (MQ=255)

                                                                                   aaaaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaa     <  2:4131169/93‑1 (MQ=255)

                                                                                    aaaTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGaaa    <  2:7615835/93‑1 (MQ=255)

                                                                                      aTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAt  >  1:5657021/1‑93 (MQ=255)

                                                                                      aTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAt  >  1:3715279/1‑93 (MQ=255)

                                                                                      aTGAGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAt  <  1:7748307/93‑1 (MQ=255)

AAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAAT  >  NC_000913/4095022‑4095200

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: `‑`