BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,420	A→C	5.4%	L214L (CTT→CTG) ‡	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,420	0	A	C	5.4%	334.9 / 8.3	112	L214L (CTT→CTG) ‡	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base A (57/49); new base C (2/4); total (59/53)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.55e-01

Only 100 of 117 total aligned reads displayed.

GTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAG  >  NC_000913/858271‑858569

gTCCAGCCCTGGAATGTCCGAAACACGGACCATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGa                                                                                                                                                       >  1:139506/1‑150 (MQ=255)

 accAGCCCTGGAATGTCGGAAACACGGAACATATGGTGGTGTTGTCCACCGTGTACTCTTCAACATCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACATTGGCGAAAACGCGCTGATACGCTGGaa                                                                                                                                                      <  2:1201078/149‑1 (MQ=255)

   cAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGaatc                                                                                                                                                    >  2:318850/1‑148 (MQ=255)

   cAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagc                                                                                                                                                    >  1:1095829/1‑150 (MQ=255)

   cAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagc                                                                                                                                                    >  2:1155692/1‑150 (MQ=255)

    aGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagca                                                                                                                                                   <  1:206922/150‑1 (MQ=255)

    aGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagca                                                                                                                                                   >  2:132748/1‑150 (MQ=255)

     gCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcag                                                                                                                                                  <  2:523338/150‑1 (MQ=255)

      cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                 >  2:20712/1‑150 (MQ=255)

      cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                 >  2:385622/1‑150 (MQ=255)

      cccTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAagcagc                                                                                                                                                 >  2:465984/1‑150 (MQ=255)

              tGTCCGAAACACTTAACATATGGTCCTGTTGTCAACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGCCTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGAAACGCTGGAAGCAGCCAGAGTTc                                                                                                                                         >  2:1614086/1‑150 (MQ=255)

               gTCCTAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                        <  2:1175634/150‑1 (MQ=255)

               gTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTcc                                                                                                                                        <  2:1559153/150‑1 (MQ=255)

                 ccGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAg                                                                                                                                      <  1:276501/150‑1 (MQ=255)

                 ccGAAACACGGAACATTTGGTGGTGTTGCCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAg                                                                                                                                      <  1:1331698/150‑1 (MQ=255)

                  cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGt                                                                                                                                     >  2:759434/1‑150 (MQ=255)

                  cGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagc                                                                                                                                     <  2:1598999/150‑7 (MQ=255)

                   gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACTTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                    <  2:1420950/150‑8 (MQ=255)

                   gAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGCTGCTTCCAGCGtatcagcg                                                                                                                                    <  2:226866/150‑8 (MQ=255)

                       cacGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGcccc                                                                                                                                >  1:1485730/1‑150 (MQ=255)

                        ccGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGAGCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCg                                                                                                                               <  1:1168514/149‑1 (MQ=255)

                           gAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAg                                                                                                                            >  1:1469955/1‑150 (MQ=255)

                              cATTTGGTGGTGTTGTCCACAGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCg                                                                                                                         >  1:357448/1‑150 (MQ=255)

                               aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                        >  1:994819/1‑150 (MQ=255)

                               aTTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCgg                                                                                                                        >  1:1183950/1‑150 (MQ=255)

                                tttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGt                                                                                                                       <  1:59591/150‑1 (MQ=255)

                                 ttGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTa                                                                                                                      <  1:1138310/150‑1 (MQ=255)

                                 ttGGTGGTGTTGTCCACCGTGTACTCTTAAACTTCCGTACGATGTGCCGCATGACAGGGATCATTTTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGATGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCTGAGTCGGTa                                                                                                                      <  1:1261654/150‑1 (MQ=255)

                                  tggtggTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTAt                                                                                                                     >  1:1067686/1‑150 (MQ=255)

                                   ggtggtGTTGTCGACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATg                                                                                                                    <  1:1137841/150‑1 (MQ=255)

                                     tggtgTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATgcg                                                                                                                  >  2:1027619/1‑150 (MQ=255)

                                       gtgtTGTCCACCGTGGACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGGTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCt                                                                                                                <  1:1172848/150‑1 (MQ=255)

                                           tgtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCGAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCaa                                                                                                            <  2:642631/150‑1 (MQ=255)

                                           tgtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCaa                                                                                                            >  2:1244306/1‑150 (MQ=255)

                                            gtCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAg                                                                                                           >  1:454049/1‑150 (MQ=255)

                                               cACCGTGTACTGTTCAACTTCCGTACGATTTGCCGCATGCCAGGGATCATATTGAGCCAGATTGCCATTAAAGACTTTGGCGAAAACGCGCTGATTCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                        <  1:1737672/150‑1 (MQ=255)

                                               cACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAcc                                                                                                        <  2:207784/150‑1 (MQ=255)

                                                 ccGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCtt                                                                                                      >  2:1499304/1‑150 (MQ=255)

                                                    tgtACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTgg                                                                                                   <  1:385622/150‑1 (MQ=255)

                                                     gtACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGa                                                                                                  <  1:486530/150‑1 (MQ=255)

                                                      tACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAg                                                                                                 <  1:1155692/150‑1 (MQ=255)

                                                         tctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTc                                                                                              <  2:1339213/150‑1 (MQ=255)

                                                          ctTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCg                                                                                             >  2:1257859/1‑150 (MQ=255)

                                                           ttCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTACGGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGt                                                                                            <  1:924529/150‑1 (MQ=255)

                                                                cTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCTTTCCgg                                                                                       >  1:166479/1‑150 (MQ=255)

                                                                 ttCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCggg                                                                                      <  1:1551812/150‑1 (MQ=255)

                                                                 ttCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCggg                                                                                      >  1:1734094/1‑150 (MQ=255)

                                                                     gTACGATGTGCCGCTTGGGAGGTATCATATTGCGCCATATTGCCATTAAAGTCGTTGGCGAAATCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGAGCTTTGGAGGTCGTTCCGGGCGGa                                                                                  <  2:1558399/150‑1 (MQ=255)

                                                                       aCGATGTGCCGCCTGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcg                                                                                <  2:1701947/150‑1 (MQ=255)

                                                                         gATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGGCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGAcggc                                                                              <  2:1095829/150‑1 (MQ=255)

                                                                             tgCCGCATGCCAGGGATCATATTGCGCCAGATTGACATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGAcggcatcg                                                                          <  2:907376/150‑5 (MQ=255)

                                                                             ggCCGCACGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGAcggcggcg                                                                          <  2:320507/149‑1 (MQ=255)

                                                                               ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGATGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                        >  1:945627/1‑150 (MQ=255)

                                                                               ccGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGa                                                                        >  1:799307/1‑150 (MQ=255)

                                                                                   aTGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACg                                                                    <  2:1433659/150‑1 (MQ=255)

                                                                                      ccAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTc                                                                 >  2:1553644/1‑150 (MQ=255)

                                                                                         gggATCATATTGCGCCAGATTGCCCTTAAAGACGTTGGCGAACACGCGCTGATACGCTGGAAGCAGCCAGCGGTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAgg                                                              <  1:234037/150‑1 (MQ=255)

                                                                                          ggATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAggg                                                             <  2:357448/150‑1 (MQ=255)

                                                                                            aTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCTGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTa                                                           <  2:933673/150‑1 (MQ=255)

                                                                                            aTCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTa                                                           >  2:1530549/1‑150 (MQ=255)

                                                                                               atatTGCGACAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGCAAGAt                                                        >  2:1312742/1‑150 (MQ=255)

                                                                                                tatTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATg                                                       >  1:472634/1‑150 (MQ=255)

                                                                                                 atTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGa                                                      >  2:1251704/1‑150 (MQ=255)

                                                                                                 atTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGa                                                      >  2:1480581/1‑150 (MQ=255)

                                                                                                     cgcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCt                                                  >  2:422922/1‑150 (MQ=255)

                                                                                                      gcCAGATTGCCATTAAAGACGTTGGCGACAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTc                                                 <  2:241148/150‑1 (MQ=255)

                                                                                                      gcCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCTGACGGCGGCGGATACGGTCAGGGTAGATGACGCTc                                                 >  1:1626453/1‑150 (MQ=255)

                                                                                                       ccAGATTGCCATTAACGACGTGGGCGAAAACGCGCGGATACGCTGGAAGCAGCCAGCGTTCCAGGGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCGca                                                <  1:1146456/150‑1 (MQ=255)

                                                                                                       ccAGATTGCCATTAAAGCCGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTca                                                <  1:1396123/150‑1 (MQ=255)

                                                                                                       ccAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTca                                                >  1:1769110/1‑150 (MQ=255)

                                                                                                        aaGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCGCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGGCAGGGTAGATGACGCTcac                                               <  2:865694/149‑1 (MQ=255)

                                                                                                            ttGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGCTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATcc                                           >  1:464046/1‑150 (MQ=255)

                                                                                                               ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                        <  1:427610/150‑1 (MQ=255)

                                                                                                               ccATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCggg                                        >  2:609628/1‑150 (MQ=255)

                                                                                                                cATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGt                                       >  2:1035517/1‑150 (MQ=255)

                                                                                                                 aTTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGtt                                      >  1:1222405/1‑150 (MQ=255)

                                                                                                                  ttAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGGCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTa                                     <  1:759434/150‑1 (MQ=255)

                                                                                                                    aaaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGATaaa                                   >  1:697145/1‑150 (MQ=255)

                                                                                                                     aaGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAAc                                  >  2:1380382/1‑150 (MQ=255)

                                                                                                                          gTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCCTTCCAGTGCCCCGGAGTCGGTATGCACTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTg                             >  2:795418/1‑150 (MQ=255)

                                                                                                                           ttGGCGAAAACGCGCTGATACGCGGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGc                            <  1:1678478/150‑1 (MQ=255)

                                                                                                                              gCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCttt                         >  1:823018/1‑150 (MQ=255)

                                                                                                                               cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                        >  1:1114659/1‑150 (MQ=255)

                                                                                                                               cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                        <  2:1183950/150‑1 (MQ=255)

                                                                                                                               cGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGCACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTc                        >  2:1077907/1‑150 (MQ=255)

                                                                                                                                 aaaaCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCa                      >  1:1761730/1‑150 (MQ=255)

                                                                                                                                  aaaCGCGCAGATACGCTGGAAGCAGCCAGCGGTCCAGGGCCCCGGCGGCGGGATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCAt                     <  1:318850/150‑1 (MQ=255)

                                                                                                                                    aCGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCa                   <  1:1416127/150‑1 (MQ=255)

                                                                                                                                     cgcgCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCAc                  <  2:1439246/150‑1 (MQ=255)

                                                                                                                                       cgcTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACtt                <  2:594855/150‑1 (MQ=255)

                                                                                                                                          tggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCa             >  1:254104/5‑150 (MQ=255)

                                                                                                                                           ggaaCGCTGGCTGCTTCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCaa            >  1:1296014/4‑150 (MQ=255)

                                                                                                                                            aTACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGGTCCGGGCGGCCGGCGGCGGAAACGGTCAGGGTAGATGACGCTCACATCCGGGTGACACCTTTGCTTTCCATCACTTTCaaa           >  1:709491/1‑150 (MQ=255)

                                                                                                                                              aCGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATg         <  1:1212155/150‑1 (MQ=255)

                                                                                                                                                gCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTc       >  1:1709710/1‑150 (MQ=255)

                                                                                                                                                 cTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTcc      <  2:1818472/150‑1 (MQ=255)

                                                                                                                                                   ggAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCcac    >  1:1712902/1‑150 (MQ=255)

                                                                                                                                                   ggAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCAATCACTTTCAAATGTcccc    >  2:221875/1‑148 (MQ=255)

                                                                                                                                                     aaGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAg  <  2:1769110/150‑1 (MQ=255)

GTCCAGCCCTGGAATGTCCGAAACACGGAACATTTGGTGGTGTTGTCCACCGTGTACTCTTCAACTTCCGTACGATGTGCCGCATGCCAGGGATCATATTGCGCCAGATTGCCATTAAAGACGTTGGCGAAAACGCGCTGATACGCTGGAAGCAGCCAGCGTTCCAGTGCCCCGGAGTCGGTATGCGCTCCAAGACCTTTGGAGGTCGTTCCGGGCGGACGGCGGCGGATACGGTCAGGGTAGATGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAG  >  NC_000913/858271‑858569

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