BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,081,057	A→G	6.6%	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	6.6%	345.8 / 12.1	123	V69A (GTG→GCG) ‡	sbmC	DNA gyrase inhibitor
Reads supporting (aligned to +/- strand): ref base A (52/61); new base G (4/4); total (58/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-02

Only 100 of 125 total aligned reads displayed.

CAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCC  >  NC_000913/2080909‑2081205

cAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCac                                                                                                                                                     <  2:1859499/150‑1 (MQ=255)

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

       cTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgtg                                                                                                                                              <  1:1964089/150‑1 (MQ=255)

       cTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgtg                                                                                                                                              >  2:337259/1‑150 (MQ=255)

       cTATTAAAGAACTGATACCAGGGTTAAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGgtg                                                                                                                                              >  2:549199/1‑150 (MQ=255)

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

         aTTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATATCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTc                                                                                                                                            <  1:333616/150‑1 (MQ=255)

          ttAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCGTGACGGTGCCg                                                                                                                                           <  1:190831/150‑1 (MQ=255)

          ttAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCGGTAATTTCTGTCAGAATGACGCCCTCACTGTTTGCGGGATGCTTAAAGTAACCCGGCACCGACGTGACGGTGCCg                                                                                                                                           <  1:1537856/150‑1 (MQ=255)

          ttAAAGAACTGATACCAGGGTTTAGCAAAAGCATCACCGACTACACCAGCTACCGCCACCGCATACTGACCACCTGTAATTGCTGGCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCg                                                                                                                                           <  1:290428/150‑1 (MQ=255)

            aaaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCa                                                                                                                                         >  2:1213925/1‑150 (MQ=255)

             aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAg                                                                                                                                        <  1:864479/150‑1 (MQ=255)

             aaGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAg                                                                                                                                        <  2:10310/150‑1 (MQ=255)

              aGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAATGTAACCCGGCACCGTCACGACGGTGTCGCAGc                                                                                                                                       >  1:373796/1‑150 (MQ=255)

               gAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCCCGCCCGTGGCggcgcg                                                                                                                                      >  1:580000/1‑145 (MQ=255)

                    gATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCATAATGACGCCCTCACTGTGTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAAtt                                                                                                                                 >  2:1405703/1‑150 (MQ=255)

                        ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCGGTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTTAGTTTTc                                                                                                                             >  2:1583861/1‑150 (MQ=255)

                        ccAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACAACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTc                                                                                                                             >  1:623017/1‑150 (MQ=255)

                        ccAGGGTTAAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTc                                                                                                                             <  1:536416/150‑1 (MQ=255)

                         cAGGGTTTAGCAAAATCATCACCGACTACCCGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                            <  2:178637/150‑1 (MQ=255)

                         cAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                            <  2:220461/150‑1 (MQ=255)

                         cAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                            <  2:774855/150‑1 (MQ=255)

                         cAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCACATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcg                                                                                                                            <  1:1843268/150‑1 (MQ=255)

                          aGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTcgg                                                                                                                           <  1:1408684/150‑1 (MQ=255)

                             gTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCGTCGCGGGTTTCGGGAGGCGTAAGGTAACCCGGCCCCGACCGGCGGGGGGCGCCGCGGATTTTTTcggcgg                                                                                                                        >  2:1523256/1‑150 (MQ=255)

                              tttAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCggg                                                                                                                       <  1:1221331/150‑1 (MQ=255)

                                  gCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGCCCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTATAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTATTTTTTCGGCGGGGGtt                                                                                                                   >  2:875751/1‑150 (MQ=255)

                                  gCAAAATCATCACAGACTACACGAGCAACCGCCACCGCATACTGACCACCTGTAAATTCTGTCAGAATGACGCCCTCACTGTTATCGGGAAGCGTAAAGTAACCCGGAACCGTCACGACGGTGTCGCAGCGTATTTTTTCGGCGGATGtt                                                                                                                   >  1:1292268/1‑150 (MQ=255)

                                  gCAAAATCAGCAGCGACTACACGAGCTACCGCCACAGCATAATGACCACCGGTAAATTCTGTCAGAATGCCGCCCTCACTGTTTTCGGGAAGCGGAAAGTACCCCGGCACCGGCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGtt                                                                                                                   >  2:133304/1‑150 (MQ=255)

                                    aaaaTCATCACCGACTACACGAGCTACCGCCACCGCATACTGCCCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTc                                                                                                                 >  1:593552/1‑150 (MQ=255)

                                     aaaTCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCa                                                                                                                <  2:1876758/150‑1 (MQ=255)

                                      aaTCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCAt                                                                                                               <  1:1982575/150‑1 (MQ=255)

                                        tcatcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCt                                                                                                             >  1:127337/1‑150 (MQ=255)

                                        tcatcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCt                                                                                                             <  1:255504/150‑1 (MQ=255)

                                         catcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGTTGTTTCATCTg                                                                                                            <  1:1169389/150‑1 (MQ=255)

                                          atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                           >  1:1637822/1‑150 (MQ=255)

                                          atcaCCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTgg                                                                                                           <  2:525086/150‑1 (MQ=255)

                                             gccGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTACCCCGGCACCGTCACGACGGTGTCGCAGCGTATTTTTTCGGCGGGTGTTTCCTCTGGAtt                                                                                                        >  1:1015367/2‑150 (MQ=255)

                                             aCCGACTACCCGAGCTACCGCCACCGCAAACTGCCCCCCTGTAATTTCTGTCAGAATGACGCCAGCACTGTTTTCGGGAATAGAAAAGTAACACGGCACCGTCACGACGGTATCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGAtt                                                                                                        >  2:293126/1‑150 (MQ=255)

                                                gACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTc                                                                                                     <  1:428147/150‑1 (MQ=255)

                                                    ccacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTaa                                                                                                 <  2:1351990/149‑1 (MQ=255)

                                                     cacGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCATCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAAt                                                                                                <  1:187067/150‑1 (MQ=255)

                                                       cgcgCTACCGCCACCGCATACTGACGACCAGTAATTTGTGTCAGAATGACGCCCTCACTGTGTTCGGGAAGCGTAAAGTAACCCGGCGCCGTCAGGAGGGTGTCGCCGCGGAATTTTTCGGCGGGTGTTTCATCGGGATTGTCGTAAttg                                                                                              >  2:249054/4‑148 (MQ=255)

                                                                     cGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGCCAGCAACCCACTc                                                                                >  2:1696461/1‑150 (MQ=255)

                                                                     cGCATACTGAACACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCTTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGATTCATCTGGATTGTCGTAATAGACAGCAACCCACTc                                                                                <  1:764308/150‑1 (MQ=255)

                                                                      gCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTcc                                                                               <  1:1027927/150‑1 (MQ=255)

                                                                      gCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTcc                                                                               <  1:1651933/150‑1 (MQ=255)

                                                                      gCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTcc                                                                               >  1:1445371/1‑150 (MQ=255)

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

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

                                                                           cTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCTTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGGCGTAATAGACAGCAACCCACTCCTTCgg                                                                          <  1:1049047/150‑1 (MQ=255)

                                                                           cTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAAGCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGGTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCgg                                                                          <  1:354653/150‑1 (MQ=255)

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

                                                                           cTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCgg                                                                          >  1:873626/1‑150 (MQ=255)

                                                                            ttaccaccTGTAATTTCTGTCAGCAAGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGGCACGAAGGTGTCGGAGGGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGc                                                                         <  2:109385/148‑1 (MQ=255)

                                                                             taccaccTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGGCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGca                                                                        <  1:1588889/149‑1 (MQ=255)

                                                                             gACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGca                                                                        >  2:1694786/1‑150 (MQ=255)

                                                                             gACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGca                                                                        >  2:341878/1‑150 (MQ=255)

                                                                             gACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGca                                                                        >  1:1416227/1‑150 (MQ=255)

                                                                                 accTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAat                                                                    <  2:1744119/150‑1 (MQ=255)

                                                                                   cTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAATGTAACCCGGCACCGTCACGACGGTGTCGCAGCGGAATTTTTCGGCGGGTGTTTCGTCTGGATTGTCGTTATAGACAGCAACCCACTCCTTCGGCAGAatat                                                                  >  1:585474/1‑150 (MQ=255)

                                                                                    tGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCTGCACCGTCACGACGGTGTCGCAGCGTAATGTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTACTTCGGAACAATAtt                                                                 >  2:337217/1‑150 (MQ=255)

                                                                                     gTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGCCGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttt                                                                >  2:1136094/1‑150 (MQ=255)

                                                                                     gTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCACCCCACTCCGTCGGCACAATAttt                                                                >  2:1106495/1‑150 (MQ=255)

                                                                                     gTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttt                                                                >  2:598286/1‑150 (MQ=255)

                                                                                      tAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTTAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAtttt                                                               <  1:1656343/150‑1 (MQ=255)

                                                                                       aaTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttttt                                                              <  1:1940911/150‑1 (MQ=255)

                                                                                       aaTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATAttttt                                                              <  2:478435/150‑1 (MQ=255)

                                                                                        aTTTCTGTCAGAATGACGCCCTCACTGTTTTCGGTAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTg                                                             <  1:511892/150‑1 (MQ=255)

                                                                                            cTGTCAAAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTAt                                                         >  1:1580816/1‑150 (MQ=255)

                                                                                              gTCAGAATGCCGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCt                                                       >  2:1770494/1‑150 (MQ=255)

                                                                                                cAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCACCCCACTCCTTCGGCACAATATTTTTGCTATCTAc                                                     >  2:1715422/1‑150 (MQ=255)

                                                                                                cAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAc                                                     >  1:56100/1‑150 (MQ=255)

                                                                                                 aGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTACCCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTGTCATCTGGATTGTCGTAATAGCCAGCAACCCCCTCCTTCGGCACAATATTTGTGCTATCTAcc                                                    >  2:1529804/1‑150 (MQ=255)

                                                                                                 aGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAcc                                                    >  2:44223/1‑150 (MQ=255)

                                                                                                 aGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAcc                                                    >  2:37754/1‑150 (MQ=255)

                                                                                                  gAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAccc                                                   >  2:1944017/1‑150 (MQ=255)

                                                                                                  gAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTAccc                                                   <  1:1798714/150‑1 (MQ=255)

                                                                                                   aaTGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCca                                                  <  1:1361840/150‑1 (MQ=255)

                                                                                                        cGCCCTCACTGTTTTCGGTAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATGGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCAcatca                                             <  2:230317/150‑1 (MQ=255)

                                                                                                             tCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAAc                                        <  1:803390/150‑1 (MQ=255)

                                                                                                             tCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAAc                                        <  1:1930030/150‑1 (MQ=255)

                                                                                                              cACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGCCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATTATCAACt                                       <  1:1445532/150‑1 (MQ=255)

                                                                                                                 tGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCt                                    >  2:1211720/1‑150 (MQ=255)

                                                                                                                  gTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTc                                   <  2:922404/150‑1 (MQ=255)

                                                                                                                       cGGGATGCGTAAAGTAACACGGCTCAGTCACGACGGTGTCGCAGCGTAATTTTTCGGGGGGTGTATCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGc                              >  2:1280308/1‑150 (MQ=255)

                                                                                                                        gggAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGcc                             <  2:1690766/150‑1 (MQ=255)

                                                                                                                            aGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctttc                         <  2:682892/150‑1 (MQ=255)

                                                                                                                             gCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctttct                        <  1:1365728/150‑1 (MQ=255)

                                                                                                                               gTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCctttcttt                      <  1:539128/150‑1 (MQ=255)

                                                                                                                                 aaaGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTAc                    >  2:1834211/1‑150 (MQ=255)

                                                                                                                                   aGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCg                  <  1:13835/150‑1 (MQ=255)

                                                                                                                                    gTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGt                 >  2:225547/1‑150 (MQ=255)

                                                                                                                                     tAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTc                <  1:1385989/150‑1 (MQ=255)

                                                                                                                                      aaCCCGGCACCGGCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCt               <  2:1599566/150‑1 (MQ=255)

                                                                                                                                       ctgcgacACCGTCGTGACGGTGTCGCCGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTg              >  2:905738/7‑150 (MQ=255)

                                                                                                                                        tgcgacACCGTCGTGACGGTGCCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGt             >  1:1647840/6‑150 (MQ=255)

                                                                                                                                        cgcgacACCGTCGTGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGGTTCATCTTGAATGTCGTAATAGCCAGCAACCCACACCTTCCGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCATTCTTTACCCTCTGt             >  1:426375/6‑150 (MQ=255)

                                                                                                                                        cccGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGt             >  2:446100/1‑150 (MQ=255)

                                                                                                                                        cccGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCGTCGGCACAATATTTTTGCTATCTACCAACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGt             >  2:1233188/1‑150 (MQ=255)

CAAGAGACTATTAAAGAACTGATACCAGGGTTTAGCAAAATCATCACCGACTACACGAGCTACCGCCACCGCATACTGACCACCTGTAATTTCTGTCAGAATGACGCCCTCACTGTTTTCGGGAAGCGTAAAGTAACCCGGCACCGTCACGACGGTGTCGCAGCGTAATTTTTCGGCGGGTGTTTCATCTGGATTGTCGTAATAGACAGCAACCCACTCCTTCGGCACAATATTTTTGCTATCTACCCACATCATCAACTGCTCAAAGCCTTTCTTTACCGTCTGTTCCCACGGGCC  >  NC_000913/2080909‑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: `‑`