Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,495,814 C→T 5.5% W132* (TGG→TGA ahr ← broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,495,8140CT 5.5% 235.5 / 2.7 74W132* (TGG→TGAahrbroad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base C (33/36);  new base T (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGA  >  NC_000913/4495667‑4495962
                                                                                                                                                   |                                                                                                                                                    
ccccAGCCCGCCAATACCAATTACCCCACCGCGGCTGGTAGCAGTGATATGGTGCATCAACAGGGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGcc                                                                                                                                                     <  1:183227/149‑1 (MQ=255)
ccccAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGc                                                                                                                                                      >  7:97705/1‑148 (MQ=255)
ccccAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGc                                                                                                                                                      <  3:262850/148‑1 (MQ=255)
 cccAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGc                                                                                                                                                      >  4:108686/1‑147 (MQ=255)
  ccAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAg                                                                                                                                                   >  1:276983/1‑149 (MQ=255)
  ccAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAg                                                                                                                                                   <  1:170590/149‑1 (MQ=255)
         gCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGCCCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCcgcac                                                                                                                                            >  3:258497/1‑149 (MQ=255)
             aTACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCTGACTCGATATCAATATTTTCTGGCAGTGGAATCACACATTGACAGTCCGCACGCaa                                                                                                                                        >  8:201708/1‑149 (MQ=255)
             aTACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCaa                                                                                                                                        >  8:166476/1‑149 (MQ=255)
             aTACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCaa                                                                                                                                        <  7:166476/149‑1 (MQ=255)
                    ttACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCg                                                                                                                                 <  3:230661/149‑1 (MQ=255)
                       ccccAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCTACTTCTCGGc                                                                                                                               >  3:265062/1‑148 (MQ=255)
                        cccAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCa                                                                                                                              >  6:117259/1‑148 (MQ=255)
                         ccAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCaaa                                                                                                                            <  1:250029/149‑1 (MQ=255)
                                ggCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGc                                                                                                                          >  6:265283/1‑144 (MQ=255)
                                ggCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGc                                                                                                                          <  5:265283/144‑1 (MQ=255)
                                 gCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCAAGCAACTTCTCGGCAAAGCCACcgc                                                                                                                    >  4:240794/1‑149 (MQ=255)
                                     gTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGAtt                                                                                                                <  7:126179/149‑1 (MQ=255)
                                        gCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATGCCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCa                                                                                                              >  5:299094/1‑148 (MQ=255)
                                        gCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCACCTTCTCGGCAAAGCCCCCGCGATTCa                                                                                                              >  2:113884/1‑148 (MQ=255)
                                          aGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATaa                                                                                                           >  4:385961/1‑149 (MQ=255)
                                            tGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAt                                                                                                          <  6:235581/148‑1 (MQ=255)
                                               tatGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      >  1:192131/1‑149 (MQ=255)
                                                 tGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCgg                                                                                                    <  3:108686/149‑1 (MQ=255)
                                                 tGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCgg                                                                                                    <  4:258497/149‑1 (MQ=255)
                                                       aTCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGgcaccgc                                                                                              <  8:382122/149‑1 (MQ=255)
                                                         cAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGgcaccgcac                                                                                            >  8:204223/1‑149 (MQ=255)
                                                           aCAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      <  7:56001/137‑1 (MQ=255)
                                                           aCAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      >  8:56001/1‑137 (MQ=255)
                                                              gTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGc                                                                                       <  5:117259/149‑1 (MQ=255)
                                                              gTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGc                                                                                       <  7:204223/149‑1 (MQ=255)
                                                                  tttAAAGACCGTGATACCTCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATaa                                                                                                           <  4:36603/125‑1 (MQ=255)
                                                                  tttAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATaa                                                                                                           >  3:36603/1‑125 (MQ=255)
                                                                     aaaGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTtctc                                                                                                                                  <  7:274944/99‑1 (MQ=255)
                                                                     aaaGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTtctc                                                                                                                                  >  8:274944/1‑99 (MQ=255)
                                                                         aCCGTGATACCGCGGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGTCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      >  3:238539/1‑123 (MQ=255)
                                                                         aCCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGTCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      <  4:233335/123‑1 (MQ=255)
                                                                         aCCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGTCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      >  3:233335/1‑123 (MQ=255)
                                                                         aCCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGTCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAAtcgtc                                                                                                      <  4:238539/123‑1 (MQ=255)
                                                                         aCCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGAt                                                                            <  4:3766/149‑1 (MQ=255)
                                                                          ccGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATc                                                                           <  3:385961/149‑1 (MQ=255)
                                                                          ccGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATc                                                                           <  8:422120/149‑1 (MQ=255)
                                                                           cGTGATGCCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCt                                                                          <  6:299094/149‑1 (MQ=255)
                                                                           cGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCt                                                                          <  1:37733/149‑1 (MQ=255)
                                                                                      cgcACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAg                                                                                 >  8:321261/1‑131 (MQ=255)
                                                                                      cgcACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAg                                                                                 <  7:321261/131‑1 (MQ=255)
                                                                                       gcACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTaa                                                              >  8:31218/1‑149 (MQ=255)
                                                                                       gcACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTaa                                                              >  2:395643/1‑149 (MQ=255)
                                                                                             caGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAgg                                                        <  2:347502/149‑1 (MQ=255)
                                                                                                 ggCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTc                                                    <  2:216958/149‑1 (MQ=255)
                                                                                                     cggcggACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGTCAAAGCCACCGCGATTCATAATCGTCGGCACCGACCCTTGCTCGCAGTTGATCTGATTACCGATAATACAGGCGTCGCAg                                                >  5:244911/1‑149 (MQ=255)
                                                                                                     cggcggACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAg                                                <  1:395643/149‑1 (MQ=255)
                                                                                                        cggACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGa                                             <  7:31218/149‑1 (MQ=255)
                                                                                                          gACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGAc                                            >  5:144090/1‑148 (MQ=255)
                                                                                                            cTCGATATAAATATTTGCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGGCGCAGTGACcac                                         <  4:25727/149‑1 (MQ=255)
                                                                                                             tCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAg                                                         >  8:426580/1‑132 (MQ=255)
                                                                                                             tCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAg                                                         <  7:426580/132‑1 (MQ=255)
                                                                                                               gATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTAc                                                                    >  7:262128/1‑119 (MQ=255)
                                                                                                               gATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTAc                                                                    <  8:262128/119‑1 (MQ=255)
                                                                                                                 tatCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTa                                    >  7:382770/1‑149 (MQ=255)
                                                                                                                   tCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTAc                                   <  6:144090/148‑1 (MQ=255)
                                                                                                                     aaTATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTAcgcg                                >  4:52887/1‑149 (MQ=255)
                                                                                                                     aaTATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCAAGCAACTTCTCTGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTAcgc                                 <  3:240794/148‑1 (MQ=255)
                                                                                                                             cTGGGAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATACTCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCa                         <  4:238453/148‑1 (MQ=255)
                                                                                                                                gCAGTGGAATCACCCATTGCCAGTCCGCACGCTACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGc                       <  4:265062/147‑1 (MQ=255)
                                                                                                                                gCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCCCCTTGCTCGCAGTTGATCTGCTTACCGCTAATACATGCGTAGCAGTGACCACAGCTACGCGACGTACAGcc                      >  6:232897/1‑148 (MQ=255)
                                                                                                                                gCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAg                        <  2:276983/146‑1 (MQ=255)
                                                                                                                                   gTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCCCCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCCATc                  >  7:241524/1‑149 (MQ=255)
                                                                                                                                     ggAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATccc                <  8:97705/149‑1 (MQ=255)
                                                                                                                                           aCCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGACACg          >  4:104844/1‑149 (MQ=255)
                                                                                                                                           aCCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGACACg          >  7:42861/1‑149 (MQ=255)
                                                                                                                                             ccATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACa                                                          <  8:190327/99‑1 (MQ=255)
                                                                                                                                             ccATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACa                                                          >  7:190327/1‑99 (MQ=255)
                                                                                                                                              cATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGACACGCt        >  6:317040/1‑148 (MQ=255)
                                                                                                                                                ttGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGac             >  8:30112/1‑141 (MQ=255)
                                                                                                                                                ttGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGac             <  7:30112/141‑1 (MQ=255)
                                                                                                                                                ttGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGACACGCTgac     >  1:39589/1‑149 (MQ=255)
                                                                                                                                                 tGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCCGTTGATCTGATTACCGCTAATACAGGAGTCGCGGTGACCACAGCTACGCGCCGTCGAGCCAATCCCGACACGCTgccc    >  2:358202/1‑149 (MQ=255)
                                                                                                                                                   ccAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCCCAGCTACGCGCCGTCCAGCCAATCCCGACACGCTgaccga  >  4:162878/1‑149 (MQ=255)
                                                                                                                                                   ccAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATACATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATCCAGGCGTCGCAGTGACCACAGCTACGCGCCGGCCAGCCAATCCCGACACGCTgaccga  >  4:290690/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCCAGCCCGCCAATACCAATTACCCCAACGCGGCTGGTAGCAGTGATATGGTGCATCAACAGTGGTTTAAAGACCGTGATACCGCCGCACAACAGCGGCCCGGCGGACTCGATATCAATATTTTCTGGCAGTGGAATCACCCATTGCCAGTCCGCACGCAACTTCTCGGCAAAGCCACCGCGATTCATAATCGTCGGCACCGCACCTTGCTCGCAGTTGATCTGATTACCGCTAATACAGGCGTCGCAGTGACCACAGCTACGCGCCGTCCAGCCAATCCCGACACGCTGACCGA  >  NC_000913/4495667‑4495962

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