Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,100,879 A→G 5.5% G130G (GGT→GGC gnd ← 6‑phosphogluconate dehydrogenase, decarboxylating

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,100,8790AG 5.5% 275.1 / 5.5 91G130G (GGT→GGCgnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base G (3/2);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

TGCGCCATCGGCACCAATATAGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCA  >  NC_000913/2100735‑2101026
                                                                                                                                                |                                                                                                                                                   
tGCGCCATCGGCACCAATATAGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGc                                                                                                                                                 >  5:294801/1‑149 (MQ=255)
  cgcCATCGGCACCAATATAGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCa                                                                                                                                               <  1:161266/149‑1 (MQ=255)
       tCGGCACCAATATAGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATATTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCAAGaaa                                                                                                                                           >  2:110209/1‑148 (MQ=255)
               aaTATAGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTAc                                                                                                                                  >  1:105950/1‑149 (MQ=255)
                    aGGTACCGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATTATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGAt                                                                                                                              >  4:124159/1‑148 (MQ=255)
                    aGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGAt                                                                                                                              >  6:149579/1‑148 (MQ=255)
                      gTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATg                                                                                                                             >  5:213520/1‑147 (MQ=255)
                       tAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATg                                                                                                                             >  1:102498/1‑146 (MQ=255)
                       tAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGaa                                                                                                                           <  2:102498/148‑1 (MQ=255)
                             aTGGTTCACCGTCTTCATCTACTGCGTCGATTTTGGTAAGGATCGGTGCTACCTATTCCTAGGCTTCTTTCTGGCGACAAGGCATAATAGAAGGACCTTTCAGCGCCCGCTCTTCACCGCCAGAAACACCGGTACCGATGAAGGTAAAg                                                                                                                    <  1:72493/149‑1 (MQ=255)
                              tGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTaa                                                                                                                      <  2:119167/146‑1 (MQ=255)
                              tGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTaa                                                                                                                      >  1:119167/1‑146 (MQ=255)
                               ggTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGcc                                                                                                                  >  4:136388/1‑149 (MQ=255)
                                   cACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCctc                                                                                                               <  8:260521/148‑1 (MQ=255)
                                        tCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTTGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTctgctg                                                                                                          >  8:167911/1‑148 (MQ=255)
                                        tCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTctgctg                                                                                                          >  2:92843/1‑148 (MQ=255)
                                         cTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGa                                                                                                         >  7:196056/1‑148 (MQ=255)
                                          ttCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGa                                                                                                         >  8:337909/1‑147 (MQ=255)
                                            gaGCTACGGCGGCGTTTTTGGTCAGGATCGGGGCTACCAAATCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGaaa                                                                                                       <  6:294801/146‑1 (MQ=255)
                                               cTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCa                                                                                                  >  4:103571/1‑149 (MQ=255)
                                                     cggcgATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTa                                                                                            <  1:52773/149‑1 (MQ=255)
                                                     aggcgATTGTGGGCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGACAACAGGCATAATAGAAGGACCTTTCAGCGCCCCATCTTCCCCGACAGAAAAAACGGTACCGATGATTTTAAAGCCATCTTCTTCAATCTCAAGAtt                                                                                             >  4:275083/2‑148 (MQ=255)
                                                      ggcgATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCACTCTGCTGAAAGATCACGATTac                                                                                           >  2:275924/1‑149 (MQ=255)
                                                           ttttGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTAcgacga                                                                                      <  7:167911/149‑1 (MQ=255)
                                                               ggTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATa                                                                                   <  7:298756/148‑1 (MQ=255)
                                                               ggTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATa                                                                                   <  5:149579/148‑1 (MQ=255)
                                                               ggTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATa                                                                                   <  7:337909/148‑1 (MQ=255)
                                                                gTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAgt                                                                                 <  8:196056/149‑1 (MQ=255)
                                                                gTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAgt                                                                                 <  2:105950/149‑1 (MQ=255)
                                                                gTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAg                                                                                  <  6:307032/148‑1 (MQ=255)
                                                                   aGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCg                                                                                                                                      >  8:138175/1‑93 (MQ=255)
                                                                      aTCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTg                                                                           <  3:133019/149‑1 (MQ=255)
                                                                         ggTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGgaa                                                                        <  1:92843/149‑1 (MQ=255)
                                                                             cTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGgaagaa                                                                     <  6:213520/148‑1 (MQ=255)
                                                                               aCCAATTCTTAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgt                                                                  >  4:109582/1‑149 (MQ=255)
                                                                               aCCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgt                                                                  >  7:366786/1‑149 (MQ=255)
                                                                               aCCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgt                                                                  >  1:373735/1‑149 (MQ=255)
                                                                                ccAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgtg                                                                 >  1:243790/1‑149 (MQ=255)
                                                                                 cAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgtgt                                                                >  5:29142/1‑149 (MQ=255)
                                                                                 cAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGgtgt                                                                <  8:218358/149‑1 (MQ=255)
                                                                                   aTTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCGCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTa                                                                                            >  4:271596/1‑119 (MQ=255)
                                                                                   aTTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCGCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTa                                                                                            >  4:276149/1‑119 (MQ=255)
                                                                                   aTTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCGCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTa                                                                                            <  3:276149/119‑1 (MQ=255)
                                                                                   aTTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCGCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTa                                                                                            <  3:271596/119‑1 (MQ=255)
                                                                                       aTAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTaccac                                                          >  5:210661/1‑149 (MQ=255)
                                                                                         aGGCTTCTTTCTGGCCACCGGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCAc                                                                                                 <  5:139113/108‑1 (MQ=255)
                                                                                         aGGCTTCTTTCTGGCCACCGGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCAc                                                                                                 >  6:139113/1‑108 (MQ=255)
                                                                                           gCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCAt                                                       >  1:98486/1‑148 (MQ=255)
                                                                                            cttcttTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGgaa                                                                        <  8:231259/130‑1 (MQ=255)
                                                                                            cttcttTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGgaa                                                                        >  7:231259/1‑130 (MQ=255)
                                                                                             ttctttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCaa                                                    >  2:192407/1‑149 (MQ=255)
                                                                                              tctttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCaa                                                    >  2:323734/1‑148 (MQ=255)
                                                                                               ctttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTAGGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAat                                                   <  1:110209/148‑1 (MQ=255)
                                                                                               ctttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAgt                                                                                 >  6:49022/1‑118 (MQ=255)
                                                                                               ctttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAgt                                                                                 <  5:49022/118‑1 (MQ=255)
                                                                                               ctttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAat                                                   <  6:210661/148‑1 (MQ=255)
                                                                                                 ttctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgat                                                >  2:283475/1‑149 (MQ=255)
                                                                                                  tctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgat                                                <  3:103571/148‑1 (MQ=255)
                                                                                                   ctGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatg                                               <  2:243790/148‑1 (MQ=255)
                                                                                                    tGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACCCAGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTTTCCTGGAAGAAGGTGTTACCACAATCAatgatgat                                             >  2:122629/1‑149 (MQ=255)
                                                                                                    tGGCCACCAGGCAAAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAACCACCGGTACCGATGAAGTTAAAGCCCGCTGATGAAACCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACAATCAatcatgag                                             >  1:266639/1‑148 (MQ=255)
                                                                                                      gCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTCCGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatgat                                             <  1:275924/147‑1 (MQ=255)
                                                                                                      gCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAatgatgat                                             <  6:352221/147‑1 (MQ=255)
                                                                                                        caccaGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGtct                                         >  5:104127/1‑149 (MQ=255)
                                                                                                         accaGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGtctc                                        >  8:70589/1‑149 (MQ=255)
                                                                                                                ataataGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATc                                 <  6:29142/149‑1 (MQ=255)
                                                                                                                ataataGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATc                                 >  8:93550/1‑149 (MQ=255)
                                                                                                                 taataGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATc                                 >  1:91846/1‑148 (MQ=255)
                                                                                                                       aaGGACCTTTCAGCGCCCCCTCTTCGCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGatat                          >  8:99138/1‑149 (MQ=255)
                                                                                                                       aaGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGatat                          >  6:296414/1‑149 (MQ=255)
                                                                                                                       aaGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGatat                          <  8:366786/149‑1 (MQ=255)
                                                                                                                        aGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATg                         <  3:136388/149‑1 (MQ=255)
                                                                                                                        aGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATg                         <  7:665/149‑1 (MQ=255)
                                                                                                                            ccTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTTTCCTTTATCGAGATATGGttt                     >  6:305338/1‑149 (MQ=255)
                                                                                                                            ccTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGttt                     <  2:373735/149‑1 (MQ=255)
                                                                                                                             cTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGttt                     >  6:197825/1‑148 (MQ=255)
                                                                                                                             cTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGttt                     >  8:81420/1‑148 (MQ=255)
                                                                                                                              tttCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTcc                                                                             >  2:373437/1‑91 (MQ=255)
                                                                                                                                tccgcgcCCCCTCTTCACCGCCAGAAACACAGGGACCGAGGAAGTTAACGCCCTCTTCTGAAATCTCAAGATTACTCCGAATAGTGTCATGGAAGAATGGGTTACCACCATCAATGATGATGTCTCCTTTATCTAGATATGGTTTGa                   <  3:275083/144‑1 (MQ=255)
                                                                                                                                  aGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGa               <  7:70589/149‑1 (MQ=255)
                                                                                                                                     gcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGgaa                                                                        >  8:304336/1‑89 (MQ=255)
                                                                                                                                     gcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGgaa                                                                        <  7:304336/89‑1 (MQ=255)
                                                                                                                                     gcCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATc            <  2:98486/149‑1 (MQ=255)
                                                                                                                                       ccccTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCaa          <  1:283475/149‑1 (MQ=255)
                                                                                                                                        cccTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCAGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGAGATGGTTTGAGGGGATCaa          >  1:144680/1‑148 (MQ=255)
                                                                                                                                         ccTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATa        <  6:104127/149‑1 (MQ=255)
                                                                                                                                           tctTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATag       >  4:255222/1‑148 (MQ=255)
                                                                                                                                             atCACCGCAATAAACACCTTTAGCTATGAAGTTAAAGACCTCTGCTGAAAGCTCACGATTTCGACGAATAGTGTCGTGGAAGAATGTGTTAGCACCTTCAATGATTATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATagca     <  1:122629/147‑1 (MQ=255)
                                                                                                                                               cACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATagca     <  3:124159/146‑1 (MQ=255)
                                                                                                                                               cACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATagca     <  7:81420/146‑1 (MQ=255)
                                                                                                                                               cACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAgcagca  >  2:30915/1‑149 (MQ=255)
                                                                                                                                                aCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATag       >  8:67926/1‑143 (MQ=255)
                                                                                                                                                aCCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATag       <  7:67926/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGCGCCATCGGCACCAATATAGGTAACGCATGGTTCACCGTCTTCAGCTACGGCGGCGATTTTGGTCAGGATCGGTGCTACCAATTCATAGGCTTCTTTCTGGCCACCAGGCATAATAGAAGGACCTTTCAGCGCCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCA  >  NC_000913/2100735‑2101026

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