BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,910,295	A→T	17.1%	D7V (GAC→GTC)	yebQ →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,910,295	0	A	T	17.1%	309.2 / 60.3	135	D7V (GAC→GTC)	yebQ	putative transporter
Reads supporting (aligned to +/- strand): ref base A (49/62); new base T (9/14); total (59/76)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

Only 100 of 137 total aligned reads displayed.

ACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCC  >  NC_000913/1910196‑1910394

aCAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGa                                                                                                     <  1:210710/100‑1 (MQ=255)

    ttCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc                                                                                                 <  2:302809/100‑1 (MQ=255)

    ttCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTAACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGc                                                                                                 >  2:1501125/1‑100 (MQ=255)

     tCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGcc                                                                                                <  2:21002/100‑1 (MQ=255)

         cGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGcc                                                                                            <  1:1406058/100‑1 (MQ=255)

         cGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGcc                                                                                            <  1:2963568/100‑1 (MQ=255)

          gACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaac                                                                                           <  2:2488878/100‑4 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATGATACCCCTGTTGATTATTCACGAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaact                                                                                          <  2:890382/100‑5 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaact                                                                                          <  1:3570851/100‑5 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaact                                                                                          <  1:78703/100‑5 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaact                                                                                          <  2:2342333/100‑5 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaact                                                                                          <  2:1643213/100‑5 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaact                                                                                          <  1:836096/100‑5 (MQ=255)

           aCATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCAt                                                                                          <  1:1603406/100‑1 (MQ=255)

             aTCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTg                                                                                        <  1:2129760/100‑1 (MQ=255)

              tCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaactttt                                                                                       <  1:2907812/100‑8 (MQ=255)

              tCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGc                                                                                       <  1:94419/100‑1 (MQ=255)

              tCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGc                                                                                       >  1:3611209/1‑100 (MQ=255)

               ccACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGcc                                                                                      <  1:1787379/100‑1 (MQ=255)

               ccACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGcc                                                                                      <  1:511714/100‑1 (MQ=255)

                cACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGccc                                                                                     <  1:1644490/100‑1 (MQ=255)

                cACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGccc                                                                                     >  1:1549603/1‑100 (MQ=255)

                 aCCGTGAATATCCGTGCTAGAATCATACCCCTCTCGATTATTCACCAAAGCTATAAAATCCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaacttttggc                                                                                    <  1:2901795/100‑11 (MQ=255)

                 aCCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGTCGGCCTgaacttttggc                                                                                    <  1:2413237/100‑11 (MQ=255)

                 aCCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGcccc                                                                                    >  2:2770646/1‑100 (MQ=255)

                    gTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGc                                                                                 >  1:1080768/1‑100 (MQ=255)

                     tGAATATCAGTTCTATAAGCATACCCCTTTTGAGGATTCACCAAAGATTTAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCg                                                                                <  1:1209462/100‑1 (MQ=255)

                       aaTATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGAt                                                                              <  2:3381885/100‑1 (MQ=255)

                         tatCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATAc                                                                            <  1:3172859/100‑1 (MQ=255)

                         tatCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCAAGCGATAc                                                                            >  2:1549428/1‑100 (MQ=255)

                          atCAGTGCTAGAATCATATCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACg                                                                           >  1:2596761/1‑100 (MQ=255)

                          atCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACg                                                                           >  1:2998429/1‑100 (MQ=255)

                           tCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACgg                                                                          <  1:3280334/100‑1 (MQ=255)

                             aGTGCTAGAATCATACCCCTGTTGATTATTCACGAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGCCGGACTGCCATTGCCCCAGCGATACGGTg                                                                        >  2:2181180/1‑100 (MQ=255)

                              gTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGc                                                                       >  1:2772552/1‑100 (MQ=255)

                              gTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGc                                                                       >  1:2250343/1‑100 (MQ=255)

                                gCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGa                                                                     <  1:328670/100‑1 (MQ=255)

                                 cTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAAATCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGat                                                                    <  2:2553405/100‑1 (MQ=255)

                                  tAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGata                                                                   <  2:3302517/100‑1 (MQ=255)

                                             cccTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATtgt                                                        <  2:615899/100‑1 (MQ=255)

                                                  ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                   <  2:406425/100‑1 (MQ=255)

                                                  ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                   <  2:2378833/100‑1 (MQ=255)

                                                  ttGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTg                                                   <  1:2314895/100‑1 (MQ=255)

                                                   tGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTgg                                                  >  1:1938085/1‑100 (MQ=255)

                                                   tGATTATTCACCAAAGAAATAAAATTGATATGCAAAAAGTTAAGGCCGACGGCCTGACAATGCCCCAGAGAGAAGAGGCGATATTAAACATGCTGATTgg                                                  >  2:323336/1‑100 (MQ=255)

                                                   ngATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATtgtgcatgg                                                  >  2:68192/2‑95 (MQ=255)

                                                     attattCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTa                                                <  2:504068/100‑1 (MQ=255)

                                                      ttattCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTAt                                               >  1:3099550/1‑100 (MQ=255)

                                                        attCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTACGATATTAACCATTGTGCTTGGTa                                                >  2:2649390/1‑97 (MQ=255)

                                                        attCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTCCCCCAGCGATACGGGGCGATATTAACCATTGTGATTGGTAttt                                             >  2:1645081/1‑100 (MQ=255)

                                                          tCACCAAAGATATAAACTTCCTATGCCAAAAGTTCAGGCCGACGACCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCg                                           <  2:2131074/100‑1 (MQ=255)

                                                           cACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGa                                          >  1:1683245/1‑100 (MQ=255)

                                                           cACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGa                                          <  2:311538/100‑1 (MQ=255)

                                                            aCCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGAt                                         <  2:2772534/100‑1 (MQ=255)

                                                             ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTcccaa                                        >  2:2073262/1‑96 (MQ=255)

                                                             ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATg                                        <  1:1609349/100‑1 (MQ=255)

                                                             ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATg                                        <  2:2178262/100‑1 (MQ=255)

                                                             ccAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATg                                        >  1:3333137/1‑100 (MQ=255)

                                                              cAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATgg                                       <  1:1428108/100‑1 (MQ=255)

                                                                 aGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCg                                    <  1:3342180/100‑1 (MQ=255)

                                                                   atatAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTc                                  >  2:1766895/1‑100 (MQ=255)

                                                                    tataAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTcc                                 <  1:2445898/100‑1 (MQ=255)

                                                                      taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCtt                               >  1:1852758/1‑100 (MQ=255)

                                                                      taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCtt                               <  2:1852749/100‑1 (MQ=255)

                                                                      taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCtt                               >  1:2783487/1‑100 (MQ=255)

                                                                      taAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCtt                               >  1:2504319/1‑100 (MQ=255)

                                                                        aaaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCtt                               >  2:76478/1‑98 (MQ=255)

                                                                         aaTTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGAc                            >  2:1443513/1‑100 (MQ=255)

                                                                            tCCTATGCCAAAAGTTCAGGCCGACGGCCTTCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACgg                          <  2:1588331/99‑1 (MQ=255)

                                                                            tCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgc                         <  1:1468302/100‑1 (MQ=255)

                                                                             ccTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcg                        <  2:1949960/100‑1 (MQ=255)

                                                                             ccTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcg                        >  1:1829048/1‑100 (MQ=255)

                                                                             ccTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcg                        <  2:1059574/100‑1 (MQ=255)

                                                                              cTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcgc                       <  2:2998406/100‑1 (MQ=255)

                                                                              cTATGCCAAAAGGTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGgcgc                       <  1:1966707/100‑1 (MQ=255)

                                                                                aTGCCAAAAGTTCAGGCCGACGGCCTGCCAGTGCCCTAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCAAGGACGgccatc                     <  2:894051/100‑5 (MQ=255)

                                                                                  gCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATc                   >  2:2789616/1‑100 (MQ=255)

                                                                                  gCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATc                   >  2:336799/1‑100 (MQ=255)

                                                                                   ccAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCg                  >  1:621151/1‑100 (MQ=255)

                                                                                    ggcaatggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGc                 >  2:1139965/9‑100 (MQ=255)

                                                                                    cAAAAGTTCAGGCCGACGGCGTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGc                 >  2:2222564/1‑100 (MQ=255)

                                                                                    cAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGc                 <  1:2005574/100‑1 (MQ=255)

                                                                                    cAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGc                 <  1:462940/100‑1 (MQ=255)

                                                                                       aggTTCGGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGcc                <  2:637310/96‑1 (MQ=255)

                                                                                       aatggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCaa              >  1:673571/6‑100 (MQ=255)

                                                                                       aatggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCaa              >  1:2664940/6‑100 (MQ=255)

                                                                                        atggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAAc             >  1:2321258/5‑100 (MQ=255)

                                                                                        atggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAAc             >  1:3188446/5‑100 (MQ=255)

                                                                                        aGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAAc             <  1:2668654/100‑1 (MQ=255)

                                                                                        aGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTGCGATGGCCGTCCTTGACGGCGCAATCGCCAAc             <  1:3451556/100‑1 (MQ=255)

                                                                                          ggcAGGCCGTCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGt           >  1:1318454/3‑100 (MQ=255)

                                                                                            cAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCg         >  2:355550/1‑100 (MQ=255)

                                                                                            cAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCg         >  1:3495566/1‑100 (MQ=255)

                                                                                            cAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCg         >  2:2660061/1‑100 (MQ=255)

                                                                                             aGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGc        <  2:3611182/100‑1 (MQ=255)

                                                                                              ggCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTc          >  2:1312053/1‑97 (MQ=255)

                                                                                               gCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGccc      <  1:855812/100‑1 (MQ=255)

                                                                                               gCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGccc      <  2:2596743/100‑1 (MQ=255)

                                                                                                 cGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTg    >  2:1570481/1‑100 (MQ=255)

                                                                                                   aCGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGcc  <  2:604538/100‑1 (MQ=255)

ACAATTCTGCGACATCCACCGTGAATATCAGTGCTAGAATCATACCCCTGTTGATTATTCACCAAAGATATAAAATTCCTATGCCAAAAGTTCAGGCCGACGGCCTGCCATTGCCCCAGCGATACGGTGCGATATTAACCATTGTGATTGGTATTTCGATGGCCGTCCTTGACGGCGCAATCGCCAACGTCGCCCTGCC  >  NC_000913/1910196‑1910394

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