BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,096	A→G	L251P (CTC→CCC)	ylaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,096	0	A	G	98.6%	890.1 / ‑2.8	289	L251P (CTC→CCC)	ylaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (127/158); total (130/159)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

ATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477006‑477185

 tAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCg                                                                                        >  2:4790541/1‑93 (MQ=255)

  aGCATCCGGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCgg                                                                                       <  1:1525928/93‑1 (MQ=255)

     aTCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAACACGCAGCACAAAAATCGCTGCAAGCAGGCCAATCACCAAACGGGCCGGTaa                                                                                    >  1:5477242/1‑93 (MQ=255)

      tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  2:6923807/1‑93 (MQ=255)

       ccTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAcc                                                                                  <  1:1003117/93‑1 (MQ=255)

       ccTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAcc                                                                                  <  1:6733703/93‑1 (MQ=255)

       ccTCCAGCCGATGATGCGGTTACTGAATACGGCACAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAcc                                                                                  <  2:6817938/93‑1 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCa                                                                                 >  1:2699160/1‑93 (MQ=255)

         tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaa                                                                                <  1:6446415/93‑1 (MQ=255)

         tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaa                                                                                >  1:370096/1‑93 (MQ=255)

          gCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaaa                                                                               >  1:5038928/1‑93 (MQ=255)

           cAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAAt                                                                              <  1:7800754/93‑1 (MQ=255)

            aGCCGATGATGCGGTGACTGAATACGGCACAAAATACGCAGCACAACCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATa                                                                             <  1:558237/93‑1 (MQ=255)

             gCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaa                                                                            <  1:2255641/93‑1 (MQ=255)

             gCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaa                                                                            >  2:4465830/1‑93 (MQ=255)

              ccGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGTCGGTAACCAAATaaa                                                                           <  2:3324864/93‑1 (MQ=255)

              ccGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaaa                                                                           <  2:123689/93‑1 (MQ=255)

               cGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAg                                                                          <  2:7626370/93‑1 (MQ=255)

                 atgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGAc                                                                        >  2:3159041/1‑93 (MQ=255)

                  tgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACt                                                                       >  2:5312039/1‑93 (MQ=255)

                  tgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACt                                                                       <  2:4715092/93‑1 (MQ=255)

                  tgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACt                                                                       <  1:4764296/93‑1 (MQ=255)

                   gatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACtt                                                                      <  2:1061281/93‑1 (MQ=255)

                   gatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACtt                                                                      <  1:5406414/93‑1 (MQ=255)

                   gatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACtt                                                                      <  2:637578/93‑1 (MQ=255)

                    atgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTg                                                                     <  2:2184557/93‑1 (MQ=255)

                     tgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGc                                                                    <  2:174436/93‑1 (MQ=255)

                      gCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGcc                                                                   <  2:783197/93‑1 (MQ=255)

                         gTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGGCGAt                                                                <  2:5045784/93‑1 (MQ=255)

                         gTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGAt                                                                <  2:2207581/93‑1 (MQ=255)

                         gTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGAt                                                                >  1:7482994/1‑93 (MQ=255)

                          tGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATg                                                               <  1:1670767/93‑1 (MQ=255)

                          tGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATg                                                               >  2:354075/1‑93 (MQ=255)

                           gACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGc                                                              >  2:4888915/1‑93 (MQ=255)

                             cTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCa                                                            <  2:6196334/93‑1 (MQ=255)

                             cTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCa                                                            <  2:5879066/93‑1 (MQ=255)

                              tGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAg                                                           >  2:3618817/1‑93 (MQ=255)

                              tGAATACGGCGCAAAATAAGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAg                                                           >  1:6290157/1‑93 (MQ=255)

                                aaTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGcc                                                         >  1:5772560/1‑93 (MQ=255)

                                aaTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGcc                                                         <  1:7542161/93‑1 (MQ=255)

                                aaTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGcc                                                         <  2:77854/93‑1 (MQ=255)

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        <  1:4401438/93‑1 (MQ=255)

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        >  1:7934212/1‑93 (MQ=255)

                                   aCGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCttt                                                      <  2:6622588/93‑1 (MQ=255)

                                       cgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGc                                                  >  2:5780835/1‑93 (MQ=255)

                                          aaaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAAc                                               >  2:2506361/1‑93 (MQ=255)

                                          aaaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAAc                                               <  1:3455522/93‑1 (MQ=255)

                                          aaaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAAc                                               >  2:1833605/1‑93 (MQ=255)

                                             aTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtt                                            <  1:3379280/93‑1 (MQ=255)

                                               aCGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtttt                                          >  2:7736516/1‑93 (MQ=255)

                                                cGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                         <  1:6193084/93‑1 (MQ=255)

                                                cGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                         <  1:1402650/93‑1 (MQ=255)

                                                cGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                         <  1:3543943/93‑1 (MQ=255)

                                                  cagcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGt                                       <  1:3936683/93‑1 (MQ=255)

                                                    gcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTg                                     >  1:3634334/1‑93 (MQ=255)

                                                    gcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTg                                     >  1:7184909/1‑93 (MQ=255)

                                                     cacaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGa                                    >  1:3292356/1‑93 (MQ=255)

                                                     cacaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGa                                    <  1:1768074/93‑1 (MQ=255)

                                                     cacaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGa                                    <  1:2483828/93‑1 (MQ=255)

                                                     cacaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGa                                    >  1:3088707/1‑93 (MQ=255)

                                                       caAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAg                                  >  2:2791700/1‑93 (MQ=255)

                                                        aaaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGc                                 >  1:6985161/1‑93 (MQ=255)

                                                         aaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGGGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGTAgtc                                >  1:4606798/1‑91 (MQ=255)

                                                          aCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGccc                               >  1:6939979/1‑93 (MQ=255)

                                                          aCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGccc                               >  1:6516740/1‑93 (MQ=255)

                                                           cATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCa                              >  1:5360307/1‑93 (MQ=255)

                                                           cATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCa                              <  2:6051340/93‑1 (MQ=255)

                                                           cATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCa                              <  1:7392052/93‑1 (MQ=255)

                                                           cATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCa                              <  2:6246259/93‑1 (MQ=255)

                                                             tCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtg                            >  2:2118584/1‑93 (MQ=255)

                                                              cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  1:1249713/93‑1 (MQ=255)

                                                              cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  2:1277205/93‑1 (MQ=255)

                                                                cTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                         <  1:986790/93‑1 (MQ=255)

                                                                cTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                         <  1:4229994/93‑1 (MQ=255)

                                                                cTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                         >  1:3829746/1‑93 (MQ=255)

                                                                 tGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                        >  1:7628854/1‑93 (MQ=255)

                                                                 tGCCAGCAGGCCAATCACCAACCCGGGCGGCAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                        >  2:3828731/1‑93 (MQ=255)

                                                                  gCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatg                       <  1:274323/93‑1 (MQ=255)

                                                                   ccAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatga                      >  2:297375/1‑93 (MQ=255)

                                                                    cagcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTgatgat                     <  2:7185565/93‑1 (MQ=255)

                                                                      gcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCg                   >  2:938867/1‑93 (MQ=255)

                                                                      gcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCg                   >  2:759143/1‑93 (MQ=255)

                                                                      gcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCg                   <  2:438493/93‑1 (MQ=255)

                                                                      gcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCg                   <  2:5587505/93‑1 (MQ=255)

                                                                      gcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCg                   >  2:1783068/1‑93 (MQ=255)

                                                                         ggCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                <  2:4625625/93‑1 (MQ=255)

                                                                          gCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAat               <  1:886556/93‑1 (MQ=255)

                                                                           ccAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAata              >  2:1603542/1‑93 (MQ=255)

                                                                            cAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat             >  1:5223919/1‑93 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            >  2:4187900/1‑93 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            >  1:2325524/1‑93 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            <  1:2062168/93‑1 (MQ=255)

                                                                              aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTc           >  1:2894095/1‑93 (MQ=255)

                                                                              aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTc           >  1:2443731/1‑93 (MQ=255)

                                                                               tCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCa          <  2:4100301/93‑1 (MQ=255)

                                                                                cACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAt         <  2:7548603/93‑1 (MQ=255)

                                                                                 aCCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        <  1:6778316/93‑1 (MQ=255)

                                                                                  ccAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtct       <  2:6509769/93‑1 (MQ=255)

                                                                                   cAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtctc      <  1:4244497/93‑1 (MQ=255)

                                                                                      cccGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCgg   >  2:1674543/1‑93 (MQ=255)

ATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477006‑477185

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