BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,442,106	A→G	S28S (AGT→AGC)	b2325 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,442,106	0	A	G	99.6%	733.8 / ‑5.1	227	S28S (AGT→AGC)	b2325	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (106/120); total (106/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 232 total aligned reads displayed.

GCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCACTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTT  >  NC_000913/2442016‑2442195

gCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCt                                                                                         <  1:608286/93‑1 (MQ=255)

 cTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTg                                                                                        <  2:216684/93‑1 (MQ=255)

 cTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTg                                                                                        >  1:6404083/1‑93 (MQ=255)

  tAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTgg                                                                                       >  1:6647561/1‑93 (MQ=255)

   aaCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGc                                                                                      >  2:4226681/1‑93 (MQ=255)

   aaCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGc                                                                                      >  1:2267249/1‑93 (MQ=255)

   aaCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGc                                                                                      >  1:405192/1‑93 (MQ=255)

   aaCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGc                                                                                      >  1:3748845/1‑93 (MQ=255)

    aCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCa                                                                                     >  2:2607673/1‑93 (MQ=255)

    aCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCa                                                                                     >  2:5712956/1‑93 (MQ=255)

    aCGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCa                                                                                     >  1:2086944/1‑93 (MQ=255)

     cGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCaa                                                                                    <  1:3243472/93‑1 (MQ=255)

         cGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaa                                                                                <  1:4239494/93‑1 (MQ=255)

         cGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaa                                                                                <  1:2896562/93‑1 (MQ=255)

          gTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaac                                                                               >  1:1137332/1‑93 (MQ=255)

          gTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAaacaac                                                                               <  2:4447748/93‑1 (MQ=255)

           tATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACt                                                                              <  2:4133523/93‑1 (MQ=255)

                aCCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatca                                                                         >  2:6662553/1‑93 (MQ=255)

                aCCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatca                                                                         <  2:5879566/93‑1 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        <  2:1675093/93‑1 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  2:2017317/1‑93 (MQ=255)

                 ccGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACtcatcat                                                                        >  2:6091681/1‑93 (MQ=255)

                    ccTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCAc                                                                     <  2:1776580/93‑1 (MQ=255)

                    ccTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCAc                                                                     <  2:3532808/93‑1 (MQ=255)

                    ccTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCAc                                                                     >  2:3410332/1‑93 (MQ=255)

                      tGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTg                                                                   >  1:4498889/1‑93 (MQ=255)

                      tGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTg                                                                   <  2:3149502/93‑1 (MQ=255)

                       gAGCGGA‑TTGTCATTCCCACTTTCCAAGTCTGCCATTGCTAATTTCAGGTTGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGc                                                                 <  1:1554517/93‑1 (MQ=255)

                          cGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCAt                                                               >  2:4385761/1‑93 (MQ=255)

                          cGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCAt                                                               <  2:496636/93‑1 (MQ=255)

                           ggAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATg                                                              >  1:373726/1‑93 (MQ=255)

                           ggAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATg                                                              >  1:3225532/1‑93 (MQ=255)

                             aGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGt                                                            <  1:3482740/93‑1 (MQ=255)

                              gTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTc                                                           >  2:4740865/1‑93 (MQ=255)

                              gTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTc                                                           >  1:5799545/1‑93 (MQ=255)

                               tGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCa                                                          >  2:414043/1‑93 (MQ=255)

                                 gTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAAc                                                        <  1:3565109/93‑1 (MQ=255)

                                  tcatcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCCCGCAAATACCCGCTGGCAAACAACTCATCATCCCTGGCATGGTCAAcc                                                       >  1:4833258/1‑93 (MQ=255)

                                  tcatcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAAcc                                                       <  1:3710855/93‑1 (MQ=255)

                                   catcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATTCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCa                                                      <  2:3956902/93‑1 (MQ=255)

                                    atcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCAt                                                     >  1:6458900/1‑93 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    >  1:5769307/1‑93 (MQ=255)

                                     tcaCCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATa                                                    <  1:3009470/93‑1 (MQ=255)

                                      caccacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATAc                                                   <  1:5000650/93‑1 (MQ=255)

                                      caccacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATAc                                                   <  1:5399863/93‑1 (MQ=255)

                                       accacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATAcc                                                  >  1:1414643/1‑93 (MQ=255)

                                        ccacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCg                                                 <  1:4581775/93‑1 (MQ=255)

                                         cacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGt                                                <  2:2062324/93‑1 (MQ=255)

                                         cacTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGt                                                >  2:1232277/1‑93 (MQ=255)

                                          acTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTc                                               >  2:1385751/1‑93 (MQ=255)

                                          acTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTc                                               >  2:1149066/1‑93 (MQ=255)

                                           cTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCg                                              >  1:6685479/1‑93 (MQ=255)

                                            tttCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGCCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGa                                             <  2:4695680/93‑1 (MQ=255)

                                              tCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGAtt                                           <  2:3317372/93‑1 (MQ=255)

                                               ccAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTa                                          >  1:654761/1‑93 (MQ=255)

                                                 aGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAAt                                        <  2:2411391/93‑1 (MQ=255)

                                                   ttCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGc                                      >  1:3983310/1‑93 (MQ=255)

                                                     cTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCa                                    <  1:4314256/93‑1 (MQ=255)

                                                      tGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAg                                   >  1:5681595/1‑93 (MQ=255)

                                                        cGATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGaa                                 <  2:1871856/93‑1 (MQ=255)

                                                         gATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAAt                                <  1:2716843/93‑1 (MQ=255)

                                                         gATGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAAt                                <  2:2999733/93‑1 (MQ=255)

                                                           tGGCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGc                              >  2:788678/1‑93 (MQ=255)

                                                            ggCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCg                             <  1:1771174/93‑1 (MQ=255)

                                                             gCTAATGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGt                            <  2:856843/93‑1 (MQ=255)

                                                             gCTAATGTCAGGTAGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGt                            <  1:973276/93‑1 (MQ=255)

                                                                aaTGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGAt                         <  2:4460933/93‑1 (MQ=255)

                                                                 aTGTCAGGTGGCCCCGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGAtt                        <  2:4849453/93‑1 (MQ=255)

                                                                 aTGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGAtt                        >  1:183463/1‑93 (MQ=255)

                                                                  tGTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTc                       >  1:1787796/1‑93 (MQ=255)

                                                                   gTCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCa                      >  1:1027065/1‑93 (MQ=255)

                                                                    tCAGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCaa                     <  2:2098578/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   >  2:6514060/1‑93 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   <  2:6433948/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   >  1:5432517/1‑93 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   <  2:3884154/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   <  2:4103108/93‑1 (MQ=255)

                                                                      aGGTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAAc                   <  2:724101/93‑1 (MQ=255)

                                                                       ggTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACt                  >  2:1062055/1‑93 (MQ=255)

                                                                       ggTGGCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACt                  >  1:4883331/1‑93 (MQ=255)

                                                                          ggCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCcg               <  2:6025160/93‑1 (MQ=255)

                                                                           gCCACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCcgc              <  1:6314200/93‑1 (MQ=255)

                                                                            ccACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCcgcg             <  1:2053001/93‑1 (MQ=255)

                                                                             cACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGa            >  2:5588289/1‑93 (MQ=255)

                                                                             cACGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGa            <  2:6212325/93‑1 (MQ=255)

                                                                              aCGCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGat           <  2:370271/93‑1 (MQ=255)

                                                                                gCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         <  1:4272759/93‑1 (MQ=255)

                                                                                gCAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatca         <  1:5414105/93‑1 (MQ=255)

                                                                                 cAAATACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGatcat        <  1:4588608/93‑1 (MQ=255)

                                                                                   aaTACCCGCTGTCAAACACCTCATCATCCCTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtct      <  1:6637106/93‑1 (MQ=255)

                                                                                    aTACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtctt     <  2:5294810/93‑1 (MQ=255)

                                                                                    aTACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtctt     >  1:5381307/1‑93 (MQ=255)

                                                                                     tACCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttt    >  1:4406451/1‑92 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   >  1:865789/1‑93 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   <  2:3141040/93‑1 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   >  1:5196385/1‑93 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   >  1:2513614/1‑93 (MQ=255)

                                                                                      aCCCGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCAtcttct   <  1:4374812/93‑1 (MQ=255)

                                                                                       cccGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATcttctt  <  2:506932/93‑1 (MQ=255)

                                                                                       cccGCTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCAGCGATCATcttatt  >  1:2582139/1‑93 (MQ=255)

GCTAACGGTCGTATGCACCGCCTGAGCGGAGTGGTCATCACCACTTTCCAGTTCTGCGATGGCTAATGTCAGGTGGCCACGCAAATACCCACTGGCAAACAACTCATCATCACTGGCATGGTCAACCATACCGTCGATTAATGCCAGAATGCGTGATTCAAACTCCGCGATCATCTTCTT  >  NC_000913/2442016‑2442195

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