BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,073,511	A→G	L231P (CTG→CCG)	yggF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,073,511	0	A	G	99.4%	552.4 / ‑4.7	173	L231P (CTG→CCG)	yggF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (106/66); total (107/66)
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.03e-01

Only 100 of 184 total aligned reads displayed.

TCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATC  >  NC_000913/3073421‑3073599

tCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGc                                                                                        >  1:2972126/1‑93 (MQ=255)

 cAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCt                                                                                       >  1:3162471/1‑93 (MQ=255)

 cAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCt                                                                                       >  2:1974500/1‑93 (MQ=255)

  aaCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTATCCTGGCAAAAGTCAATCGGCTc                                                                                      <  2:442584/93‑1 (MQ=255)

  aaCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTc                                                                                      <  1:3863936/93‑1 (MQ=255)

  aaCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTc                                                                                      <  1:3265383/93‑1 (MQ=255)

  aaCCCCCACTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTGTAGCCTGGCAAAAGGCAATCGGCTc                                                                                      <  1:3863937/93‑1 (MQ=255)

   aCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCg                                                                                     >  2:2764902/1‑93 (MQ=255)

   aCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCg                                                                                     >  2:787214/1‑93 (MQ=255)

   aCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCg                                                                                     <  1:2365787/93‑1 (MQ=255)

      cccATTGCCTTGCAACGCTTGCGCTCCTGCTCTTCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGcc                                                                                  >  2:3655987/1‑93 (MQ=255)

      cccATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGcc                                                                                  <  1:2917142/93‑1 (MQ=255)

         aTTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGc                                                                               <  2:2436857/93‑1 (MQ=255)

          ttGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCa                                                                              >  1:21420/1‑93 (MQ=255)

          ttGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCa                                                                              >  1:3037890/1‑93 (MQ=255)

          ttGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCa                                                                              >  2:604697/1‑93 (MQ=255)

           tGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCat                                                                             <  1:583303/93‑1 (MQ=255)

            gCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCata                                                                            <  2:436115/93‑1 (MQ=255)

             ccTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTATCCTGGCAAAAGTCAATCGGCTCGGCCTGCatat                                                                           <  2:1553580/93‑1 (MQ=255)

             ccTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCatat                                                                           >  2:3141535/1‑93 (MQ=255)

              cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATc                                                                          >  1:2123667/1‑93 (MQ=255)

              cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATc                                                                          >  2:737625/1‑93 (MQ=255)

              cTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATc                                                                          <  2:1361036/93‑1 (MQ=255)

               ttGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTAAATCGGCTCAGCCTGCATATcc                                                                         >  2:876385/1‑93 (MQ=255)

                   aaCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGc                                                                     >  1:2774162/1‑93 (MQ=255)

                   aaCGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGc                                                                     <  1:2837141/93‑1 (MQ=255)

                     cGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCACAAGTCAATCGGCTCGGCCTGCCTATCCCCGccc                                                                   >  2:1744788/1‑92 (MQ=255)

                     cGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCt                                                                   <  2:2587578/93‑1 (MQ=255)

                        ttGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAAt                                                                <  2:3104391/93‑1 (MQ=255)

                        ttGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAAt                                                                <  2:3923213/93‑1 (MQ=255)

                         tGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATg                                                               >  1:4323986/1‑93 (MQ=255)

                           cgcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTTCATATCCCCGCCTAATGCt                                                             >  1:3647038/1‑93 (MQ=255)

                            gcTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCtt                                                            <  1:2612413/93‑1 (MQ=255)

                             cTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCttt                                                           >  1:693982/1‑93 (MQ=255)

                              tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTa                                                          >  2:352979/1‑93 (MQ=255)

                              tCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTa                                                          >  1:2719936/1‑93 (MQ=255)

                                cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAAc                                                        >  1:3172492/1‑93 (MQ=255)

                                cTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATACCCGCCTAATGCTTTAAc                                                        <  2:2706926/93‑1 (MQ=255)

                                 tGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACg                                                       <  2:2358006/93‑1 (MQ=255)

                                 tGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACg                                                       <  1:4065725/93‑1 (MQ=255)

                                 tGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACg                                                       <  1:115560/93‑1 (MQ=255)

                                  gCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACgg                                                      >  2:3736124/1‑93 (MQ=255)

                                    tctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGca                                                    >  1:1392335/1‑93 (MQ=255)

                                    tctGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGca                                                    >  1:4698918/1‑93 (MQ=255)

                                      tGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATACCCGCCTAATGCTTTAACGGcaca                                                  >  1:3409245/1‑93 (MQ=255)

                                       gCAATCTGCCTATTTGCCGTGTAATCCCCTTTAGCTTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                 <  2:3918515/93‑1 (MQ=255)

                                       gCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                 >  1:461894/1‑93 (MQ=255)

                                       ccAACCTGCCGATTTTCCGTGTAATCCCCTTTACCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAg                                                 <  2:2712543/92‑1 (MQ=255)

                                         aaTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGc                                               >  2:901008/1‑93 (MQ=255)

                                         aaTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGc                                               >  1:110082/1‑93 (MQ=255)

                                          aTCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCa                                              >  2:4253201/1‑93 (MQ=255)

                                           tCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAg                                             >  2:3358877/1‑93 (MQ=255)

                                            cTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGa                                            >  1:3429971/1‑93 (MQ=255)

                                            cTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACCGCACAGGCAGa                                            >  1:1488588/1‑93 (MQ=255)

                                             tGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGGATATCCCCGCCTAATGCTTTCACGGCACAGGCAGaa                                           >  1:3380340/1‑93 (MQ=255)

                                             tGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaa                                           >  1:1204985/1‑93 (MQ=255)

                                              gCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTTACGGCACAGGCAGaaa                                          >  2:3959158/1‑93 (MQ=255)

                                              gCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaaa                                          <  1:426344/93‑1 (MQ=255)

                                              gCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGaaa                                          <  2:401863/93‑1 (MQ=255)

                                                cGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATc                                        >  2:1765712/1‑93 (MQ=255)

                                                cGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATc                                        >  2:4652942/1‑93 (MQ=255)

                                                 gATTTTCCGTGTAATCCCCTTTATCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCa                                       <  1:1374365/93‑1 (MQ=255)

                                                 gATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAAtgc                                       >  1:3599722/1‑91 (MQ=255)

                                                 gATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCa                                       >  1:2742798/1‑93 (MQ=255)

                                                  aTTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCAc                                      <  2:2423615/93‑1 (MQ=255)

                                                   ttttCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACg                                     >  1:1351324/1‑93 (MQ=255)

                                                   ttttCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACg                                     >  1:952211/1‑93 (MQ=255)

                                                        cGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTc                                <  1:4650652/93‑1 (MQ=255)

                                                         gtgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCt                               >  1:848592/1‑93 (MQ=255)

                                                         gtgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCt                               >  2:434600/1‑93 (MQ=255)

                                                          tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACCGGCAGAAATCACGCCTTCTg                              >  1:1892242/1‑93 (MQ=255)

                                                          tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTg                              >  2:544367/1‑93 (MQ=255)

                                                          tgtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTg                              >  1:4552696/1‑93 (MQ=255)

                                                           gtAATCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTgg                             <  2:1885816/93‑1 (MQ=255)

                                                             aaTCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGgcg                           >  1:1025879/1‑93 (MQ=255)

                                                               tCCCCTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGgcgcc                         >  1:4040905/1‑92 (MQ=255)

                                                                ccccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAc                        >  2:3484472/1‑93 (MQ=255)

                                                                ccccTTTAGCCTGGAAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATAACACCTTCTGGCGCAc                        >  2:3808419/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                       >  1:868064/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                       >  2:236375/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                       >  1:3578606/1‑93 (MQ=255)

                                                                 cccTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAcc                       <  1:3026904/93‑1 (MQ=255)

                                                                   cTTTATCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgc                     <  2:4132457/93‑1 (MQ=255)

                                                                   cTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgc                     >  1:1984050/1‑93 (MQ=255)

                                                                   cTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgc                     >  2:1695325/1‑93 (MQ=255)

                                                                   cTTTAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgc                     >  2:4458947/1‑93 (MQ=255)

                                                                     ttAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgccg                   <  1:3075483/93‑1 (MQ=255)

                                                                     ttAGCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCAccgccg                   >  1:2888729/1‑93 (MQ=255)

                                                                        gCCTGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATg                >  2:4580336/1‑93 (MQ=255)

                                                                           tGGCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtg             >  2:1013403/1‑93 (MQ=255)

                                                                            ggCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGgtgt            >  1:2022955/1‑93 (MQ=255)

                                                                             gCAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTa           <  2:2066779/93‑1 (MQ=255)

                                                                              cAAAAGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTac          <  1:4566554/93‑1 (MQ=255)

                                                                               aaaaGTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTaca         >  2:3041183/1‑93 (MQ=255)

                                                                                   gTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcac     >  2:2556031/1‑93 (MQ=255)

                                                                                   gTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcac     <  1:1951651/93‑1 (MQ=255)

                                                                                   gTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcac     <  1:2202888/93‑1 (MQ=255)

                                                                                   gTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcac     <  2:717411/93‑1 (MQ=255)

                                                                                   gTCAATCGGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcac     >  2:3241687/1‑93 (MQ=255)

                                                                                      aaTCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGAACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTacatcacatc  >  1:1638585/1‑93 (MQ=255)

TCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATC  >  NC_000913/3073421‑3073599

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