BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,475	G→A	15.1%	L61L (TTG→TTA)	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,475	0	G	A	15.1%	444.2 / 55.4	177	L61L (TTG→TTA)	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base G (77/72); new base A (11/16); total (88/89)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 196 total aligned reads displayed.

CATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113377‑4113574

cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                    <  2:623023/100‑1 (MQ=255)

cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                    >  2:1004854/1‑100 (MQ=255)

 aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                   <  1:3678582/100‑1 (MQ=255)

 aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                   <  2:1755691/100‑1 (MQ=255)

  tAATGACGCTGACCTGACGTTAATTCATATTGATGATGGCTTGAGCGAGTTGTATCCGCGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  >  2:1333496/1‑100 (MQ=255)

  tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  >  2:588684/1‑100 (MQ=255)

   aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  <  1:3786155/99‑1 (MQ=255)

   aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                                                                                                 >  1:910755/1‑100 (MQ=255)

    aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  <  2:3831246/98‑1 (MQ=255)

    aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                                                                                                <  1:3009629/100‑1 (MQ=255)

        cGCTCACCTGACGTTAATTCATATTGATGATGGCTTCAGAGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGACGATATTCTTCAATTGTTgaagaa                                                                                             <  2:1239115/99‑1 (MQ=255)

        cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAATCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                            <  2:3141986/100‑1 (MQ=255)

         gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                           <  1:2709143/100‑1 (MQ=255)

          cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                          <  2:2480430/100‑2 (MQ=255)

           tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                         <  2:543377/100‑1 (MQ=255)

           tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                         <  1:1729797/100‑3 (MQ=255)

           tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTGCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                         <  1:2703317/100‑1 (MQ=255)

           tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTACCCTGCAACAGAAGATATTCATCAACAATTGAAGAATatc                                                                                         <  2:3771064/100‑3 (MQ=255)

             aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                       <  2:3149115/100‑1 (MQ=255)

              ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                      >  2:2393292/1‑100 (MQ=255)

              ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                      >  2:3925713/1‑100 (MQ=255)

              ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                      <  1:4143226/100‑6 (MQ=38)

               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATGCGTGTATCTACTGCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                     <  2:2935022/100‑7 (MQ=37)

               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                     <  2:762820/100‑1 (MQ=255)

               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                     <  1:1879110/100‑7 (MQ=38)

               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                     <  1:4040163/100‑7 (MQ=38)

                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATatcttctg                                                                                    <  1:3125819/100‑8 (MQ=255)

                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                    >  2:3370773/1‑100 (MQ=255)

                 gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgt                                                                                   <  1:2956699/100‑9 (MQ=37)

                  aCGTTAATTCATCTTGATGATGCATTAAGCGATTTGTATCCGTGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgt                                                                                   <  2:184044/99‑9 (MQ=21)

                   cGTTAGTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                 >  1:4192818/1‑100 (MQ=255)

                   cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                 >  1:1525065/1‑100 (MQ=255)

                       aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAg                                                                             >  1:1019956/1‑100 (MQ=255)

                        aTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                            <  2:3621722/100‑1 (MQ=255)

                          tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATCTTCTTCAATTGTTTAAGAATAAGTCGGATAACCAGCt                                                                           >  2:3657800/1‑99 (MQ=255)

                           cATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGt                                                                         <  1:2505311/100‑1 (MQ=255)

                            atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                        <  1:1326971/100‑1 (MQ=255)

                                 gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                   >  2:2092903/1‑100 (MQ=255)

                                 gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                   >  1:2359259/1‑100 (MQ=255)

                                 gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                   <  2:634019/100‑1 (MQ=255)

                                 gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                   >  1:2740497/1‑100 (MQ=255)

                                  atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                  <  1:3821581/100‑1 (MQ=255)

                                  atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                  <  2:2960649/100‑1 (MQ=255)

                                     atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                               <  1:2884074/100‑1 (MQ=255)

                                     atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                               >  2:318741/1‑100 (MQ=255)

                                        gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                            >  1:2654061/1‑100 (MQ=255)

                                         cTTAAGCGAGTTGTATCCGGGTATCTCCTTCGCTGCACCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGCATGAGAAGCTGTATAATCTTACGaaa                                                           >  2:3697349/1‑100 (MQ=255)

                                          ttAAGCGCGCTGTCTCCGGGTATCTACTTCCTTGCAACAGAACATATTCTTCAATTGTTGATGAATAAGACGGATAACAAGCTGTATAAACCGACg                                                              >  2:3244834/1‑96 (MQ=255)

                                           tAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                         >  2:3740952/1‑100 (MQ=255)

                                             aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                       >  2:3490925/1‑100 (MQ=255)

                                             aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                       >  2:2171053/1‑100 (MQ=255)

                                              gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                      >  2:40179/1‑100 (MQ=255)

                                               cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                     >  1:3060177/1‑100 (MQ=255)

                                                 aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                   <  1:974927/100‑1 (MQ=255)

                                                 aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                   >  1:466318/1‑100 (MQ=255)

                                                    tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                <  1:2157860/100‑1 (MQ=255)

                                                     gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                               >  2:2698521/1‑100 (MQ=255)

                                                       aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                               >  2:3020176/1‑98 (MQ=255)

                                                       aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                             >  2:2109506/1‑100 (MQ=255)

                                                        tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                            >  1:2573210/1‑100 (MQ=255)

                                                          cGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGTATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGACGaa                                          <  2:1401928/100‑1 (MQ=255)

                                                            ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                        >  1:2276750/1‑100 (MQ=255)

                                                            ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                        >  1:3725563/1‑100 (MQ=255)

                                                            ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                        <  1:3046735/100‑1 (MQ=255)

                                                             gTATCTACTTCCCTGCAACAGAAGAGATTCTTCAATTGGTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                        <  2:2868270/99‑1 (MQ=255)

                                                             gTATCTACTTCCCGGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                         <  1:3968097/98‑1 (MQ=255)

                                                              ttctgaCTTCCCTACATCCTAATATATTCTTCAATTGTATCAAAATAATTTGGATAACAAGTTGTATAGACTGACGAAACTTATCCACTGGCCGaaa                                         <  2:530062/92‑2 (MQ=255)

                                                              tATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                      >  2:945729/1‑100 (MQ=255)

                                                              tATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                      >  1:193877/1‑100 (MQ=255)

                                                               aTCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaa                                     <  2:22702/100‑1 (MQ=255)

                                                                  tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                  >  1:754734/1‑100 (MQ=255)

                                                                    cTTCCCTGCAACAGTAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                >  2:1527981/1‑100 (MQ=255)

                                                                    cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                <  1:1473545/100‑1 (MQ=255)

                                                                     ttCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                               <  2:3003905/100‑1 (MQ=255)

                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              >  2:3125532/1‑100 (MQ=255)

                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              >  1:3347632/1‑100 (MQ=255)

                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              <  2:2407255/100‑1 (MQ=255)

                                                                       cccTGCAACAGAAGATATTCTTCAATTGGTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                             >  1:1263999/1‑100 (MQ=255)

                                                                        ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                            >  2:120995/1‑100 (MQ=255)

                                                                        ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                            <  2:1907563/100‑1 (MQ=255)

                                                                         cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           >  1:1035555/1‑100 (MQ=255)

                                                                         cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           >  2:4132075/1‑100 (MQ=255)

                                                                         cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           <  1:3490942/100‑1 (MQ=255)

                                                                          tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                          >  1:2234719/1‑100 (MQ=255)

                                                                            cgatctgacttatTCTTCAACAATGGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGa                        <  1:2442568/90‑1 (MQ=255)

                                                                                ccgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                    >  1:1048972/7‑100 (MQ=38)

                                                                                 gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                   <  1:546136/100‑1 (MQ=255)

                                                                                 gaagaTATTCTTCAAGTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGTCGAAAAATATTCAATAGCCGAAGACAAAACTGCGTATTGAAcgcg                   >  2:561178/1‑100 (MQ=255)

                                                                                 cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTTAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                   >  2:338448/6‑100 (MQ=37)

                                                                                  gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  2:2758659/5‑100 (MQ=38)

                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  <  1:3312647/100‑1 (MQ=255)

                                                                                   agaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 <  2:263984/100‑1 (MQ=255)

                                                                                    cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                >  1:2928318/3‑100 (MQ=255)

                                                                                          cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGc          >  1:388946/1‑100 (MQ=255)

                                                                                           ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGTGTATTGAACGCGGAGAAATGcc         <  1:846000/100‑1 (MQ=255)

                                                                                           ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc         >  2:3942503/1‑100 (MQ=255)

                                                                                              aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa      <  1:3740973/100‑1 (MQ=255)

                                                                                              aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa      <  2:1313657/100‑1 (MQ=255)

                                                                                                  gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAaca  >  2:430354/1‑100 (MQ=255)

                                                                                                  gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAaca  >  2:910526/1‑100 (MQ=255)

CATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113377‑4113574

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