Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 691,390 G→C 6.5% A132G (GCC→GGC)  ybeX ← putative ion transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,3900GC 6.5% 471.7 / 19.6 138A132G (GCC→GGC) ybeXputative ion transport
Reads supporting (aligned to +/- strand):  ref base G (57/72);  new base C (4/5);  total (61/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.14e-01

Only 100 of 140 total aligned reads displayed.

AATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTT  >  NC_000913/691242‑691538
                                                                                                                                                    |                                                                                                                                                    
aaTGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAgg                                                                                                                                                      >  2:514099/1‑149 (MQ=255)
    gTCACTAGACCGGAAACCCCCCCGAATTCGTCAATAACGATCGCCATGTGGTAACGTTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGGGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTc                                                                                                                                                  >  4:901433/1‑149 (MQ=255)
     tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                 >  2:784683/1‑149 (MQ=255)
     tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                 <  1:797759/149‑1 (MQ=255)
     tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                 <  2:212040/149‑1 (MQ=255)
     tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                 >  3:1016533/1‑149 (MQ=255)
     tCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACATAACACTTTGTCCATGCTGAAGGCTTCa                                                                                                                                                 <  2:612736/149‑1 (MQ=255)
      cACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAg                                                                                                                                                <  1:845374/149‑1 (MQ=255)
         cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCAt                                                                                                                                             <  4:795884/149‑1 (MQ=255)
         cAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAt                                                                                                                                             <  4:805242/149‑1 (MQ=255)
              cGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTg                                                                                                                                        >  4:113300/1‑149 (MQ=255)
               ggAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTgc                                                                                                                                       <  3:605647/149‑1 (MQ=255)
                   aCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCAt                                                                                                                                   <  4:1010887/149‑1 (MQ=255)
                   aCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCAt                                                                                                                                   >  2:343157/1‑149 (MQ=255)
                   aCCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCAt                                                                                                                                   <  2:265524/149‑1 (MQ=255)
                     cccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATc                                                                                                                                            >  2:622335/1‑138 (MQ=255)
                     cccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATc                                                                                                                                            <  1:622335/138‑1 (MQ=255)
                     cccACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATaa                                                                                                                                 >  3:334631/1‑149 (MQ=255)
                       cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAAc                                                                                                                               >  2:584275/1‑149 (MQ=255)
                       cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                            <  3:961757/136‑2 (MQ=255)
                       cACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                            >  4:961757/1‑135 (MQ=255)
                               tCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCaa                                                                                                                       <  1:689550/149‑1 (MQ=255)
                               tCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCaa                                                                                                                       <  1:784683/149‑1 (MQ=255)
                                  tCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTc                                                                                                                    >  1:13651/1‑149 (MQ=255)
                                   cAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                            <  2:167665/124‑2 (MQ=255)
                                   cAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAtg                                                                                                                                            >  1:167665/1‑123 (MQ=255)
                                      tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTc                                                                                                                <  1:125097/149‑1 (MQ=255)
                                      tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTc                                                                                                                <  4:39144/149‑1 (MQ=255)
                                      tAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTc                                                                                                                <  4:241113/149‑1 (MQ=255)
                                          gATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                                                                                            <  1:303373/149‑1 (MQ=255)
                                          gATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                                                                                            <  3:524406/149‑1 (MQ=255)
                                          gATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCa                                                                                                            >  4:113923/1‑149 (MQ=255)
                                           aTCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCAt                                                                                                           <  4:188037/149‑1 (MQ=255)
                                            tCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATc                                                                                                          <  1:363214/149‑1 (MQ=255)
                                            tCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATc                                                                                                          >  4:812020/1‑149 (MQ=255)
                                             cGCCAGGTGGTATCGCTGAGAGCGAAACTCTTTGAGCATCCGGTCTACGCGCTTACTTTCGGGAACGACAACCGCCTGGCGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATATACGGCAGCACGTCTTTCGCCATCa                                                                                                         <  2:166703/149‑1 (MQ=255)
                                               ccATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAt                                                                                                                                             >  3:958575/1‑111 (MQ=255)
                                               ccATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGCCTTCAGCAt                                                                                                                                             <  4:958575/111‑1 (MQ=255)
                                                cATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGTAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGaa                                                                                                      <  2:250458/149‑1 (MQ=255)
                                                     ggTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCt                                                                                                 <  1:645999/149‑1 (MQ=255)
                                                         aCGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCaa                                                                                             <  1:244892/149‑1 (MQ=255)
                                                           gCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtg                                                                                           <  3:70658/149‑1 (MQ=255)
                                                             tGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAAtgtg                                                                                         <  1:156495/149‑1 (MQ=255)
                                                              gagaGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGa                                                                                        <  1:241044/149‑1 (MQ=255)
                                                               agagCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGAt                                                                                       <  1:304514/149‑1 (MQ=255)
                                                                   cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                   <  2:985977/149‑1 (MQ=255)
                                                                   cGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCttt                                                                                   <  2:411411/149‑1 (MQ=255)
                                                                    gAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTg                                                                                  <  4:296257/149‑1 (MQ=255)
                                                                       aCTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCt                                                                               <  2:1040379/149‑1 (MQ=255)
                                                                         tctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTc                                                                             <  1:33368/149‑1 (MQ=255)
                                                                          ctTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTCACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCg                                                                            <  3:863315/149‑1 (MQ=255)
                                                                           tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                           >  1:506379/1‑149 (MQ=255)
                                                                           tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                           >  1:586545/1‑149 (MQ=255)
                                                                           tttCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                           <  1:703501/149‑1 (MQ=255)
                                                                              cAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTaa                                                                        >  1:554225/1‑149 (MQ=255)
                                                                                 cATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAACGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCa                                                                     >  4:727797/1‑149 (MQ=255)
                                                                                   tCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCAcc                                                                   >  2:746666/1‑149 (MQ=255)
                                                                                     cGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCgg                                                                 >  2:421002/1‑149 (MQ=255)
                                                                                       gTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                           <  2:842302/137‑1 (MQ=255)
                                                                                       gTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGc                                                                           >  1:842302/1‑137 (MQ=255)
                                                                                               gcTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAg                                                       >  1:18606/1‑149 (MQ=255)
                                                                                                cTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGa                                                        <  4:475111/147‑1 (MQ=255)
                                                                                                cTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGa                                                        >  3:475111/1‑147 (MQ=255)
                                                                                                 ttACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTg                                                     >  4:851924/1‑149 (MQ=255)
                                                                                                    cTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                                  <  4:334631/149‑1 (MQ=255)
                                                                                                     tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAGTCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCg                                                 >  4:869693/1‑149 (MQ=255)
                                                                                                     tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa                                                             >  2:896729/1‑137 (MQ=255)
                                                                                                     tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGaaa                                                             <  1:896729/137‑1 (MQ=255)
                                                                                                     tttCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCg                                                 >  1:171031/1‑149 (MQ=255)
                                                                                                         aGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACt                                             <  3:812020/149‑1 (MQ=255)
                                                                                                          ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                            >  4:548097/1‑133 (MQ=255)
                                                                                                          ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAAc                                                            <  3:548097/133‑1 (MQ=255)
                                                                                                          ggAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTc                                            >  3:228912/1‑149 (MQ=255)
                                                                                                           gAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCg                                           >  4:845529/1‑149 (MQ=255)
                                                                                                              cGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatg                                        >  1:370882/1‑149 (MQ=255)
                                                                                                               gACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatga                                       <  4:403272/149‑1 (MQ=255)
                                                                                                                aCAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatgat                                      <  3:428010/149‑1 (MQ=255)
                                                                                                                 cAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGatgatg                                     <  3:40886/149‑1 (MQ=255)
                                                                                                                 cAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGCCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGatgatg                                     <  3:727797/149‑1 (MQ=255)
                                                                                                                  aaCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                                  >  1:608232/1‑135 (MQ=255)
                                                                                                                  aaCCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGc                                                  <  2:608232/135‑1 (MQ=255)
                                                                                                                             gTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACAt                         >  1:350463/1‑149 (MQ=255)
                                                                                                                              tAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGTGTGGACTCGATGATGACATCAAGACAtt                        >  2:149064/1‑149 (MQ=255)
                                                                                                                               aaCACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTc                       <  1:894808/149‑1 (MQ=255)
                                                                                                                               aaCACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTc                       <  2:747646/149‑1 (MQ=255)
                                                                                                                                  acTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTc                    <  2:831819/149‑1 (MQ=255)
                                                                                                                                    tttGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCgatgat                                      <  4:320248/129‑1 (MQ=255)
                                                                                                                                          gaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACa                                  >  3:870377/2‑127 (MQ=255)
                                                                                                                                          gaTGCTGAAGCCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACa                                  <  4:870377/126‑1 (MQ=255)
                                                                                                                                          cATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt            >  1:973502/1‑149 (MQ=255)
                                                                                                                                          cATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt            <  3:113300/149‑1 (MQ=255)
                                                                                                                                          cATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCt            >  1:1028002/1‑149 (MQ=255)
                                                                                                                                            tGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTgg          <  3:113923/149‑1 (MQ=255)
                                                                                                                                                gAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGc      >  4:709394/1‑149 (MQ=255)
                                                                                                                                                gAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGc      >  1:787129/1‑149 (MQ=255)
                                                                                                                                                gAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGc      <  3:33666/149‑1 (MQ=255)
                                                                                                                                                 aaGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTc                       <  3:17055/131‑1 (MQ=255)
                                                                                                                                                 aaGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCg     >  1:121884/1‑149 (MQ=255)
                                                                                                                                                  aGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGt    >  4:76533/1‑149 (MQ=255)
                                                                                                                                                    gCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGttt  <  1:421002/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATGGTCACCAGACCGGAAACCCCACCGAATTCGTCAATAACGATCGCCATGTGGTAACGCTGAGAGCGAAACTCTTTCAGCATCCGGTCTACGCGCTTACTTTCAGGAACGACAACCGCCTGACGTAACACTTTGTCCATGCTGAAGGCTTCAGCATCGCTGCGCATAAACGGCAGCAAGTCTTTCGCCATCAGAATCCCTTCAATGTGATCTTTGTCTTCGCTAATCACCGGGAAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTT  >  NC_000913/691242‑691538

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