BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,339,161	T→C	D87G (GAC→GGC)	gyrA ←	DNA gyrase subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,339,161	0	T	C	99.2%	420.8 / ‑4.8	119	D87G (GAC→GGC)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (51/67); total (51/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

Only 100 of 122 total aligned reads displayed.

CATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACC  >  NC_000913/2339024‑2339292

cATGGGCAATTTTCGTCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCa                                                                                                                                   >  2:209910/1‑140 (MQ=255)

cATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCa                                                                                                                                   >  2:976620/1‑140 (MQ=255)

 aTGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAGTGGATGGGCCATGAGGACGATCGTGCCAt                                                                                                                                  >  1:574577/1‑140 (MQ=255)

 aTGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACC‑AAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAAAGGCGGCGCCATGCGGACGATCGTGCCATa                                                                                                                                 >  1:764492/1‑140 (MQ=255)

      cAATTTTCGCCAGACGGATTTCCGTATAACGCGTTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATTTAACGCAGAGAGAATGCCTGCGCCATGCGGACGATCGTGCCATAGAcc                                                                                                                             >  2:750030/1‑140 (MQ=255)

      cAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAcc                                                                                                                             >  2:602400/1‑140 (MQ=255)

      cAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAcc                                                                                                                             <  2:145545/140‑1 (MQ=255)

       aaTTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccg                                                                                                                            >  2:244515/1‑140 (MQ=255)

       aaTTTTCGCCAGACGGATTTCCATATAACCCATTGCCACCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccg                                                                                                                            <  1:928611/140‑1 (MQ=255)

        tttttCGCCAGACGGATTTCCGTATACCGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccgc                                                                                                                           <  1:1413490/139‑1 (MQ=255)

        aTTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGCATGGCTGCGCCATGCGGACGAGCGTGCCATAgctgtc                                                                                                                           >  1:596039/1‑135 (MQ=255)

        aTTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccgc                                                                                                                           >  1:1623725/1‑140 (MQ=255)

         ttttCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccgcc                                                                                                                          <  1:1067377/140‑1 (MQ=255)

          tttCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAccgccg                                                                                                                         <  2:402597/140‑1 (MQ=255)

           ttCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGa                                                                                                                        >  1:1913339/1‑140 (MQ=255)

            tCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAg                                                                                                                       <  1:173502/140‑1 (MQ=255)

             cGCCAGACGGATTTCCGTATAATGCATTGCCGCCGCATAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGt                                                                                                                      <  1:836360/140‑1 (MQ=255)

               ccAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCa                                                                                                                    >  1:1533343/1‑140 (MQ=255)

                cAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCAc                                                                                                                   <  2:1913339/140‑1 (MQ=255)

                 aGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCAcc                                                                                                                  <  1:268278/140‑1 (MQ=255)

                  gACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCa                                                                                                                 <  1:1859629/140‑1 (MQ=255)

                   aCGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAg                                                                                                                       >  2:1817052/1‑133 (MQ=255)

                   aCGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCAt                                                                                                                <  1:244515/140‑1 (MQ=255)

                   aCGGATTCCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATACCGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCAt                                                                                                                <  2:1869521/140‑1 (MQ=255)

                       aTTTCCGTTTAACGCAAGGCCGCCGCAGAGAGGCCATCGCTAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATgggg                                                                                                            <  2:886106/140‑1 (MQ=255)

                       aTTTCCGTATAACGCATTGCCGGCGCAGAGTAGCCGTCGTTAGAACCGGAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATgggg                                                                                                            <  2:681334/140‑1 (MQ=255)

                       aTTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATgggg                                                                                                            <  2:28314/140‑1 (MQ=255)

                       aTTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATgggg                                                                                                            <  1:1323392/140‑1 (MQ=255)

                       aTTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATgggg                                                                                                            <  1:976620/140‑1 (MQ=255)

                        tttCCGTATAACGCGTTGCCGCCGCTGAGCCGCCGTCAATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           <  1:750030/140‑1 (MQ=255)

                        tttCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           <  1:1705647/140‑1 (MQ=255)

                        tttCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           <  1:1595660/140‑1 (MQ=255)

                        tttCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           <  1:390194/140‑1 (MQ=255)

                        tttCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           <  2:1031520/140‑1 (MQ=255)

                         acaggTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGAt                                                                                                          <  1:1189970/136‑1 (MQ=255)

                              tataACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTAt                                                                                                     <  1:1793542/140‑1 (MQ=255)

                              tataACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGAGCGTGCCACAGAGCGACGAGTCACCATGGGGATGGTAt                                                                                                     >  1:569499/1‑140 (MQ=255)

                                   cGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           >  2:746842/1‑129 (MQ=255)

                                   cGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGa                                                                                                           <  1:746842/129‑1 (MQ=255)

                                   cGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTAcc                                                                                                <  2:323102/140‑1 (MQ=255)

                                   cGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTAcc                                                                                                >  2:1829995/1‑140 (MQ=255)

                                     cATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGGCGAATCACCATGGCGATGGTATTTACCGa                                                                                              >  1:570011/1‑140 (MQ=255)

                                      tttGCCGCCGCAGAGTCGCCGTCGATAGAGCGGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGTGCCATCCGGACGTTCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTTCCGAt                                                                                             <  2:1649016/139‑1 (MQ=255)

                                      aTTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGAt                                                                                             <  1:602400/140‑1 (MQ=255)

                                      aTTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGAt                                                                                             >  1:971288/1‑140 (MQ=255)

                                        tGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGGTCGTGCCAGAGATCGCGGAGTCACCATGGGGATGGTATTTACCGATTa                                                                                           >  2:540709/1‑140 (MQ=255)

                                        tGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTAACGATTa                                                                                           >  2:1781414/1‑140 (MQ=255)

                                        tGCCGCCGCAGAGACGCCGTCGATAGAACCGAAGTTACCCTGACCGCCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCCGACGGACGATCGTGTCATAGACCCCCGAGTCACCATGGGGAAGGTAGTTACCgctgt                                                                                           >  2:1189970/1‑136 (MQ=255)

                                          ccgccgCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACg                                                                                         <  2:971288/140‑1 (MQ=255)

                                          ccgccgCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGAACGCCGAGGCACGATCGGGATGCTATTTACCGATTACg                                                                                         >  1:714710/1‑140 (MQ=255)

                                           cgccgcAGAGTTGCCCTCGATAGAACCGAAGTTCCCCTAAACGTCTAGCAGCATATAACCCAGCGAGAATGGGTGCGCCATGCGGACGATCGTGCCAAAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGt                                                                                        <  2:940897/140‑1 (MQ=255)

                                           cgccgcAGAGTCGCCTTCGATAGAATCGACGTTACCCAGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGt                                                                                        <  1:590623/140‑1 (MQ=255)

                                           cgccgcAGAGTCGCCGTCGATATAATCGACGTTACCCCGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCAGCATGGGGATGGTATTTACCGATTACGt                                                                                        <  1:576415/140‑1 (MQ=255)

                                             ccgcAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCa                                                                                      >  1:1125318/1‑140 (MQ=255)

                                                     tCGCCGTCGATAGAACCTAAGTCACCCTGACGGTCTACCGGCATATAATGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGac                                                                              <  2:716948/140‑1 (MQ=255)

                                                     tCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGac                                                                              <  2:1683554/140‑1 (MQ=255)

                                                      cGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGaca                                                                             <  1:460314/140‑1 (MQ=255)

                                                        ccGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACg                                                                           >  1:1096947/1‑140 (MQ=255)

                                                        ccGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACg                                                                           >  1:139346/1‑140 (MQ=255)

                                                          gTCGATAGAACCGAAGTCACCCCGACCGTCTCCCACCATCTAACGCAGTGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACggg                                                                         <  2:569499/140‑1 (MQ=255)

                                                           tCGATAGAACCGAAGTTTCCCTAACCGTCTACCAGCATAGAACGCAGCAAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGc                                                                        <  2:792135/140‑1 (MQ=255)

                                                           tCGATAGAACCGAAGCTACCCTCACTGTCCACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGc                                                                        <  1:540709/140‑1 (MQ=255)

                                                            cGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTAACGATTACGTCACCAACGACACGGGCa                                                                       <  1:1781414/140‑1 (MQ=255)

                                                               tAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGAt                                                                    <  1:1829995/140‑1 (MQ=255)

                                                                           tACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCt                                                        <  1:1456360/140‑1 (MQ=255)

                                                                              ccTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCACCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGAtt                                                                   >  1:1444496/1‑126 (MQ=255)

                                                                              ccTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCACCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGAtt                                                                   <  2:1444496/126‑1 (MQ=255)

                                                                                  aCCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTAGCGATTACGGCACAAACCACACGGGCAGATTTTTTATAGGCTTTGTTcc                                                 >  2:717964/1‑140 (MQ=255)

                                                                                      tCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTc                                             <  1:1369725/140‑1 (MQ=255)

                                                                                      tCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGAATACGTCACAAACAACAGGGGCAGATCTTTTATAGGCTTTGTTCCAGTc                                             >  2:682698/1‑140 (MQ=255)

                                                                                       cTACCAGCATATAACTCAGCGAGAATGGCTGCGCCATTCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCa                                            <  1:623047/140‑1 (MQ=255)

                                                                                       cTACCAGCATATAACGCAGCGAAAATGGCTGCGCCATACGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCa                                            <  1:678903/140‑1 (MQ=255)

                                                                                        tACCAGCATATAACGCATCGAGTATGGCTGCTCCATGCGTACGATCGTTCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCAt                                           <  2:574577/140‑1 (MQ=255)

                                                                                              cATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTa                                     >  2:974709/1‑140 (MQ=255)

                                                                                                       cAGCGAGAATGGCTGCGCCATGTGGACGATCATGCGATATACCGCCGAATCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCa                            <  2:570011/140‑1 (MQ=255)

                                                                                                        aGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCAt                           >  1:1110603/1‑140 (MQ=255)

                                                                                                        aGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCAt                           >  2:47184/1‑140 (MQ=255)

                                                                                                        aGCGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCGTTGGTCGAGTCATGGCCTAGTACGTTCAt                           >  1:614638/1‑140 (MQ=255)

                                                                                                          cGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGTCAGATTTTTTATAGGCGCTGTTCCAATCATACCCTAGTAGGTTCATgg                         >  2:903983/1‑140 (MQ=255)

                                                                                                          cGAGAATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATgg                         >  1:1793778/1‑140 (MQ=255)

                                                                                                           gagaATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATgtc                        >  2:1762478/1‑138 (MQ=255)

                                                                                                           gagaATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGc                        <  1:1852894/140‑1 (MQ=255)

                                                                                                            agaATGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCAt                           <  2:1445153/136‑1 (MQ=255)

                                                                                                                tGGCTGCGCCATGCGGACGATCTTGCCATCGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTaaa                   <  1:717964/140‑1 (MQ=255)

                                                                                                                tGGCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTaaa                   <  2:1664575/140‑1 (MQ=255)

                                                                                                                gacaggcgcCATGCGGATGATCTTGCCATAGTCCGCCGAGTCACCATGTGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTaaa                   <  2:593248/135‑1 (MQ=255)

                                                                                                                  gCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGt                 <  2:1533343/140‑1 (MQ=255)

                                                                                                                  gCTGCGCCATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGt                 >  1:308524/1‑140 (MQ=255)

                                                                                                                        ccATGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGa           <  1:47184/140‑1 (MQ=255)

                                                                                                                        ccATGCGGACGATCGAGCCATAGACCGCCTAGTCACCGTGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGa           <  1:911813/140‑1 (MQ=255)

                                                                                                                          aTGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACg         <  2:1530377/140‑1 (MQ=255)

                                                                                                                          aTGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACg         <  1:223687/140‑1 (MQ=255)

                                                                                                                          aTGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTCCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACg         <  2:856574/140‑1 (MQ=255)

                                                                                                                           tGCGGACGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCGTTGTTCCAGTCAGTTCATAGTACGTTCATGGAGTAAAGTGCGCGACgg        >  2:911813/1‑140 (MQ=255)

                                                                                                                            gCGGACGATCGTGCCCTAGACCGCCCAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGgt       <  1:682698/140‑1 (MQ=255)

                                                                                                                                aCGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACTTTCATGCCGTAAAGTACGCGACggctgtc   >  1:593248/1‑135 (MQ=255)

                                                                                                                                aCGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGa           <  1:361784/132‑1 (MQ=255)

                                                                                                                                aCGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGa           >  2:361784/1‑132 (MQ=255)

                                                                                                                                aCGATCGTGCCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTAc   <  1:262225/140‑1 (MQ=255)

                                                                                                                                 cGATCGTGCCATAGACCCCCGATTCACCATGAGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTAcc  <  2:585807/140‑1 (MQ=255)

CATGGGCAATTTTCGCCAGACGGATTTCCGTATAACGCATTGCCGCCGCAGAGTCGCCGTCGATAGAACCGAAGTTACCCTGACCGTCTACCAGCATATAACGCAGCGAGAATGGCTGCGCCATGCGGACGATCGTGTCATAGACCGCCGAGTCACCATGGGGATGGTATTTACCGATTACGTCACCAACGACACGGGCAGATTTTTTATAGGCTTTGTTCCAGTCATTGCCTAGTACGTTCATGGCGTAAAGTACGCGACGGTGTACC  >  NC_000913/2339024‑2339292

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