BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,572,726	T→C	P341P (CCT→CCC)	yjiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,572,726	0	T	C	99.3%	487.6 / ‑4.6	151	P341P (CCT→CCC)	yjiU	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (82/68); total (82/69)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 158 total aligned reads displayed.

CATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCA  >  NC_000913/4572637‑4572815

cATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCaaa                                                                                        <  1:2964881/93‑1 (MQ=255)

cATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCaaa                                                                                        >  2:771307/1‑93 (MQ=255)

cATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCaaa                                                                                        >  1:1259648/1‑93 (MQ=255)

cATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCaaa                                                                                        >  1:1969869/1‑93 (MQ=255)

 aTTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAAc                                                                                       <  2:1719187/93‑1 (MQ=255)

  ttATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACg                                                                                      >  2:2530782/1‑93 (MQ=255)

   tATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACgg                                                                                     <  1:4099681/93‑1 (MQ=255)

    aTCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGa                                                                                    >  1:1938782/1‑93 (MQ=255)

    aTCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGa                                                                                    <  1:1770072/93‑1 (MQ=255)

     tCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGAc                                                                                   >  1:4617390/1‑93 (MQ=255)

     tCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGAc                                                                                   <  2:2239590/93‑1 (MQ=255)

      cTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAGTGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACt                                                                                  >  2:4709432/1‑93 (MQ=255)

       ttAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTa                                                                                 >  1:2506422/1‑93 (MQ=255)

       ttAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTa                                                                                 >  2:2721247/1‑93 (MQ=255)

         aaGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAAc                                                                               <  1:2683922/93‑1 (MQ=255)

           gCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAg                                                                             <  1:435925/93‑1 (MQ=255)

            ccTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGc                                                                            <  1:2002169/93‑1 (MQ=255)

            ccTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGc                                                                            >  2:4449191/1‑93 (MQ=255)

            ccTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGc                                                                            <  2:1675530/93‑1 (MQ=255)

             cTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGcc                                                                           >  1:788244/1‑93 (MQ=255)

             cTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGcc                                                                           <  2:1124782/93‑1 (MQ=255)

               gAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTg                                                                         >  1:311429/1‑93 (MQ=255)

                 aGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGa                                                                       >  2:3870804/1‑93 (MQ=255)

                  gtgtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGaa                                                                      >  1:102583/1‑93 (MQ=255)

                   tgtTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGaaa                                                                     <  1:4430764/93‑1 (MQ=255)

                      ttAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAacac                                                                  <  2:2067293/93‑1 (MQ=255)

                        aaTGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAg                                                                >  2:3090137/1‑93 (MQ=255)

                        aaTGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAg                                                                >  2:754652/1‑93 (MQ=255)

                          tGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAggg                                                              >  1:1586600/1‑93 (MQ=255)

                            aaaaCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGc                                                            <  1:2213113/93‑1 (MQ=255)

                            aaaaCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGc                                                            <  2:3199883/93‑1 (MQ=255)

                            aaaaCGGCGCCCCCCTTAAGCTTTCTGAAATGTTCTTCTTCCGGCAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAACCACAggcaa                                                            <  2:1420662/93‑4 (MQ=255)

                             aaaCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGcg                                                           >  1:1039986/1‑92 (MQ=255)

                             aaaCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCt                                                           >  1:2745483/1‑93 (MQ=255)

                              aaCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTg                                                          >  1:3197403/1‑93 (MQ=255)

                               aCGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTgg                                                         >  1:1791076/1‑93 (MQ=255)

                                cggcggCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGt                                                        <  1:4729437/93‑1 (MQ=255)

                                cggcggCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGt                                                        <  2:588325/93‑1 (MQ=255)

                                 ggcggcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGtt                                                       >  1:1779640/1‑93 (MQ=255)

                                  gcggcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTc                                                      >  1:2932616/1‑93 (MQ=255)

                                   cggcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCa                                                     <  1:422052/93‑1 (MQ=255)

                                   cggcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCa                                                     <  2:3400731/93‑1 (MQ=255)

                                    ggcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAg                                                    <  1:628459/93‑1 (MQ=255)

                                     gcCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAgg                                                   <  1:2882144/93‑1 (MQ=255)

                                      cccccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggt                                                  >  1:3920455/1‑93 (MQ=255)

                                      cccccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggt                                                  <  1:3173133/93‑1 (MQ=255)

                                         ccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttgg                                               <  1:911170/93‑1 (MQ=255)

                                         ccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttgg                                               >  1:578237/1‑93 (MQ=255)

                                         ccTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAgggtgg                                               >  2:3552299/1‑93 (MQ=255)

                                          cTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggt                                              <  2:841443/93‑1 (MQ=255)

                                          cTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggt                                              >  2:4006897/1‑93 (MQ=255)

                                           ttAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggtt                                             <  2:3536178/93‑1 (MQ=255)

                                           ttAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggtt                                             >  1:3874401/1‑93 (MQ=255)

                                           ttAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggttggtt                                             <  2:1863194/93‑1 (MQ=255)

                                           ttAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAggctggtt                                             >  1:4141868/1‑93 (MQ=255)

                                             aaGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTAc                                           <  2:1206483/93‑1 (MQ=255)

                                               gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTTGAAACACAGGGGCTGGTTCAGGTTGGTTACCa                                         >  1:817170/1‑93 (MQ=255)

                                               gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTTCCa                                         >  2:1260170/1‑93 (MQ=255)

                                               gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCa                                         >  2:3780393/1‑93 (MQ=255)

                                               gCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCa                                         <  2:231223/93‑1 (MQ=255)

                                                cTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAg                                        <  2:4122283/93‑1 (MQ=255)

                                                cTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAg                                        >  2:2255957/1‑93 (MQ=255)

                                                 tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgg                                       >  1:142235/1‑93 (MQ=255)

                                                 tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgg                                       >  1:1601992/1‑93 (MQ=255)

                                                 tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgg                                       <  1:383475/93‑1 (MQ=255)

                                                 tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgg                                       >  1:18148/1‑93 (MQ=255)

                                                 tttCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgg                                       <  1:76070/93‑1 (MQ=255)

                                                   tCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAgggg                                     >  1:1463395/1‑93 (MQ=255)

                                                      gAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTg                                  >  2:3357195/1‑93 (MQ=255)

                                                      gAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTg                                  >  2:1483780/1‑93 (MQ=255)

                                                       aaaTGTTGTTGTTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTgg                                 >  2:4562565/1‑93 (MQ=255)

                                                        aaTTTTTTTGTTCCGGGAATTCATGCGGCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGCCTGGTTCAGGTTGGTTACCAGGGGCTGga                                <  1:3566860/93‑1 (MQ=255)

                                                          tgttgttgtTCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGgaga                              <  2:2725083/93‑1 (MQ=255)

                                                           gttgttgttCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGaa                             >  2:2315783/1‑93 (MQ=255)

                                                             tgttgttCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGaaaa                           <  1:2051638/93‑1 (MQ=255)

                                                                 ntnCCGGNAATTCATGCGTCGNCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACa                       <  1:3713016/92‑1 (MQ=255)

                                                                 gttCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACa                       >  1:3473305/1‑93 (MQ=255)

                                                                 gttCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACa                       >  2:2129013/1‑93 (MQ=255)

                                                                   tCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGTTGGAGAAAATACATa                     <  2:4111470/93‑1 (MQ=255)

                                                                   tCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATa                     <  1:4063729/93‑1 (MQ=255)

                                                                   tCCGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATa                     <  1:2792576/93‑1 (MQ=255)

                                                                    ccGGGAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATaa                    <  2:173260/93‑1 (MQ=255)

                                                                      gggAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAaga                  >  2:1569311/1‑93 (MQ=255)

                                                                      gggAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAaga                  <  2:806472/93‑1 (MQ=255)

                                                                       ggAATTCATGCGACGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGGAAATACATAagag                 >  2:1538136/1‑93 (MQ=255)

                                                                        gAATTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGc                >  1:507988/1‑93 (MQ=255)

                                                                        gAAGTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGc                <  2:1537329/93‑1 (MQ=255)

                                                                         aaTTCATGCGTCGCCCCAAACGGACTCACAGCCTGGAAACACAGGGGCTGGTCCAGGTTGGTTCCCAGGGGCTGGCGCAACAACAAAAAAGcc               >  1:1155092/1‑93 (MQ=255)

                                                                         aaTTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGcc               >  1:1155089/1‑93 (MQ=255)

                                                                         aaTTCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGcc               >  1:4560668/1‑93 (MQ=255)

                                                                           ttCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGcccc             >  2:1667605/1‑93 (MQ=255)

                                                                           ttCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGcccc             >  1:4728916/1‑93 (MQ=255)

                                                                            tCATGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCt            <  1:1264386/93‑1 (MQ=255)

                                                                              aTGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTg          >  1:3724382/1‑93 (MQ=255)

                                                                              aTGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTg          <  1:417883/93‑1 (MQ=255)

                                                                               tGCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGc         <  2:4321698/93‑1 (MQ=255)

                                                                                gCGTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCa        <  2:3058376/93‑1 (MQ=255)

                                                                                  gTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACt      >  1:3462524/1‑93 (MQ=255)

                                                                                  gTCGCCCCAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACt      <  1:13121/93‑1 (MQ=255)

                                                                                      ccccAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGACCCTTGCACTGGCa  >  2:885733/1‑93 (MQ=255)

CATTATCTTAAGCCTGAAGTGTTTAATGAAAACGGCGGCCCCCTTAAGCTTTCTGAAATGTTGTTGTTCCGGGAATTCATGCGTCGCCCTAAACGGACTAACAGCCTGGAAACACAGGGGCTGGTTCAGGTTGGTTACCAGGGGCTGGAGAAAATACATAAGAGCCCCTTGCACTGGCA  >  NC_000913/4572637‑4572815

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