BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,804,224	A→G	T157T (ACT→ACC)	rfaG ←	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,804,224	0	A	G	99.7%	903.4 / ‑5.3	290	T157T (ACT→ACC)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (173/116); total (173/117)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 325 total aligned reads displayed.

GCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCGCTCGCTC  >  NC_000913/3804134‑3804313

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGttt                                                                                       <  1:4636797/93‑1 (MQ=255)

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGttt                                                                                       >  2:4148199/1‑93 (MQ=255)

  gATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGATCAGGTTCGGttt                                                                                       >  2:4472086/1‑93 (MQ=255)

   aTAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTg                                                                                      >  2:1143605/1‑93 (MQ=255)

    tAAATTTCACGGGTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGa                                                                                     <  1:6326754/93‑1 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  <  2:7517257/93‑1 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  <  2:6213647/93‑1 (MQ=255)

       aTTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATaa                                                                                  >  2:3780320/1‑93 (MQ=255)

        tttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAAt                                                                                 >  2:5879942/1‑93 (MQ=255)

        tttCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACAGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAAt                                                                                 <  2:6084022/93‑1 (MQ=255)

          tCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGc                                                                               <  1:106833/93‑1 (MQ=255)

            aCGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCtt                                                                             >  2:5306157/1‑93 (MQ=255)

             cGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTc                                                                            <  1:667999/93‑1 (MQ=255)

             cGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTc                                                                            >  2:5495664/1‑93 (MQ=255)

             cGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTc                                                                            <  1:5263915/93‑1 (MQ=255)

             cGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTc                                                                            <  2:726942/93‑1 (MQ=255)

              ggCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCt                                                                           <  2:3684734/93‑1 (MQ=255)

              ggCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCt                                                                           <  2:1778671/93‑1 (MQ=255)

                cTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTgg                                                                         >  1:1541509/1‑93 (MQ=255)

                 tGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGa                                                                        >  1:5389046/1‑93 (MQ=255)

                  gTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATACTGCTTCTGGaa                                                                       >  1:4786126/1‑93 (MQ=255)

                  gTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaa                                                                       >  1:474908/1‑93 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      >  1:7243322/1‑93 (MQ=255)

                   tttGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGaaa                                                                      >  2:3412696/1‑93 (MQ=255)

                     tGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATc                                                                    <  2:3027190/93‑1 (MQ=255)

                     tGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATc                                                                    <  2:2644595/93‑1 (MQ=255)

                        gATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGc                                                                 <  2:7836962/93‑1 (MQ=255)

                         aTTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCg                                                                >  2:6317009/1‑93 (MQ=255)

                         aTTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCg                                                                >  2:1609739/1‑93 (MQ=255)

                          tttGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGa                                                               <  2:4671345/93‑1 (MQ=255)

                            tGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGAtt                                                             <  1:5566896/93‑1 (MQ=255)

                               tCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGc                                                          <  2:4997443/93‑1 (MQ=255)

                                cACTGTGTTTTATGTCCAGATAAATAACGGGAGGAAGAATTTGCAAACGTTCAGGTTCGGTTTGATAATGCTTCTGAAAATCGGCGATTTgat                                                         <  2:7046126/93‑3 (MQ=255)

                                cACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCt                                                         <  1:3481881/93‑1 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        >  1:6387763/1‑93 (MQ=255)

                                 aCTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        >  2:2816608/1‑93 (MQ=255)

                                 aCTGTATTTTCTGTACGGATAACTACCGGGACGAAGAATTTGAAACCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCtt                                                        <  1:7770362/93‑1 (MQ=255)

                                  cTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTa                                                       >  1:1398943/1‑93 (MQ=255)

                                  cTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTa                                                       >  2:4904351/1‑93 (MQ=255)

                                  cTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTa                                                       <  2:1065235/93‑1 (MQ=255)

                                  cTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTa                                                       <  1:3890976/93‑1 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTag                                                      <  2:5617554/93‑2 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTAt                                                      <  2:6873964/93‑1 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTAt                                                      <  1:2703033/93‑1 (MQ=255)

                                   tGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGGTTGATAATGCTTCTGGAAATCGGCGATTTGCTTAt                                                      >  1:3841609/1‑93 (MQ=255)

                                    gTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATc                                                     >  2:5753208/1‑93 (MQ=255)

                                     tATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCg                                                    >  1:2383220/1‑93 (MQ=255)

                                      aTTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCgg                                                   <  1:3213063/93‑1 (MQ=255)

                                       ttttCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGATTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGt                                                  <  2:6647482/93‑1 (MQ=255)

                                          tCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTTCTTATCGGTCAg                                               >  2:1988543/1‑93 (MQ=255)

                                           cTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGc                                              >  2:6444985/1‑93 (MQ=255)

                                            tGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGca                                             >  2:7642341/1‑93 (MQ=255)

                                               ccGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatca                                          <  2:821445/93‑1 (MQ=255)

                                               ccGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatca                                          <  1:6024755/93‑1 (MQ=255)

                                                cGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGcatcat                                         >  2:127689/1‑93 (MQ=255)

                                                  gATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATaa                                       >  1:5138167/1‑93 (MQ=255)

                                                  gATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATaa                                       >  1:5830147/1‑93 (MQ=255)

                                                  gATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATaa                                       >  1:7596910/1‑93 (MQ=255)

                                                    tAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGt                                     >  1:119300/1‑93 (MQ=255)

                                                    tAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGt                                     >  1:7361272/1‑93 (MQ=255)

                                                     aaaTACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGtt                                    >  2:7079148/1‑93 (MQ=255)

                                                         aCCGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgt                                <  2:1464137/93‑1 (MQ=255)

                                                          ccGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtc                               <  2:5346664/93‑1 (MQ=255)

                                                           cGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtcg                              >  2:2553298/1‑93 (MQ=255)

                                                           cGGGAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTtcgtcg                              >  1:6631903/1‑93 (MQ=255)

                                                            gggAGGAAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGa                             >  1:7553801/1‑93 (MQ=255)

                                                             ggaggaATAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAt                            <  1:3677611/93‑1 (MQ=255)

                                                             ggaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAt                            <  2:2362586/93‑1 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           <  1:5538970/93‑1 (MQ=255)

                                                              gaggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAtt                           >  1:5004180/1‑93 (MQ=255)

                                                                ggaAGAATTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTa                         >  2:1899448/1‑93 (MQ=255)

                                                                 gaagaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAc                        >  2:4844253/1‑93 (MQ=255)

                                                                 gaagaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAc                        >  2:246928/1‑93 (MQ=255)

                                                                 gaagaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAc                        >  1:4573490/1‑93 (MQ=255)

                                                                   agaaTTTGAGAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  1:5913293/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  2:6160891/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  1:4433893/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      <  2:6731862/93‑1 (MQ=255)

                                                                   agaaTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTAccc                      >  1:2148097/1‑93 (MQ=255)

                                                                      aTTTGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGc                   >  2:7141039/1‑93 (MQ=255)

                                                                        ttGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTc                 >  2:3592737/1‑93 (MQ=255)

                                                                         tGAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCg                >  1:1871546/1‑93 (MQ=255)

                                                                          gAAAACGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGa               >  2:3194367/1‑93 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGCAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaaa             >  1:3991459/1‑93 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaaa             >  1:6903820/1‑93 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGaaa             <  2:3293672/93‑1 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCATCGATTTACCCTGCTCGaaa             <  2:1571815/93‑1 (MQ=255)

                                                                            aaaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATATTAAGTTTCGTCGATTTACCCTGCTCGaaa             >  1:6730495/1‑93 (MQ=255)

                                                                             aaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAg            <  2:6234414/93‑1 (MQ=255)

                                                                             aaCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAg            >  1:3607789/1‑93 (MQ=255)

                                                                              aCGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGt           >  1:801904/1‑93 (MQ=255)

                                                                               cGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtc          <  1:1934650/93‑1 (MQ=255)

                                                                               cGTTCAGGTTCGGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtc          >  2:3845201/1‑93 (MQ=255)

                                                                                gttcaggttcggTTTTATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcg         >  2:7230551/1‑93 (MQ=255)

                                                                                gttcaggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcg         >  1:773923/1‑93 (MQ=255)

                                                                                   caggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctc      >  2:2148139/1‑93 (MQ=255)

                                                                                   caggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctc      <  2:594517/93‑1 (MQ=255)

                                                                                   caggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGtcgctc      >  1:6557639/1‑93 (MQ=255)

                                                                                   caggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGAGTTACTCTGCTCGAAAGtcgctc      >  1:7044831/1‑93 (MQ=255)

                                                                                     ggttcggTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATGTACCCTGCTCGAAAGtcgctcgc    >  2:4289756/1‑93 (MQ=255)

GCGATAAATTTCACGGCTGTTTGGGATTTGCTCACTGTATTTTCTGTCCGGATAAATACCGGGAGGAAGAATTTGAAAACGTTCAGGTTCAGTTTGATAATGCTTCTGGAAATCGGCGATTTGCTTATCGGTCAGCATCATAAGTTTCGTCGATTTACCCTGCTCGAAAGTCGCTCGCTC  >  NC_000913/3804134‑3804313

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