BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,959,281	G→A	E330E (GAG→GAA)	rep →	rep helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,959,281	0	G	A	99.5%	799.3 / ‑5.0	218	E330E (GAG→GAA)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (112/104); total (113/105)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

Only 100 of 221 total aligned reads displayed.

GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAA  >  NC_000913/3959191‑3959370

  gTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGa                                                                                       >  1:477603/1‑93 (MQ=255)

   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:2544758/1‑93 (MQ=255)

   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      <  2:3540831/93‑1 (MQ=255)

   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:3514747/1‑93 (MQ=255)

   tCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGAt                                                                                      >  1:2528838/1‑93 (MQ=255)

    cTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATc                                                                                     >  2:4459604/1‑93 (MQ=255)

      gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGc                                                                                   >  1:947240/1‑93 (MQ=255)

      gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGc                                                                                   >  1:4885507/1‑93 (MQ=255)

       ttCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGcc                                                                                  >  1:2928046/1‑93 (MQ=255)

         ctcCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCa                                                                                >  2:1089190/1‑93 (MQ=255)

           ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATc                                                                              >  2:3431735/1‑93 (MQ=255)

             gAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCAc                                                                            >  1:3078695/1‑93 (MQ=255)

               aCTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACtt                                                                          >  1:5105157/1‑93 (MQ=255)

                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         >  2:2513834/1‑93 (MQ=255)

                cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTc                                                                         >  2:566319/1‑93 (MQ=255)

                  gggTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGt                                                                       >  1:3043888/1‑93 (MQ=255)

                    gTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCa                                                                     >  1:5999390/1‑93 (MQ=255)

                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  2:5867940/93‑1 (MQ=255)

                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  1:463899/93‑1 (MQ=255)

                     ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCaa                                                                    <  1:4877585/93‑1 (MQ=255)

                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  <  1:5563972/93‑1 (MQ=255)

                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  >  1:4467734/1‑93 (MQ=255)

                       aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATa                                                                  <  2:5900178/93‑1 (MQ=255)

                          gcgcGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATaaaa                                                               >  1:1037537/1‑93 (MQ=255)

                          gcgcGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATaaaa                                                               <  2:624090/93‑1 (MQ=255)

                            gcgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACg                                                             >  2:319306/1‑93 (MQ=255)

                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            <  2:1764367/93‑1 (MQ=255)

                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            <  2:3624879/93‑1 (MQ=255)

                             cgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGc                                                            >  1:3102720/1‑93 (MQ=255)

                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           <  1:3837008/93‑1 (MQ=255)

                              ggAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCa                                                           >  1:1964978/1‑93 (MQ=255)

                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:2131001/93‑1 (MQ=255)

                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:3978443/93‑1 (MQ=255)

                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          >  1:4751568/1‑93 (MQ=255)

                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  2:5043216/93‑1 (MQ=255)

                               gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAg                                                          <  1:482350/93‑1 (MQ=255)

                                  cTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAc                                                       >  1:773437/1‑93 (MQ=255)

                                  cTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAc                                                       <  2:3444080/93‑1 (MQ=255)

                                  cTAAAAGAATTAAGCGCGAATAACGAAGAACATGAGGCTGAGGGCGTTAATGGCGAACTGATGCCCCATCACTTCGTCAATAAAACGCAgcaa                                                       >  2:800717/1‑90 (MQ=255)

                                     aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaaa                                                    >  1:4390450/1‑93 (MQ=255)

                                     aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaaa                                                    >  2:4042052/1‑93 (MQ=255)

                                     aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACaaa                                                    <  1:5577308/93‑1 (MQ=255)

                                      aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAg                                                   >  1:6102937/1‑93 (MQ=255)

                                      aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATAACTTCGTCAATAAAACGCAGTACaaaa                                                   >  2:1529083/1‑92 (MQ=255)

                                       atcaTTGAGCGCGGATAGCGCAGTGCATGCGTCTGTGCGCGTTACTGGCAAACTGATCGCCCATCACTTCCTCAATAAAACGCAGTACAAAgc                                                  <  2:597851/90‑2 (MQ=255)

                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  <  2:305991/93‑1 (MQ=255)

                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  <  1:4509227/93‑1 (MQ=255)

                                       aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGa                                                  >  2:5285495/1‑93 (MQ=255)

                                        gTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAt                                                 >  1:1059270/1‑93 (MQ=255)

                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACCAAGAtt                                                >  2:1346169/1‑93 (MQ=255)

                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                >  1:677663/1‑93 (MQ=255)

                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                >  2:208659/1‑93 (MQ=255)

                                         tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGAtt                                                <  1:5064109/93‑1 (MQ=255)

                                          aTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTa                                               <  2:931703/93‑1 (MQ=255)

                                          aTTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTa                                               >  1:5234186/1‑93 (MQ=255)

                                           ttAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAc                                              <  2:614465/93‑1 (MQ=255)

                                           ttAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTAc                                              <  1:1930576/93‑1 (MQ=255)

                                              aGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGcc                                           >  1:5790405/1‑93 (MQ=255)

                                                  cgAATAACGACGAACATGAGGCTGACCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  1:2570518/93‑1 (MQ=255)

                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCAttt                                       >  2:4258213/1‑92 (MQ=255)

                                                  cgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTc                                       <  1:4613198/93‑1 (MQ=255)

                                                    aaTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCtt                                     <  1:4132438/93‑1 (MQ=255)

                                                     aTAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCttt                                    <  1:5914247/93‑1 (MQ=255)

                                                      tAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTa                                   <  1:5397594/93‑1 (MQ=255)

                                                       aaCGAAGAACATGAGGCTGAGCGCGTTACTGGCCAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTAt                                  <  1:5935652/93‑1 (MQ=255)

                                                          gaagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgc                               >  2:2656834/1‑93 (MQ=255)

                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              >  2:282422/1‑93 (MQ=255)

                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              <  2:2177350/93‑1 (MQ=255)

                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              <  1:2666532/93‑1 (MQ=255)

                                                           aagaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATcgcg                              <  2:2459970/93‑1 (MQ=255)

                                                            agaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCgg                             >  2:4817365/1‑93 (MQ=255)

                                                            agaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTAACAAGATTACGCCATTCTTTATCGCgg                             >  2:3733857/1‑93 (MQ=255)

                                                             gaaCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGt                            <  2:1451426/93‑1 (MQ=255)

                                                               aCATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTaa                          >  2:2102552/1‑93 (MQ=255)

                                                                cATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAc                         >  2:2382915/1‑93 (MQ=255)

                                                                cATGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAc                         >  1:4093706/1‑93 (MQ=255)

                                                                 aTGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAcc                        <  1:442862/93‑1 (MQ=255)

                                                                 aTGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAAcc                        <  2:3598109/93‑1 (MQ=255)

                                                                  tGAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCa                       >  2:4107612/1‑93 (MQ=255)

                                                                   gAGGCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCAt                      <  2:4214893/93‑1 (MQ=255)

                                                                      gCTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAg                   <  1:1818019/93‑1 (MQ=255)

                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  <  1:4704098/93‑1 (MQ=255)

                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  >  2:5754171/1‑93 (MQ=255)

                                                                       cTGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGt                  <  1:3456434/93‑1 (MQ=255)

                                                                        tGAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTc                 >  2:3124887/1‑93 (MQ=255)

                                                                         gAGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcg                >  1:5055877/1‑93 (MQ=255)

                                                                          aGCGCGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgc               >  1:1231497/1‑93 (MQ=255)

                                                                           gcgcGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTCACCAACAGTcgcg              >  1:1707953/1‑93 (MQ=255)

                                                                           gcgcGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgcg              >  2:5534956/1‑93 (MQ=255)

                                                                           gcgcGTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTcgcg              >  1:1518182/1‑93 (MQ=255)

                                                                            cgcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCgg             <  1:1337902/93‑1 (MQ=255)

                                                                             gcgTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCggg            <  2:3774748/93‑1 (MQ=255)

                                                                               gTTACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtg          >  2:1006375/1‑93 (MQ=255)

                                                                                ttACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtgt         >  2:595764/1‑93 (MQ=255)

                                                                                ttACTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGgtgt         >  2:2245058/1‑93 (MQ=255)

                                                                                  aCTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGttt       <  2:4723633/93‑1 (MQ=255)

                                                                                  aCTGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGttt       >  2:6049859/1‑93 (MQ=255)

                                                                                    tGGCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGa     >  2:4646822/1‑93 (MQ=255)

                                                                                     ggCGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaa    >  1:2288396/1‑93 (MQ=255)

                                                                                       cGAACTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGaaaa  <  1:2086246/93‑1 (MQ=255)

GCGTCTGTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAA  >  NC_000913/3959191‑3959370

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