BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,825,188	T→C	M369V (ATG→GTG)	gltS ←	glutamate transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,825,188	0	T	C	83.9%	721.8 / 136.1	322	M369V (ATG→GTG)	gltS	glutamate transport
Reads supporting (aligned to +/- strand): ref base T (27/25); new base C (154/116); total (181/141)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 328 total aligned reads displayed.

AAAAATCGGCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGCAGC  >  NC_000913/3825096‑3825278

    aTCGGCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGc                                                                                        <  1:5776107/93‑1 (MQ=255)

    aTCGGCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGc                                                                                        <  1:6374463/93‑1 (MQ=255)

    aTCGGCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGc                                                                                        <  1:4060276/93‑1 (MQ=255)

      cGGCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGa                                                                                      <  1:3432722/93‑1 (MQ=255)

       ggCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGAc                                                                                     >  1:2752450/1‑93 (MQ=255)

        gCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACg                                                                                    <  2:7218834/93‑1 (MQ=255)

         cAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACgg                                                                                   >  2:2960597/1‑93 (MQ=255)

            cATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGcc                                                                                >  1:7112559/1‑93 (MQ=255)

             aTCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCa                                                                               >  1:7197405/1‑93 (MQ=255)

             aTCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCa                                                                               >  2:5533283/1‑93 (MQ=255)

             aTCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCa                                                                               >  1:1616489/1‑93 (MQ=255)

             aTCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCa                                                                               >  2:5751696/1‑93 (MQ=255)

              tCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCaa                                                                              <  1:2696159/93‑1 (MQ=255)

              tCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACACCACGTGCGACGGGCCaa                                                                              >  2:2211490/1‑93 (MQ=255)

              tCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACACCACGTGCGACGGGCCaa                                                                              >  1:3662010/1‑93 (MQ=255)

               cAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCaaa                                                                             <  2:1085566/93‑1 (MQ=255)

               cAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCaaa                                                                             <  2:190331/93‑1 (MQ=255)

                aaaTACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAg                                                                            >  1:4737018/1‑93 (MQ=255)

                aaaTACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAg                                                                            >  2:275936/1‑93 (MQ=255)

                aaaTACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACACCACGTGCGACGGGCCAAAg                                                                            <  2:5222339/93‑1 (MQ=255)

                 aaTACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGc                                                                           <  1:3540496/93‑1 (MQ=255)

                   tACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGt                                                                         <  2:2636499/93‑1 (MQ=255)

                   tACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGt                                                                         >  2:3382897/1‑93 (MQ=255)

                    aCAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGtt                                                                        <  2:4372264/93‑1 (MQ=255)

                      aaCTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCa                                                                      <  2:1504817/93‑1 (MQ=255)

                      aaCTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCa                                                                      <  2:6297434/93‑1 (MQ=255)

                       aCTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAg                                                                     >  1:4218780/1‑93 (MQ=255)

                        cTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGt                                                                    >  1:1405159/1‑93 (MQ=255)

                        cTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGt                                                                    >  2:7012451/1‑93 (MQ=255)

                         ttAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTg                                                                   >  2:4978668/1‑93 (MQ=255)

                         ttAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACACCACGTGCGACGGGCCAAAGCGTTCAGTg                                                                   <  1:3016120/93‑1 (MQ=255)

                             ttACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCg                                                               >  2:1296929/1‑93 (MQ=255)

                              tACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGc                                                              >  2:775452/1‑93 (MQ=255)

                                ccAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCt                                                            >  1:2281810/1‑93 (MQ=255)

                                  aGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGc                                                          >  1:3809683/1‑93 (MQ=255)

                                   gcgcATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCa                                                         >  2:6563466/1‑93 (MQ=255)

                                     gcATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATg                                                       <  2:6793861/93‑1 (MQ=255)

                                      cATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGt                                                      <  1:6690204/93‑1 (MQ=255)

                                      cATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGt                                                      <  1:5631684/93‑1 (MQ=255)

                                      cATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGt                                                      <  2:7156756/93‑1 (MQ=255)

                                      cATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGt                                                      <  1:923162/93‑1 (MQ=255)

                                       aTTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGtt                                                     <  2:4030957/93‑1 (MQ=255)

                                        ttGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTg                                                    <  2:6168476/93‑1 (MQ=255)

                                         tGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTgg                                                   <  1:2199651/93‑1 (MQ=255)

                                         tGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTgg                                                   <  1:2830115/93‑1 (MQ=255)

                                         tGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTgg                                                   <  2:7467516/93‑1 (MQ=255)

                                          gACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGc                                                  >  1:6050513/1‑93 (MQ=255)

                                          gACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGc                                                  >  2:6357686/1‑93 (MQ=255)

                                          gACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGc                                                  <  2:4573355/93‑1 (MQ=255)

                                          gACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGc                                                  >  2:3955437/1‑93 (MQ=255)

                                          gACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGc                                                  >  2:7259281/1‑93 (MQ=255)

                                           aCGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCg                                                 <  1:248997/93‑1 (MQ=255)

                                             gATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGAt                                               <  2:3807444/93‑1 (MQ=255)

                                              atatCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGAtt                                              <  1:1330131/93‑1 (MQ=255)

                                                  cGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCg                                          >  2:6253483/1‑93 (MQ=255)

                                                   gATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGt                                         >  2:3340493/1‑93 (MQ=255)

                                                   gATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGt                                         <  1:514880/93‑1 (MQ=255)

                                                    aTAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGtt                                        >  1:1086617/1‑93 (MQ=255)

                                                    aTAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGtt                                        <  1:6837035/93‑1 (MQ=255)

                                                     tAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTg                                       >  1:1867112/1‑93 (MQ=255)

                                                     tAAAGAACGCACCGACCATCGGCACCACCAAAAACACCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTg                                       >  2:4705768/1‑93 (MQ=255)

                                                     tAAACAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTg                                       >  2:7371773/1‑93 (MQ=255)

                                                       aaGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGc                                     >  1:5595596/1‑93 (MQ=255)

                                                       aaGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGc                                     >  1:5371852/1‑93 (MQ=255)

                                                        aGAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCg                                    >  2:2125537/1‑93 (MQ=255)

                                                        aGAACGCACCGACCATCGGCACCACCAAAAACGCCACGCGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCg                                    <  1:6330944/93‑1 (MQ=255)

                                                         gAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGt                                   <  2:5999668/93‑1 (MQ=255)

                                                         gAACGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGt                                   >  2:7101415/1‑93 (MQ=255)

                                                         gAACGCACCGACCATCGGCACCACCAAAAACACCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGt                                   <  1:4325790/93‑1 (MQ=255)

                                                          aaCGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGtt                                  <  1:7329571/93‑1 (MQ=255)

                                                           aCGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTg                                 >  1:4475214/1‑93 (MQ=255)

                                                            cGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGc                                <  2:6916818/93‑1 (MQ=255)

                                                            cGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGc                                <  2:2444095/93‑1 (MQ=255)

                                                            cGCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGc                                <  2:5494011/93‑1 (MQ=255)

                                                             gCACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCa                               >  1:3257110/1‑93 (MQ=255)

                                                              cACCGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCAc                              >  1:7447743/1‑93 (MQ=255)

                                                                ccGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCg                            >  2:3227909/1‑93 (MQ=255)

                                                                ccGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCg                            >  2:1805714/1‑93 (MQ=255)

                                                                 cGACCATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGa                           <  1:6467737/93‑1 (MQ=255)

                                                                  gACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTACACCCAg                          >  2:5611999/1‑93 (MQ=255)

                                                                    ccATCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGc                        >  2:446302/1‑93 (MQ=255)

                                                                       tCGGCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCaa                     <  2:2019383/93‑1 (MQ=255)

                                                                         ggCACCACCAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCaaaa                   >  1:1286417/1‑93 (MQ=255)

                                                                            accaccAAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACca                >  2:4939358/1‑93 (MQ=255)

                                                                              caccaAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACcaca              <  1:3040880/93‑1 (MQ=255)

                                                                               accaAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAg             >  1:5006535/1‑93 (MQ=255)

                                                                               accaAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAg             <  1:4189821/93‑1 (MQ=255)

                                                                               accaAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAg             <  2:4645274/93‑1 (MQ=255)

                                                                                 caAAAACGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTg           >  1:2132603/1‑93 (MQ=255)

                                                                                  aaaaaCGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGa          <  1:6917444/93‑1 (MQ=255)

                                                                                   aaaaCGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGAc         >  2:3944448/1‑93 (MQ=255)

                                                                                   aaaaCGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGAc         >  1:6569317/1‑93 (MQ=255)

                                                                                    aaaCGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGAcc        >  2:4138434/1‑93 (MQ=255)

                                                                                     aaCGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGAccc       <  2:1749705/93‑1 (MQ=255)

                                                                                       cGCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGc     >  2:4524005/1‑93 (MQ=255)

                                                                                        gCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGca    >  1:2010372/1‑93 (MQ=255)

                                                                                        gCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGca    >  1:1743537/1‑93 (MQ=255)

                                                                                        gCCACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGca    <  2:501624/93‑1 (MQ=255)

                                                                                         ccACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGcag   >  2:3479844/1‑93 (MQ=255)

                                                                                         ccACGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGcag   >  1:6030554/1‑93 (MQ=255)

AAAAATCGGCAACATCAAATACAACTTAATTACCAGCGCATTGACGATATCGATAAAGAACGCACCGACCATCGGCACCACCAAAAACGCCATGTGCGACGGGCCAAAGCGTTCAGTGATCGCCTGCATGTTGGCGATTGCCGTTGGCGTTGCACCGAGGCCAAAACCACAGTGACCCGCAGC  >  NC_000913/3825096‑3825278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`