BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,461,706	T→C	S225S (AGT→AGC)	b2345 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,461,706	0	T	C	99.3%	435.5 / ‑4.9	139	S225S (AGT→AGC)	b2345	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (81/57); total (81/58)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 143 total aligned reads displayed.

TTACTTATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGT  >  NC_000913/2461616‑2461795

ttACTTATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGa                                                                                         >  1:3000977/1‑93 (MQ=255)

ttACTTATCTCTAGGGTACTACGATCAGGACGACTACCACGTTCAAAGAGGTATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGGGa                                                                                         >  2:4313929/1‑93 (MQ=255)

   cTTATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGaaaa                                                                                      >  1:2450491/1‑93 (MQ=255)

   cTTATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGaaaa                                                                                      <  1:3329035/93‑1 (MQ=255)

   cTTATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGaaaa                                                                                      >  1:711148/1‑93 (MQ=255)

    ttATCTCTCGGTTACTACGATCAGGACGACTACCACATTTAAAGAGGAATTCGTATTAATGGGGATAGTCTCACCCAATATTGTAGCGAAAAc                                                                                     >  2:3607447/1‑93 (MQ=255)

    ttATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAAc                                                                                     >  1:1268366/1‑93 (MQ=255)

     tATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACg                                                                                    >  2:126336/1‑93 (MQ=255)

      aTCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGc                                                                                   >  2:779142/1‑93 (MQ=255)

         tctcGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAg                                                                                <  1:4069750/93‑1 (MQ=255)

         tctcGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAg                                                                                >  2:515913/1‑93 (MQ=255)

         tctcGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAg                                                                                >  2:592448/1‑93 (MQ=255)

           tcGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGa                                                                              <  2:1407221/93‑1 (MQ=255)

             ggTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGt                                                                            <  1:3974915/93‑1 (MQ=255)

              gTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTg                                                                           >  1:1794959/1‑93 (MQ=255)

               ttACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGc                                                                          <  2:4184499/93‑1 (MQ=255)

                tactacGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCt                                                                         >  1:1115145/1‑93 (MQ=255)

                 actacGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCtt                                                                        >  2:1287227/1‑93 (MQ=255)

                   tacGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCa                                                                      <  2:1772188/93‑1 (MQ=255)

                    acGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAg                                                                     <  1:1680251/93‑1 (MQ=255)

                     cGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGc                                                                    >  1:4642838/1‑93 (MQ=255)

                        tCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTg                                                                 <  2:1961404/93‑1 (MQ=255)

                        tCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTg                                                                 >  2:339308/1‑93 (MQ=255)

                        tCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTg                                                                 >  2:1570775/1‑93 (MQ=255)

                         cAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTgg                                                                <  2:4233672/93‑1 (MQ=255)

                          aGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGt                                                               >  2:1472109/1‑93 (MQ=255)

                           ggACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGtt                                                              >  2:2305443/1‑93 (MQ=255)

                              cgacTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGa                                                           <  2:3532939/93‑1 (MQ=255)

                              cgacTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGa                                                           >  1:2090757/1‑93 (MQ=255)

                                acTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGaaa                                                         >  2:785778/1‑93 (MQ=255)

                                  tACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGc                                                       >  1:3525162/1‑93 (MQ=255)

                                   aCCACATTCAAAGAGGAATTCGTATTAATGGTGATGGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCa                                                      <  1:4344081/93‑1 (MQ=255)

                                   aCCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCa                                                      <  2:3913419/93‑1 (MQ=255)

                                    ccACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAg                                                     <  1:701760/93‑1 (MQ=255)

                                     cacaTTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGt                                                    >  2:2980366/1‑93 (MQ=255)

                                     cacaTTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGt                                                    >  1:1469102/1‑93 (MQ=255)

                                      acaTTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTa                                                   >  1:1658392/1‑93 (MQ=255)

                                       caTTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTaa                                                  >  1:277802/1‑93 (MQ=255)

                                       caTTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTaa                                                  >  2:71961/1‑93 (MQ=255)

                                       caTTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTaa                                                  <  1:549217/93‑1 (MQ=255)

                                         ttCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAg                                                >  1:374106/1‑93 (MQ=255)

                                         ttCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAg                                                <  2:3205069/93‑1 (MQ=255)

                                         ttCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAg                                                >  1:84120/1‑93 (MQ=255)

                                            aaaGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTa                                             <  2:1773088/93‑1 (MQ=255)

                                            aaaGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTa                                             >  1:3204358/1‑93 (MQ=255)

                                            aaaGAGGAATTCGTATTAAAGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTa                                             >  2:4552048/1‑93 (MQ=255)

                                             aaGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTAt                                            >  1:1756121/1‑93 (MQ=255)

                                              agagGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATc                                           >  2:3277451/1‑93 (MQ=255)

                                              agagGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATc                                           >  1:4370163/1‑93 (MQ=255)

                                                agGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCAt                                         >  2:579304/1‑93 (MQ=255)

                                                 ggAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATg                                        >  1:1531185/1‑93 (MQ=255)

                                                 ggAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATg                                        <  1:4525465/93‑1 (MQ=255)

                                                   aaTTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGc                                      >  1:2799015/1‑93 (MQ=255)

                                                   aaTTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGc                                      <  2:3459158/93‑1 (MQ=255)

                                                     ttCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAg                                    >  1:244360/1‑93 (MQ=255)

                                                     ttCGTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAg                                    <  1:862121/93‑1 (MQ=255)

                                                        gTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAAt                                 <  2:695681/93‑1 (MQ=255)

                                                        gTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAAt                                 >  2:2175558/1‑93 (MQ=255)

                                                        gTATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAAt                                 >  1:1434365/1‑93 (MQ=255)

                                                         tATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATc                                >  2:3389699/1‑93 (MQ=255)

                                                         tATTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATc                                >  2:4010877/1‑93 (MQ=255)

                                                          aTTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCa                               >  2:3853118/1‑93 (MQ=255)

                                                          aTTAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCa                               <  2:527765/93‑1 (MQ=255)

                                                           ttAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCAt                              >  2:4314367/1‑93 (MQ=255)

                                                           ttAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCAt                              <  1:1827051/93‑1 (MQ=255)

                                                           ttAATGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCAt                              <  2:1596226/93‑1 (MQ=255)

                                                             aaTGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCAttt                            <  1:2911537/93‑1 (MQ=255)

                                                             aaTGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCAttt                            <  1:2325611/93‑1 (MQ=255)

                                                              aTGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTg                           <  1:3840176/93‑1 (MQ=255)

                                                              aTGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTg                           <  1:4045449/93‑1 (MQ=255)

                                                              aTGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTg                           >  1:4553361/1‑93 (MQ=255)

                                                               tGGTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGc                          >  2:4154065/1‑93 (MQ=255)

                                                                ggTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCa                         >  2:4202983/1‑93 (MQ=255)

                                                                 gTGATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATAATTTGCaa                        <  2:4579119/93‑1 (MQ=255)

                                                                   gATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACt                      <  1:3960136/93‑1 (MQ=255)

                                                                   gATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACt                      >  2:1022213/1‑93 (MQ=255)

                                                                   gATAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACt                      <  2:1201270/93‑1 (MQ=255)

                                                                    aTAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTg                     >  1:1089446/1‑93 (MQ=255)

                                                                    aTAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTg                     <  1:3243493/93‑1 (MQ=255)

                                                                    aTAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTg                     >  1:3726653/1‑93 (MQ=255)

                                                                     tAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTgg                    >  1:2899170/1‑93 (MQ=255)

                                                                     tAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTgg                    <  2:2645050/93‑1 (MQ=255)

                                                                     tAGTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTgg                    >  1:2583225/1‑93 (MQ=255)

                                                                       gTCTCACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGtt                  >  2:2796468/1‑93 (MQ=255)

                                                                        tctcACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTc                 <  1:352649/93‑1 (MQ=255)

                                                                         ctcACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTcc                <  1:2049231/93‑1 (MQ=255)

                                                                          tcACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCg               >  2:2420970/1‑93 (MQ=255)

                                                                          tcACCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCg               >  2:3666187/1‑93 (MQ=255)

                                                                            aCCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGat             >  1:996222/1‑93 (MQ=255)

                                                                            aCCCAATATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGat             <  1:3041878/93‑1 (MQ=255)

                                                                                aaTATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGatcatc         <  2:1736393/93‑1 (MQ=255)

                                                                                aaTATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGatcatc         >  2:2610815/1‑93 (MQ=255)

                                                                                aaTATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGatcatc         <  1:4368543/93‑1 (MQ=255)

                                                                                aaTATTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGatcatc         <  2:3652381/93‑1 (MQ=255)

                                                                                 atatTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGAtcatca        >  2:730395/1‑93 (MQ=255)

                                                                                  tatTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAg       <  1:544522/93‑1 (MQ=255)

                                                                                   atTGTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAgg      >  2:3236054/1‑93 (MQ=255)

                                                                                      gTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTg   <  2:2597368/93‑1 (MQ=255)

                                                                                      gTAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTg   >  2:2548971/1‑93 (MQ=255)

                                                                                       tAGCGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGt  >  2:3189864/1‑93 (MQ=255)

TTACTTATCTCTCGGTTACTACGATCAGGACGACTACCACATTCAAAGAGGAATTCGTATTAATGGTGATAGTCTCACCCAATATTGTAGTGAAAACGCCAGGAGTGCTTCAGCGTGGTTTGAAAGCAGTAAAGCTATCATGGCAGAATCATTTGCAACTGGTTCCGATCATCAGGTTGT  >  NC_000913/2461616‑2461795

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