BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,081,057	A→G	7.9%	V69A (GTG→GCG) ‡	sbmC ←	DNA gyrase inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,081,057	0	A	G	7.9%	279.9 / 12.8	103	V69A (GTG→GCG) ‡	sbmC	DNA gyrase inhibitor
Reads supporting (aligned to +/- strand): ref base A (49/45); new base G (3/5); total (52/51)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-02

Only 100 of 107 total aligned reads displayed.

GCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCC  >  NC_000913/2080908‑2081205

gCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCa                                                                                                                                                      >  1:661524/1‑150 (MQ=255)

  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:254706/1‑150 (MQ=255)

                          cAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                            >  1:1574417/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)

                                            tcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGa                                                                                                          >  2:540996/1‑150 (MQ=255)

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

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

                                                  aCTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCg                                                                                                    <  2:1584312/150‑1 (MQ=255)

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

                                                   cTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGt                                                                                                   <  2:1574417/150‑1 (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)

                                                                  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                                                                            >  1:1042405/1‑150 (MQ=255)

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

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

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

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

                                                                           aCTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCg                                                                           <  1:2528/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:963459/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:1167216/150‑1 (MQ=255)

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

                                                                                                            ccTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACAtcatca                                          <  1:895474/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)

                                                                                                                    ttttCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCa                                  <  1:558872/150‑1 (MQ=255)

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

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

                                                                                                                          agAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCaa                                 >  1:616982/2‑145 (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:1107409/1‑150 (MQ=255)

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

GCAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCC  >  NC_000913/2080908‑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: `‑`