Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,099,249 C→T 5.8% G122D (GGT→GAT)  ugd ← UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,099,2490CT 5.8% 432.2 / 11.0 138G122D (GGT→GAT) ugdUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (63/66);  new base T (6/2);  total (70/68)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 141 total aligned reads displayed.

CAGCGAAACGTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCAGCATCC  >  NC_000913/2099100‑2099398
                                                                                                                                                      |                                                                                                                                                     
cagcGAAACGTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAAc                                                                                                                                                       <  1:1136493/150‑1 (MQ=255)
  gcGAAACGTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCa                                                                                                                                                     <  1:689579/150‑1 (MQ=255)
     aaaCGTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACg                                                                                                                                                  >  1:778756/1‑150 (MQ=255)
        cGTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGa                                                                                                                                               <  1:101380/150‑1 (MQ=255)
         gTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGaa                                                                                                                                              >  2:496343/1‑150 (MQ=255)
              gcgcGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGtt                                                                                                                                         >  1:1170779/1‑150 (MQ=255)
              gcgcGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGtt                                                                                                                                         >  2:1403289/1‑150 (MQ=255)
                 cgTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGttgat                                                                                                                                      <  2:421837/150‑1 (MQ=255)
                  gTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGttgatt                                                                                                                                     <  1:1352372/150‑1 (MQ=255)
                  gTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGttgagt                                                                                                                                     >  2:1429332/1‑148 (MQ=255)
                   ttCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGttgattt                                                                                                                                    <  2:1500634/150‑1 (MQ=255)
                    tCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGttgatttg                                                                                                                                   >  1:406679/1‑150 (MQ=255)
                    tCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGttgatttg                                                                                                                                   >  2:1336410/1‑150 (MQ=255)
                      tGAACGCTCCCCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACACTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATTCATCGCTGCGGTAAAACCAACGGGAAACGttgatatgat                                                                                                                                 >  2:1404128/1‑150 (MQ=255)
                       gAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTtgatg                                                                                                                                <  2:179773/150‑1 (MQ=255)
                        aaCGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAATCAACGGTTCCCGttggttttaccg                                                                                                                               <  1:373454/150‑6 (MQ=255)
                        aaCGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAATCAACGGTTCCCGttggttttaccg                                                                                                                               <  1:1467480/150‑6 (MQ=255)
                        aaCGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTACAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTtgatga                                                                                                                               <  1:1617407/150‑1 (MQ=255)
                          cGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGAcc                                                                                                                             >  2:537575/1‑150 (MQ=255)
                             tCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATa                                                                                                                          <  1:452985/150‑1 (MQ=255)
                              cACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATaa                                                                                                                         >  1:625844/1‑150 (MQ=255)
                                 cGATGACAATAGGTGAAGGATGGAGATTATAGTAAAGGGCTTTACCATCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCg                                                                                                                      >  2:906964/1‑150 (MQ=255)
                                 cGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCg                                                                                                                      <  2:403852/150‑1 (MQ=255)
                                   aTGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCa                                                                                                                    <  2:103983/150‑1 (MQ=255)
                                      aCAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATaa                                                                                                                 >  2:277291/1‑150 (MQ=255)
                                      aCAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATaa                                                                                                                 >  2:50117/1‑150 (MQ=255)
                                      aCAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCAAGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAAACAACGGGAACCGTTGATTTGATGACCATAACCGCATaa                                                                                                                 >  1:464671/1‑150 (MQ=255)
                                       cAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAg                                                                                                                >  2:145143/1‑150 (MQ=255)
                                       cAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACAGCATAAg                                                                                                                >  1:1026769/1‑150 (MQ=255)
                                         aTACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGa                                                                                                              >  2:1296233/1‑150 (MQ=255)
                                          tACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGAt                                                                                                             <  1:277291/150‑1 (MQ=255)
                                          tACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGAt                                                                                                             <  1:50117/150‑1 (MQ=255)
                                           aCGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGAtt                                                                                                            <  2:85216/150‑1 (MQ=255)
                                            cGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGAttt                                                                                                           >  2:1551888/1‑150 (MQ=255)
                                            cGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGAttt                                                                                                           >  1:196541/1‑150 (MQ=255)
                                             gTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTa                                                                                                          <  2:172071/150‑1 (MQ=255)
                                              tGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTa                                                                                                          >  2:1328085/1‑150 (MQ=255)
                                                aaGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTAtct                                                                                                       <  2:599625/150‑1 (MQ=255)
                                                     tGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACt                                                                                                  <  1:496343/150‑1 (MQ=255)
                                                      ggAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTa                                                                                                 <  1:1078984/150‑1 (MQ=255)
                                                        agaTTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACg                                                                                               <  2:473853/150‑1 (MQ=255)
                                                         cgtgagggTAAAGCCCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGt                                                                                              >  2:657383/8‑150 (MQ=255)
                                                          aTTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTc                                                                                             <  2:922307/150‑1 (MQ=255)
                                                             ttCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCttt                                                                                          <  2:1301011/149‑1 (MQ=255)
                                                             aTCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCttt                                                                                          >  2:73343/1‑150 (MQ=255)
                                                             aTCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGCATTTATCTCAACTACGTCttt                                                                                          >  2:219792/1‑150 (MQ=255)
                                                                 tAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAAtt                                                                                      <  1:1206548/150‑1 (MQ=255)
                                                                 tAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAAtt                                                                                      >  2:21386/1‑150 (MQ=255)
                                                                    aGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACt                                                                                   <  2:778756/150‑1 (MQ=255)
                                                                    aGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATAGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACt                                                                                   <  1:755187/150‑1 (MQ=255)
                                                                        cTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGAtt                                                                               >  1:986902/1‑150 (MQ=255)
                                                                         tttACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTc                                                                              <  1:1395371/150‑1 (MQ=255)
                                                                         tttACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTc                                                                              <  1:849999/150‑1 (MQ=255)
                                                                         tttACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTc                                                                              <  2:16442/150‑1 (MQ=255)
                                                                           tacggtCACGGATAAATTCTGGGGAGAATATAATATTTTCAGTGTGATTTTTCTTATGCATAGCTGCGGTAAAACCAACTGGAACTTTTGATTTGATGACCATATCGGCATAAGGATTTATCTCAATTACGTCTTTAATTAATGATTCTa                                                                            <  2:565132/145‑1 (MQ=255)
                                                                             cccTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTaca                                                                          <  2:455469/150‑1 (MQ=255)
                                                                             cccTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTaca                                                                          >  1:882017/1‑150 (MQ=255)
                                                                               cTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACt                                                                        >  2:930619/1‑150 (MQ=255)
                                                                               cTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACt                                                                        >  1:79021/1‑150 (MQ=255)
                                                                                tCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTg                                                                       <  2:442726/150‑1 (MQ=255)
                                                                                  aCGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGa                                                                     <  1:437059/150‑1 (MQ=255)
                                                                                   cGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGCCCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGAt                                                                    >  2:703668/1‑150 (MQ=255)
                                                                                    ggAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATg                                                                   >  1:692718/1‑150 (MQ=255)
                                                                                     gagaAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGt                                                                  >  1:496080/1‑150 (MQ=255)
                                                                                      agaAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTa                                                                 <  2:1161075/150‑1 (MQ=255)
                                                                                        aaaTTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTAtt                                                               <  2:340535/150‑1 (MQ=255)
                                                                                        aaaTTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTAtt                                                               >  2:468738/1‑150 (MQ=255)
                                                                                         aaTTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTg                                                              >  1:1157192/1‑150 (MQ=255)
                                                                                           ttCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGaa                                                            >  2:950845/1‑150 (MQ=255)
                                                                                            tCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGaaa                                                           <  1:906964/150‑1 (MQ=255)
                                                                                               ggggAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACAATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAaataa                                                        >  1:1488440/1‑150 (MQ=255)
                                                                                                gggAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAaataat                                                       >  1:12730/1‑150 (MQ=255)
                                                                                                   agaATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAATTATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAg                                                    <  2:608439/150‑1 (MQ=255)
                                                                                                       tataatatTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGtttt                                                <  1:145143/150‑1 (MQ=255)
                                                                                                          aatatTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAgg                                             <  1:815150/150‑1 (MQ=255)
                                                                                                           atatTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGa                                            <  1:468738/150‑1 (MQ=255)
                                                                                                            tatTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGAt                                           <  2:1095840/150‑1 (MQ=255)
                                                                                                                     gTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCg                                   <  1:1551888/150‑1 (MQ=255)
                                                                                                                       gCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGt                                 >  1:296669/1‑150 (MQ=255)
                                                                                                                       gCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGt                                 >  1:258294/1‑150 (MQ=255)
                                                                                                                            atTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAg                            <  2:1136816/150‑1 (MQ=255)
                                                                                                                             tttCTTATTCATAGCTGCGGTAAAACCAACGGTAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGt                           <  1:154337/150‑1 (MQ=255)
                                                                                                                                cTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGc                        <  1:1296233/150‑1 (MQ=255)
                                                                                                                                   aTGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCgat                     <  2:24176/150‑1 (MQ=255)
                                                                                                                                    tGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCgatg                    >  1:1383943/1‑150 (MQ=255)
                                                                                                                                       aTCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGatgatg                 >  1:363588/1‑150 (MQ=255)
                                                                                                                                       aTCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGatgatg                 <  2:635618/150‑1 (MQ=255)
                                                                                                                                        tCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGAtgatga                >  1:164530/1‑150 (MQ=255)
                                                                                                                                        tCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGAtgatga                <  2:1157192/150‑1 (MQ=255)
                                                                                                                                             gCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATaa           >  1:292438/1‑150 (MQ=255)
                                                                                                                                              tcatcaaaTCAACGGTTCCCGTTGGTTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAAt          >  2:1259734/6‑150 (MQ=255)
                                                                                                                                               catcaaaTCAACGGTTCCCGTTGGTTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATc         >  2:566314/5‑150 (MQ=255)
                                                                                                                                                gTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCCCATAAGGATTTATCTCAACTACGTCTTAAATCACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCa        <  2:1070951/150‑1 (MQ=255)
                                                                                                                                                atcaaaTCAACGGTTCCCGTTGGTTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCa        >  2:898173/4‑150 (MQ=255)
                                                                                                                                                 tAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAg                            <  2:1042199/129‑1 (MQ=255)
                                                                                                                                                 tAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAg                            >  1:1042199/1‑129 (MQ=255)
                                                                                                                                                  aaaaCCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCAGc      <  2:384946/150‑1 (MQ=255)
                                                                                                                                                   aaaTCAACGGTTCCCGTTGGTTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCAGCa     >  1:576933/1‑150 (MQ=255)
                                                                                                                                                   aaaTCAACGGTTCCCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCAGCa     >  2:1479425/1‑150 (MQ=255)
                                                                                                                                                      ccAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCAGCATcc  >  2:826000/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAGCGAAACGTTCTGCGCGTTCTGAACGCTCACCGATGACAATACGTGAAGGATGGAGATTATCGTAAAGGGCTTTACCCTCACGGAGAAATTCCGGGGAGAATATAATATTTTCAGTGCGATATTTCTTATGCATCGCTGCGGTAAAACCAACGGGAACCGTTGATTTGATGACCATAACCGCATAAGGATTTATCTCAACTACGTCTTTAATTACTGATTCTACACTGGATGTATTGAAATAATTAGTTTTAGGATCATAGTCGGTTGGAGTGGCGATGATGACATAATCAGCATCC  >  NC_000913/2099100‑2099398

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: