BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,724,413	A→T	7.6%	S96C (AGC→TGC)	yfiM →	protein YfiM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,724,413	0	A	T	7.6%	431.4 / 15.8	140	S96C (AGC→TGC)	yfiM	protein YfiM
Reads supporting (aligned to +/- strand): ref base A (60/69); new base T (3/8); total (63/77)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

CACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAGC  >  NC_000913/2724264‑2724562

 acaGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAg                                                                                                                                                      <  1:810112/150‑1 (MQ=255)

  caGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGc                                                                                                                                                     <  2:1678542/150‑1 (MQ=255)

  aaGCATCAGGGGATGATCCGGGATCGCAGGGCCATGTTTGGATTTATGTTCTCTGTCAGTTTGGTGGCTGCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGc                                                                                                                                                     <  2:2364372/149‑1 (MQ=255)

   aGCATCATGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCa                                                                                                                                                    <  1:2025090/150‑1 (MQ=255)

   aGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCa                                                                                                                                                    <  2:1389041/150‑1 (MQ=255)

    gCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCAc                                                                                                                                                   <  1:348092/150‑1 (MQ=255)

    gCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCAc                                                                                                                                                   >  1:461808/1‑150 (MQ=255)

       tCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCAGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCgg                                                                                                                                                <  1:1721883/150‑1 (MQ=255)

        cAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGc                                                                                                                                               <  1:530932/150‑1 (MQ=255)

        cAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGc                                                                                                                                               <  2:1621999/150‑1 (MQ=255)

        cAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGc                                                                                                                                               <  1:1005976/150‑1 (MQ=255)

        cAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGc                                                                                                                                               <  1:2352505/150‑1 (MQ=255)

        cAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGc                                                                                                                                               >  1:1174469/1‑150 (MQ=255)

          ggggATGAGCCGGGATCGCAGTGCCTTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCTTGCACCGGcga                                                                                                                                             <  1:155073/150‑3 (MQ=255)

          ggggATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCta                                                                                                                                             <  2:489526/150‑1 (MQ=255)

             gATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATAc                                                                                                                                          <  1:1444494/150‑1 (MQ=255)

               tGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCg                                                                                                                                        >  2:1141181/1‑150 (MQ=255)

                 aGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTc                                                                                                                                      <  1:442920/150‑1 (MQ=255)

                  gCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCt                                                                                                                                     >  2:818513/1‑150 (MQ=255)

                   ccGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTg                                                                                                                                    >  1:1175737/1‑150 (MQ=255)

                          cGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTg                                                                                                                             <  1:2002100/150‑1 (MQ=255)

                             aGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGAcc                                                                                                                          <  2:876036/150‑1 (MQ=255)

                             aGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGAcc                                                                                                                          <  2:461808/150‑1 (MQ=255)

                               tGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCg                                                                                                                        <  2:1564071/150‑1 (MQ=255)

                               tGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCg                                                                                                                        >  1:2259108/1‑150 (MQ=255)

                               tGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCg                                                                                                                        <  2:541718/150‑1 (MQ=255)

                                  cATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTAGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCa                                                                                                                     >  2:1986713/1‑150 (MQ=255)

                                   aTGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCAc                                                                                                                    <  1:2154767/150‑1 (MQ=255)

                                        tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAg                                                                                                               <  2:297948/150‑1 (MQ=255)

                                        tGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAg                                                                                                               <  2:2259108/150‑1 (MQ=255)

                                         ggATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGa                                                                                                              >  1:1754618/1‑150 (MQ=255)

                                          gATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGAc                                                                                                             >  2:1843414/1‑150 (MQ=255)

                                              gATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCAt                                                                                                         >  1:620953/1‑150 (MQ=255)

                                              gATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCAt                                                                                                         >  2:1966123/1‑150 (MQ=255)

                                                  ttCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATcccc                                                                                                     >  1:1860494/1‑150 (MQ=255)

                                                     gctGTCAGGTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttc                                                                                                  <  1:1362728/149‑1 (MQ=255)

                                                      ctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttcc                                                                                                 <  2:1316169/150‑1 (MQ=255)

                                                      ctGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCcccttcc                                                                                                 >  2:69548/1‑150 (MQ=255)

                                                        gTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGCCGCATCcccttccct                                                                                               >  1:776315/1‑150 (MQ=255)

                                                           cgTTTGGGGGCGGCAAAAGAGCGTTGGGATAGCCG‑CCCCGAAGGGAGCGGCGGGAGCTGGAAGGATTTGGCCTGGGATGGCGCCGGGGCAAGCACCGGCTATACCGTCTGGCAACTGCCCCGTCACTAAAGACGCATCCCCTTCCCTTTg                                                                                            <  2:2297534/149‑1 (MQ=255)

                                                           aGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTg                                                                                            <  1:296890/150‑1 (MQ=255)

                                                           aGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTg                                                                                            >  1:19698/1‑150 (MQ=255)

                                                           aGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTg                                                                                            >  2:2260459/1‑150 (MQ=255)

                                                             tttGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCg                                                                                          <  2:2235215/150‑1 (MQ=255)

                                                               tGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCCCTAAAGACGCATCCCCTTCACTTTGCGAt                                                                                        >  2:1137446/1‑150 (MQ=255)

                                                               tGGGGGCGTCAAAACAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGAt                                                                                        <  2:2163482/150‑1 (MQ=255)

                                                                  gggCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTa                                                                                     >  2:2166456/1‑150 (MQ=255)

                                                                   ggCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAg                                                                                    <  2:911107/150‑1 (MQ=255)

                                                                      gTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCAt                                                                                 <  1:1986806/150‑1 (MQ=255)

                                                                       tCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATc                                                                                <  1:833374/150‑1 (MQ=255)

                                                                        cAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCa                                                                               <  1:2013406/150‑1 (MQ=255)

                                                                        aaaaaGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGTAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGTCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTTCGATGTAGCATCa                                                                               <  1:291968/149‑1 (MQ=255)

                                                                           aaGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAg                                                                            <  1:2150745/150‑1 (MQ=255)

                                                                                cTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCCTCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAcc                                                                       >  1:2482516/1‑150 (MQ=255)

                                                                                cTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAcc                                                                       >  2:166551/1‑150 (MQ=255)

                                                                                cTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAAcc                                                                       >  1:76612/1‑150 (MQ=255)

                                                                                  ttGGGATAGCCGCCCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCa                                                                      <  1:805208/150‑1 (MQ=255)

                                                                                   tGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGa                                                                    <  2:1816169/150‑1 (MQ=255)

                                                                                       aTAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACg                                                                <  1:587699/150‑1 (MQ=255)

                                                                                       aTAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACg                                                                <  1:579004/150‑1 (MQ=255)

                                                                                       aTAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACg                                                                <  1:1517107/150‑1 (MQ=255)

                                                                                         aGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGcc                                                              <  1:1076517/150‑1 (MQ=255)

                                                                                          gccg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGcca                                                             <  1:1698998/150‑1 (MQ=255)

                                                                                          gccg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGcca                                                             <  1:1639652/150‑1 (MQ=255)

                                                                                          gccg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGcca                                                             >  2:2328901/1‑150 (MQ=255)

                                                                                           ccg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccac                                                            >  1:1258198/1‑150 (MQ=255)

                                                                                           ccg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccac                                                            >  1:2310324/1‑150 (MQ=255)

                                                                                            cg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccacc                                                           <  1:1141181/150‑1 (MQ=255)

                                                                                            cg‑ccCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGccacc                                                           >  2:515440/1‑150 (MQ=255)

                                                                                                cccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACg                                                        >  2:2285851/1‑150 (MQ=255)

                                                                                                 ccGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACgg                                                       >  1:1012631/1‑150 (MQ=255)

                                                                                                   gAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGcc                                                     <  1:1648168/150‑1 (MQ=255)

                                                                                                     aGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCAt                                                   >  1:43959/1‑150 (MQ=255)

                                                                                                       ggAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCa                                                 <  1:660846/150‑1 (MQ=255)

                                                                                                          gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAgg                                              >  2:2079959/1‑150 (MQ=255)

                                                                                                          gCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAgg                                              >  2:2088329/1‑150 (MQ=255)

                                                                                                            ggCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTc                                            >  2:2356770/1‑150 (MQ=255)

                                                                                                            ggCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTc                                            <  2:2484157/150‑1 (MQ=255)

                                                                                                            ggCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTc                                            >  1:1011772/1‑150 (MQ=255)

                                                                                                            ggCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTc                                            >  1:357482/1‑150 (MQ=255)

                                                                                                            ggCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTc                                            >  1:2450916/1‑150 (MQ=255)

                                                                                                                   gctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATAcc                                     <  2:9095/150‑1 (MQ=255)

                                                                                                                   gctggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATAcc                                     >  2:751026/1‑150 (MQ=255)

                                                                                                                      ggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAACGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCaga                                  <  1:1144102/150‑1 (MQ=255)

                                                                                                                      ggaAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCaga                                  >  1:519969/1‑150 (MQ=255)

                                                                                                                           gATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGaagaag                             >  1:793193/1‑150 (MQ=255)

                                                                                                                            aTTTGGCCTGGGATGTCGCCTGTGCAAGCACAGGCTATACAGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAgg                            >  2:504112/1‑150 (MQ=255)

                                                                                                                                ggCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGt                        >  1:1177655/1‑150 (MQ=255)

                                                                                                                                 gCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGtt                       >  2:1980565/1‑150 (MQ=255)

                                                                                                                                    tGGGAGGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTcc                    <  2:1750560/150‑1 (MQ=255)

                                                                                                                                      ggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCAt                  <  2:1888273/150‑1 (MQ=255)

                                                                                                                                      ggATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCAt                  <  2:620953/150‑1 (MQ=255)

                                                                                                                                         tGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTcc               <  1:2216358/150‑1 (MQ=255)

                                                                                                                                           tCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTa             <  2:1175737/150‑1 (MQ=255)

                                                                                                                                            cGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt            >  2:869193/1‑150 (MQ=255)

                                                                                                                                            cGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGAAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt            >  1:1848978/1‑150 (MQ=255)

                                                                                                                                            agCCGGTGCTTGCACCGGCTATACCGTCTGGCAACTGCCCCGTCACTAAAGACGCATCCCCTTCCCTTTCCGATGTAGCATCAAAGACACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTAt            >  1:1275542/2‑150 (MQ=255)

                                                                                                                                                    cAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAATGCTGTTTCCATTCCTATTGATTTCa    <  2:119629/150‑1 (MQ=255)

                                                                                                                                                      aGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTACAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAGc  <  1:1102607/150‑1 (MQ=255)

                                                                                                                                                      aGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAGc  >  1:2393782/1‑150 (MQ=255)

CACAGCATCAGGGGATGAGCCGGGATCGCAGTGCCATGTTTGGATTGATGTTCTCTGTCAGTTTGGGGGCGTCAAAAGAGCTTTGGGATAGCCG‑CCCCGAAGGGAGCGGCTGGAGCTGGAAGGATTTGGCCTGGGATGTCGCCGGTGCAAGCACCGGCTATACCGTCTGGCAACTGACCCGTCACTAAAGACGCATCCCCTTCCCTTTGCGATGTAGCATCAAAGAAACCAGAAACGCCACCACGGCCATCAAGGTCACATACCAGAAGAAGGCTGTTTCCATTCCTATTGATTTCAGC  >  NC_000913/2724264‑2724562

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