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.G98.2% 405.7 / ‑0.3 109coding (174/3513 nt)hsdRtype I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (64/43);  total (66/43)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 131 total aligned reads displayed.

GGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTG  >  NC_000913/4586449‑4586723
                                                                                                                                             |                                                                                                                                       
gggCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACaa                                                                                                                                          <  1:500540/140‑1 (MQ=255)
gggCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACaa                                                                                                                                          <  1:456689/140‑1 (MQ=255)
  gCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgg                                                                                                                                       >  1:1242341/1‑140 (MQ=255)
    aTCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAgggg                                                                                                                                     <  2:879828/140‑1 (MQ=255)
     tCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggg                                                                                                                                    <  2:1607/140‑1 (MQ=255)
       gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggg                                                                                                                                  <  1:936429/140‑2 (MQ=255)
       gTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggg                                                                                                                                  <  1:820756/140‑2 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 >  1:679507/1‑138 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 >  2:793549/1‑138 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 <  2:423372/140‑3 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 <  2:224853/140‑3 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 <  1:1391546/140‑3 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 <  2:541274/140‑3 (MQ=255)
        ttGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAggggggga                                                                                                                                 <  2:58740/140‑3 (MQ=255)
            gaTCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGtt                                                                                                                              <  1:845218/140‑1 (MQ=255)
            gaTCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGtt                                                                                                                              <  2:161463/140‑1 (MQ=255)
                gTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGt                                                                                                                          >  1:554560/1‑140 (MQ=255)
                 ttATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTa                                                                                                                         >  1:68108/1‑140 (MQ=255)
                  tATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTaa                                                                                                                        <  1:985225/140‑1 (MQ=255)
                  tATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTaa                                                                                                                        >  2:1492263/1‑140 (MQ=255)
                     gATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAg                                                                                                                     >  2:123455/1‑140 (MQ=255)
                        atTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGAcc                                                                                                                  >  2:441842/1‑140 (MQ=255)
                         ttCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCa                                                                                                                 >  1:1187896/1‑140 (MQ=255)
                           cGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGa                                                                                                              >  2:1191255/1‑140 (MQ=255)
                            gTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGa                                                                                                              >  2:571015/1‑140 (MQ=255)
                               cACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGt                                                                                                           >  1:52259/1‑140 (MQ=255)
                                 ccGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGttt                                                                                                         >  1:387119/1‑140 (MQ=255)
                                 ccGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGttt                                                                                                         >  2:1353462/1‑140 (MQ=255)
                                   gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCg                                                                                                       >  1:626501/1‑140 (MQ=255)
                                   gCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCg                                                                                                       >  2:742818/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                      >  1:563969/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                      >  2:1322951/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                      >  2:1301926/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                      >  1:834189/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                      >  1:21733/1‑140 (MQ=255)
                                    ccTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGc                                                                                                      >  1:307702/1‑140 (MQ=255)
                                      tGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCtg                                                                                                    >  2:2548/1‑140 (MQ=255)
                                          tACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGc                                                                                                <  2:1289708/140‑1 (MQ=255)
                                          tACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGc                                                                                                <  2:535528/140‑1 (MQ=255)
                                               aTGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTTATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGc                                                                                           <  2:1023465/140‑1 (MQ=255)
                                               aTGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGc                                                                                           >  2:463309/1‑140 (MQ=255)
                                                tGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGcc                                                                                          >  2:1481933/1‑140 (MQ=255)
                                                tGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGcc                                                                                          >  1:1158279/1‑140 (MQ=255)
                                                  cggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaa                                                                                        >  1:1243914/1‑140 (MQ=255)
                                                  cggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaa                                                                                        >  2:1304509/1‑140 (MQ=255)
                                                   ggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaaa                                                                                       >  1:270348/1‑140 (MQ=255)
                                                   ggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaaa                                                                                       >  2:1077032/1‑140 (MQ=255)
                                                   ggcgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCaaa                                                                                       <  1:1104114/140‑1 (MQ=255)
                                                     cgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACa                                                                                     >  2:743110/1‑140 (MQ=255)
                                                     cgTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACa                                                                                     >  2:610732/1‑140 (MQ=255)
                                                                    aCAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGGTGCTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCa                                                                      >  2:185994/1‑140 (MQ=255)
                                                                    aCAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCa                                                                      <  2:1187896/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  2:1401272/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  2:694122/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  2:368480/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  2:113665/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  1:1304509/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  1:1329956/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  1:1353462/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  2:879149/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  1:323291/140‑1 (MQ=255)
                                                                       gagagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcg                                                                   <  1:954284/140‑1 (MQ=255)
                                                                        agagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgt                                                                  <  2:759094/140‑1 (MQ=255)
                                                                         gagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtc                                                                 <  2:1448390/140‑1 (MQ=255)
                                                                         gagGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGcgtc                                                                 <  1:743110/140‑1 (MQ=255)
                                                                                 tgtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTggg                                                         >  2:56964/1‑140 (MQ=255)
                                                                                  gtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACaa                                                                                                                                          <  2:104366/58‑1 (MQ=255)
                                                                                  gtCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACaa                                                                                                                                          >  1:104366/1‑58 (MQ=255)
                                                                                   tcatcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGat                                                       <  1:763082/140‑1 (MQ=255)
                                                                                      tcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatcat                                                    >  1:1462713/1‑140 (MQ=255)
                                                                                      tcaACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGatcat                                                    <  1:1481933/140‑1 (MQ=255)
                                                                                          caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCgg                                                >  2:1308158/1‑140 (MQ=255)
                                                                                          caaAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCgg                                                >  2:1375636/1‑140 (MQ=255)
                                                                                           aaaGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCggg                                               <  1:441842/140‑1 (MQ=255)
                                                                                             aGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACATGGGGGGGGAGTTGAGTAACAGACCAAGATGATTCGCTGCGGCTTCGCAAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTa                                             >  2:1346818/1‑140 (MQ=255)
                                                                                             aGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTa                                             >  2:284431/1‑140 (MQ=255)
                                                                                                 gATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTa                                         <  2:679507/140‑1 (MQ=255)
                                                                                                  aTTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTAt                                        >  2:1407827/1‑140 (MQ=255)
                                                                                                        ccGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTccg                                  <  1:571015/140‑1 (MQ=255)
                                                                                                        ccGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTccg                                  <  2:821908/140‑1 (MQ=255)
                                                                                                             ttCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGca                             >  1:914748/1‑140 (MQ=255)
                                                                                                             ttCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGca                             <  2:496737/140‑1 (MQ=255)
                                                                                                             ttCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGca                             <  1:1193064/140‑1 (MQ=255)
                                                                                                             ttCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGca                             <  1:56964/140‑1 (MQ=255)
                                                                                                               cACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGcaca                           <  1:1492263/140‑1 (MQ=255)
                                                                                                                   cAGGAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCg                       <  2:270348/140‑1 (MQ=255)
                                                                                                                     ggAGATCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGAt                     >  2:77760/1‑140 (MQ=255)
                                                                                                                      gagaTCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATg                    >  2:453008/1‑140 (MQ=255)
                                                                                                                       agaTCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATgg                   >  1:300298/1‑140 (MQ=255)
                                                                                                                        gaTCGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGc                  <  2:307702/140‑1 (MQ=255)
                                                                                                                           cGTGTTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATa               >  2:918262/1‑140 (MQ=255)
                                                                                                                              gtTGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGt            >  1:1263298/1‑140 (MQ=255)
                                                                                                                               ttGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTg           <  2:1221710/140‑1 (MQ=255)
                                                                                                                                tGATTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGa          >  1:690684/1‑140 (MQ=255)
                                                                                                                                  aTTCTCACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAg        >  2:1414321/1‑140 (MQ=255)
                                                                                                                                    tctcACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTc      >  1:228103/1‑140 (MQ=255)
                                                                                                                                    tctcACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTc      <  1:1077032/140‑1 (MQ=255)
                                                                                                                                     ctcACAAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCg     >  1:1081544/1‑140 (MQ=255)
                                                                                                                                       cacaAGGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGtt   >  2:645071/1‑140 (MQ=255)
                                                                                                                                       cacaAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTg  >  1:736526/1‑140 (MQ=255)
                                                                                                                                             |                                                                                                                                       
GGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGAGGATGTTGTCATCAACAAAGGCGATTTTGCCGAGTTCACGCAGGAGATCGTGTTGATTCTCACAAGGGGGGATGTTGAGTAACAGACCAAGATGTTTCGCTGTGGCTTCGCCAAACATACGCATTTTAATCAGCGTCGTGTTGGGATCATCCGGGTAGTTATTTTCCGCCGCACAGGCGATGGCATAAGTGAAGTCGTTG  >  NC_000913/4586449‑4586723

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