BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,081,058	C→T	7.7%	V69M (GTG→ATG) ‡	sbmC ←	DNA gyrase inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,081,058	0	C	T	7.7%	320.4 / 8.4	103	V69M (GTG→ATG) ‡	sbmC	DNA gyrase inhibitor
Reads supporting (aligned to +/- strand): ref base C (49/46); new base T (3/5); total (52/51)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

AAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCC  >  NC_000913/2080910‑2081205

aaGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCacg                                                                                                                                                    <  2:1376832/150‑1 (MQ=255)

  gagaCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCacgac                                                                                                                                                  <  1:230035/150‑1 (MQ=255)

    gaCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACgg                                                                                                                                                >  2:633623/1‑150 (MQ=255)

     aCTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTATTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgt                                                                                                                                               >  2:1553076/1‑150 (MQ=255)

       tATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgtgt                                                                                                                                             >  1:598802/1‑150 (MQ=255)

       tATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGATACCGCCACCGCATAATGACCACCTGTAATTTCTGTCAGAATGACGCACTCACTGTTTTCGGGAAGCGTAAAGTAACCCTGCACCGTCACGACGctgt                                                                                                                                             >  1:1682488/1‑150 (MQ=255)

         ttAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGAGTAAAGTACCCCGGCACCGTCACGCCGGTGTCg                                                                                                                                           >  1:913928/1‑150 (MQ=255)

          tAACGAACTGATACCACCGGATAGCAAAATCATGACCGACTACACGAGCTAGCGCCACCGAATACTGACCACCTGTAATATCTGTCAGAATGCCGCCCTCACTGTTTTCGGGAAGCGTCAAGTAACCCGACACCGGCGTGACGGTGCCgg                                                                                                                                          <  1:271902/150‑2 (MQ=255)

          tAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCgg                                                                                                                                          <  1:123312/150‑2 (MQ=255)

          tAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGc                                                                                                                                          >  1:482270/1‑150 (MQ=255)

          tAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGc                                                                                                                                          <  2:661524/150‑1 (MQ=255)

          tAAAGAACTGAAACCAGGGTTAAGCAAAATCATCACCGACTACACCAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCgg                                                                                                                                          <  1:1715433/150‑2 (MQ=255)

           aaaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGGCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGACACCGTCGTGACGGTGCCggg                                                                                                                                         <  1:123999/150‑3 (MQ=255)

            aaGAACAGATACCCGTGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCgggt                                                                                                                                        <  2:617150/150‑4 (MQ=255)

                   gATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGTGTCGCAGCGGAAtt                                                                                                                                 >  2:467796/1‑150 (MQ=255)

                     tACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTGCCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAtttt                                                                                                                               >  2:1647357/1‑150 (MQ=255)

                     tACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCATCACGACGGTGTCGCAGCGTAAtttt                                                                                                                               >  1:190728/1‑150 (MQ=255)

                       cccgggTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGA‑CACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                            <  2:44994/147‑1 (MQ=255)

                       ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTc                                                                                                                             >  2:1686813/1‑150 (MQ=255)

                       ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTc                                                                                                                             >  2:254706/1‑150 (MQ=255)

                               tAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGgtg                                                                                                                     <  2:1670271/150‑1 (MQ=255)

                                 gCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCCGCGTAATTTTTAGGCGGGTGtt                                                                                                                   >  2:1084937/1‑150 (MQ=255)

                                 gCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGtt                                                                                                                   >  1:1706289/1‑150 (MQ=255)

                                  cAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGttt                                                                                                                  <  2:1706289/150‑1 (MQ=255)

                                   aaaaTCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTc                                                                                                                 <  1:839458/150‑1 (MQ=255)

                                        catcaCCGACTACACGTGCTACCGCCACCGAATAGTGACAACCTTTAATTTGTGTCTGAATGAGGCCCTCAGTGTTTTCGGTAAGCTTAAAGTAACCCGGCACCATCACGACGGTGTCGCAGGTTAATTTTTCGGCGGGTGTTTCATCTg                                                                                                            <  2:190728/150‑1 (MQ=255)

                                            aCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGAtt                                                                                                        <  2:516348/150‑1 (MQ=255)

                                                aCTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTTTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCg                                                                                                    <  2:1227413/150‑1 (MQ=255)

                                                aCTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCg                                                                                                    >  1:1512460/1‑150 (MQ=255)

                                                 cTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGt                                                                                                   <  2:1574417/150‑1 (MQ=255)

                                                 cTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGt                                                                                                   >  2:345816/1‑150 (MQ=255)

                                                   acacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTaa                                                                                                 >  2:786208/1‑150 (MQ=255)

                                                   acacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTaa                                                                                                 >  2:895474/1‑150 (MQ=255)

                                                   aaacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCTTAAAGTAACCCGGCACCGTCACGACGGTCTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTaa                                                                                                 >  2:745399/3‑150 (MQ=255)

                                                     acGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCTCTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATa                                                                                               <  1:254336/150‑1 (MQ=255)

                                                     acGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATa                                                                                               <  2:1041977/150‑1 (MQ=255)

                                                        aGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGATGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGAc                                                                                            <  2:1618142/150‑1 (MQ=255)

                                                        aGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGAc                                                                                            <  1:701185/150‑1 (MQ=255)

                                                         gCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACa                                                                                           <  1:633623/150‑1 (MQ=255)

                                                               gCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAcc                                                                                     <  1:817059/150‑1 (MQ=255)

                                                               gCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTACATCTGGATTGTCGTAATAGACAGCAAcc                                                                                     >  2:1202635/1‑150 (MQ=255)

                                                                ccACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAccc                                                                                    <  1:1686813/150‑1 (MQ=255)

                                                                ccACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAccc                                                                                    >  2:705625/1‑150 (MQ=255)

                                                                    cGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACATCAACCCACTc                                                                                >  2:1167216/1‑150 (MQ=255)

                                                                      cATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCTGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCt                                                                              >  2:523320/1‑150 (MQ=255)

                                                                       aTACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCtt                                                                             >  1:717111/1‑150 (MQ=255)

                                                                        tACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTc                                                                            <  2:598802/150‑1 (MQ=255)

                                                                        tACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTc                                                                            >  1:1416605/1‑150 (MQ=255)

                                                                        tACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTc                                                                            >  1:1042405/1‑150 (MQ=255)

                                                                        tACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTc                                                                            <  1:254706/150‑1 (MQ=255)

                                                                        tACTGACCACCTGTAATTTCTGTCAGAATGACGACCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTATTTTTTCGGCGGGTGATTCATCTGGATTGTCGTAATAGACAGCAACCCACTACtta                                                                            >  1:1231599/1‑149 (MQ=255)

                                                                         aCTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCTGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCg                                                                           <  1:380982/150‑1 (MQ=255)

                                                                         aCTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCg                                                                           <  1:2528/150‑1 (MQ=255)

                                                                         aCTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCg                                                                           <  1:527256/150‑1 (MQ=255)

                                                                          cTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCgg                                                                          <  2:501176/150‑1 (MQ=255)

                                                                          cTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTCTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCgg                                                                          <  1:1553076/150‑1 (MQ=255)

                                                                            gACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTTGATTGTCGTAATAGAGAGCAACCCACTCCTTCGGcc                                                                        >  2:924205/1‑149 (MQ=255)

                                                                              ccaccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGcaca                                                                      >  2:909038/1‑150 (MQ=255)

                                                                              ccaccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTGTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGccca                                                                      >  1:1067986/1‑150 (MQ=255)

                                                                                 ccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAata                                                                   >  1:1127551/1‑150 (MQ=255)

                                                                                       aTTTCTGTCCTAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGTGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTg                                                             <  2:774923/150‑1 (MQ=255)

                                                                                       aTTTCTGTCAGAATGACGCCCTCACTTTTTTCTGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTg                                                             <  2:1231599/150‑1 (MQ=255)

                                                                                           cTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTAt                                                         >  1:123098/1‑150 (MQ=255)

                                                                                            tGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATc                                                        >  2:641463/1‑150 (MQ=255)

                                                                                             gTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCt                                                       >  1:954322/1‑150 (MQ=255)

                                                                                               cAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGATATCTAc                                                     >  1:1217086/1‑150 (MQ=255)

                                                                                                  aaTGACGCCCTCACTGTTTTCTGTAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCca                                                  <  2:577050/150‑1 (MQ=255)

                                                                                                    tGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCcaca                                                <  2:1512460/150‑1 (MQ=255)

                                                                                                       cGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCAcatca                                             >  2:1334147/1‑150 (MQ=255)

                                                                                                         cccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACatcatc                                           <  1:1647357/150‑1 (MQ=255)

                                                                                                          ccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACAtcatca                                          <  1:895474/150‑1 (MQ=255)

                                                                                                          ccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACAtcatca                                          <  1:1068211/150‑1 (MQ=255)

                                                                                                          ccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACAtcatca                                          <  1:1167216/150‑1 (MQ=255)

                                                                                                           cTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCaa                                         <  1:540996/150‑1 (MQ=255)

                                                                                                             cACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTTGGCACAATATTTTTGCTATCTACCCACATCATCAACt                                       >  2:387636/1‑150 (MQ=255)

                                                                                                               cTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGc                                     <  1:1450602/150‑1 (MQ=255)

                                                                                                               cTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGc                                     <  1:352046/150‑1 (MQ=255)

                                                                                                               cTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCAACACACATCATCAACTGc                                     >  2:1582937/1‑150 (MQ=255)

                                                                                                                        ggAAGCGTAAAGTAACCCGGCACCGTCACGTCGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCaa                                 <  2:616982/145‑1 (MQ=255)

                                                                                                                        ggAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCCCAATATTTTTGCTATCTACCCACATCAGCAACTGCTCAAAGcgt                            >  2:302511/1‑148 (MQ=255)

                                                                                                                        agAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCaa                                 >  1:616982/2‑145 (MQ=255)

                                                                                                                           aGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTTTCTACCCACATCATCAACTGCTCAAAGCctttc                         <  1:786208/150‑1 (MQ=255)

                                                                                                                            gCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCCTCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAACGCctttct                        >  1:1125724/1‑150 (MQ=255)

                                                                                                                              gTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctttcttt                      >  1:1084494/1‑150 (MQ=255)

                                                                                                                                aaaGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACAGCTCAAAGCCTTTCTTTAc                    <  2:1127551/150‑1 (MQ=255)

                                                                                                                                  aGTAACCCGGCACCGTCTCGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCg                  <  2:913928/150‑1 (MQ=255)

                                                                                                                                      ctgcgacACCGTCGTGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCGACCCACATCATCAACTGCTCAAAGCCGTTCTTTACCGTCTg              >  1:678813/7‑150 (MQ=255)

                                                                                                                                      ctgcgacACCGTCGTGACGGTGCCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCtt              >  1:575203/7‑149 (MQ=255)

                                                                                                                                      aCCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTg              <  2:963459/150‑1 (MQ=255)

                                                                                                                                        gcgacACCGTCGTGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGtt            >  1:987072/5‑150 (MQ=255)

                                                                                                                                        ccGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAAACCACACCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGtt            >  1:916203/1‑150 (MQ=255)

                                                                                                                                         cGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTc           <  2:1042405/150‑1 (MQ=255)

                                                                                                                                          ggCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGGAATAGACATCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACGCACATCATCAACTGCTCAAAGCCTTTCTTTTCCGTCTGttac          >  2:308013/1‑148 (MQ=255)

                                                                                                                                            cACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCa        >  1:1497424/1‑150 (MQ=255)

                                                                                                                                             aCCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCAc       >  1:790416/1‑150 (MQ=255)

                                                                                                                                             aCCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTGCCCAc       <  1:389050/150‑1 (MQ=255)

                                                                                                                                                 tgacgacgGTGTAGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGc   <  2:1682488/148‑1 (MQ=255)

                                                                                                                                                 tCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGc   <  1:705625/150‑1 (MQ=255)

                                                                                                                                                  cACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGcc  <  2:123098/150‑1 (MQ=255)

                                                                                                                                                  cACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGcc  >  2:1107409/1‑150 (MQ=255)

AAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCC  >  NC_000913/2080910‑2081205

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