Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,586,594 (G)6→7 coding (168/3513 nt) hsdR ← type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,586,5881.G93.0% 184.3 / 6.8 58coding (174/3513 nt)hsdRtype I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (27/26);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAG  >  NC_000913/4586449‑4586717
                                                                                                                                              |                                                                                                                                    
gggCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACaa                                                                                                                                       <  2:558547/140‑1 (MQ=255)
      cGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggggg                                                                                                                             >  2:28130/1‑140 (MQ=255)
      cGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggggg                                                                                                                             >  2:161910/1‑140 (MQ=255)
      cGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTCCACGCAGGAGATAGTGTTGATTCTCACAAgggggg                                                                                                                             >  2:280705/1‑140 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                           >  2:114697/1‑138 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGAGCGTGTTGATACTCACAAgggtgggg                                                                                                                           >  1:323538/1‑138 (MQ=255)
          gagaTCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggggggatg                                                                                                                         <  2:354249/140‑5 (MQ=255)
             aTCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTg                                                                                                                       <  2:598072/140‑1 (MQ=255)
                gTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGAAGATCGTCTTGGTTGTCACAAGGGGGGGATGTTGAGt                                                                                                                    >  1:327684/1‑140 (MQ=255)
                 ttATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGATTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTg                                                                                                                       <  2:659771/136‑1 (MQ=255)
                 ttATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTg                                                                                                                       >  1:659771/1‑136 (MQ=255)
                  tATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTaa                                                                                                                  <  1:536503/140‑1 (MQ=255)
                  tATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGCGGGATGTTGAGTaa                                                                                                                  <  2:417895/140‑1 (MQ=255)
                   aTGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAAc                                                                                                                 >  2:420150/1‑140 (MQ=255)
                           cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAg                                                                                                         <  2:391148/140‑1 (MQ=255)
                                aCCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAACGGGGGGATGTTGAGTAACAGACCAAGATGtt                                                                                                    >  2:221111/1‑140 (MQ=255)
                                  cGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTc                                                                                                  >  2:502023/1‑140 (MQ=255)
                                   gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCg                                                                                                 >  1:442007/1‑140 (MQ=255)
                                   gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCg                                                                                                 <  1:415796/140‑1 (MQ=255)
                                   gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCg                                                                                                 >  1:521226/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAGTTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                >  1:388831/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAGAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGTCGGGGTTGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                >  1:462243/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACGAGGGGGGGATGATGAGTAACAGACCAAGATGTTTCGc                                                                                                >  1:182573/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                >  2:113768/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                >  2:358227/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                >  1:83707/1‑140 (MQ=255)
                                          tACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGc                                                                                          <  1:458340/140‑1 (MQ=255)
                                               aTGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCttagc                                                                                     <  2:340326/140‑4 (MQ=255)
                                               aTGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGc                                                                                     <  2:428010/140‑1 (MQ=255)
                                                 gcggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCa                                                                                   >  2:446416/1‑140 (MQ=255)
                                                  cggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaa                                                                                  >  1:529539/1‑140 (MQ=255)
                                                   ggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaaa                                                                                 >  2:69719/1‑140 (MQ=255)
                                                     cgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACa                                                                               >  1:346606/1‑140 (MQ=255)
                                                     cgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAAGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACa                                                                               >  1:94576/1‑140 (MQ=255)
                                                       tAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATa                                                                             >  1:82228/1‑140 (MQ=255)
                                                        aaTTTGTGCAATACAGTAAGGACGTTGTCATCACCGAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATAc                                                                            <  2:321794/140‑1 (MQ=255)
                                                        aaTTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATAc                                                                            <  2:529539/140‑1 (MQ=255)
                                                          tttGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTCCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGc                                                                          <  2:622058/140‑1 (MQ=255)
                                                             gtgAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAtt                                                                       <  2:521226/140‑1 (MQ=255)
                                                             gtgAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCAtt                                                                       >  2:88251/1‑140 (MQ=255)
                                                                  aTACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAAt                                                                  <  1:114697/140‑1 (MQ=255)
                                                                       gagagGATTTTGTCATCAACAAAGGCGATTCTGCCTAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                             <  2:327684/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                             <  2:83707/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                             <  2:346606/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                             <  1:420150/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                             <  2:442007/140‑1 (MQ=255)
                                                                        agagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCCTCGCCAAACATACGCATTTTAATCAGcgt                                                            >  1:229253/1‑140 (MQ=255)
                                                                          agGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtcg                                                          >  1:92478/1‑140 (MQ=255)
                                                                            gATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCgtg                                                        <  1:161910/140‑1 (MQ=255)
                                                                                 tgtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCGTCGCCACACATACGCATTTTAATCAGCGTCGTGTTggg                                                   >  2:461185/1‑140 (MQ=255)
                                                                                  gtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGa                                                  >  1:580451/1‑140 (MQ=255)
                                                                                   tcatcaACAAAGGCGGTTTTGCCGAGTTCCCGCAGGATATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGat                                                 <  1:299749/140‑1 (MQ=255)
                                                                                    catcaACAAAGGCGATTTTGCCGAGTTCACGGAGCAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTAGCTTCGCCAAACATACCCATTTTAATCAACATCGTGTTGGGatc                                                <  2:462243/140‑1 (MQ=255)
                                                                                    catcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGAGGTTGAGTAACAGACCAAGATGTTTGGCTGTGCCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatc                                                >  2:111140/1‑140 (MQ=255)
                                                                                          caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCgtg                                                        >  1:102375/1‑126 (MQ=255)
                                                                                          caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCgtg                                                        <  2:102375/126‑1 (MQ=255)
                                                                                             aGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATATGCATTTTAATCAGCATCGTGTCGGGATCATCCGGGTa                                       >  1:312351/1‑140 (MQ=255)
                                                                                                      tGCCGAGTTCACGCATGAGATCCTGTTGATTCTCACATGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTc                              <  1:221111/140‑1 (MQ=255)
                                                                                                             ttCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGca                       <  1:502023/140‑1 (MQ=255)
                                                                                                                  gCAGGAGATCGTGTTGACTCTCTCAAGGGGGTGATCTTGCGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGc                  <  2:182573/140‑1 (MQ=255)
                                                                                                                    aGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGa                <  1:461185/140‑1 (MQ=255)
                                                                                                                    aGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGa                <  1:159515/140‑1 (MQ=255)
                                                                                                                       agaTCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATgg             <  1:358227/140‑1 (MQ=255)
                                                                                                                               ttGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTg     <  2:439014/140‑1 (MQ=255)
                                                                                                                                tGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGa    >  2:645786/1‑140 (MQ=255)
                                                                                                                                 gATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGaa   >  2:608967/1‑140 (MQ=255)
                                                                                                                                  aTTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAg  >  1:537885/1‑140 (MQ=255)
                                                                                                                                              |                                                                                                                                    
GGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAG  >  NC_000913/4586449‑4586717

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: