Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,137,909 A→C 5.5% intergenic (‑252/+23) yijF ← / ← gldA DUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,9090AC 5.5% 374.5 / 10.5 126intergenic (‑252/+23)yijF/gldADUF1287 family protein/glycerol dehydrogenase, NAD+ dependent; 1,2‑propanediol:NAD+ oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (53/66);  new base C (4/3);  total (57/69)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

Only 100 of 127 total aligned reads displayed.

TGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTT  >  NC_000913/4137761‑4138056
                                                                                                                                                    |                                                                                                                                                   
tGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGa                                                                                                                                                     <  7:38465/149‑1 (MQ=255)
     gTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGc                                                                                                                                                >  8:40563/1‑149 (MQ=255)
      ttacacaTTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGACTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGcc                                                                                                                                               >  6:341488/7‑149 (MQ=255)
         cAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCgg                                                                                                                                             <  7:377577/148‑1 (MQ=255)
             ttACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGtt                                                                                                                                        <  3:107972/149‑1 (MQ=255)
             ttACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGtt                                                                                                                                        <  2:316574/149‑1 (MQ=255)
              tACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGttt                                                                                                                                       >  4:342082/1‑149 (MQ=255)
                cTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTgg                                                                                                                                     <  8:352407/149‑1 (MQ=255)
                   gCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTgaa                                                                                                                                    <  5:178026/147‑3 (MQ=255)
                   gCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGt                                                                                                                                  <  2:72439/149‑1 (MQ=255)
                    cacaGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTa                                                                                                                                 <  8:160042/149‑1 (MQ=255)
                     acaGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAg                                                                                                                                <  8:133648/149‑1 (MQ=255)
                      caGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAgg                                                                                                                               <  7:40563/149‑1 (MQ=255)
                       agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTg                                                                                                                                      <  2:153770/141‑1 (MQ=255)
                       agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTg                                                                                                                                      >  1:153770/1‑141 (MQ=255)
                       agagTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGt                                                                                                                              >  7:301173/1‑149 (MQ=255)
                         agTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTa                                                                                                                            >  7:258876/1‑149 (MQ=255)
                         agTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTa                                                                                                                            <  5:201643/149‑1 (MQ=255)
                           tGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTg                                                                                                                                      <  5:343278/137‑1 (MQ=255)
                            gCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTc                                                                                                                         <  2:145789/149‑1 (MQ=255)
                             ccGGATGCGTAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTcc                                                                                                                        <  4:279210/149‑1 (MQ=255)
                             ccGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTcc                                                                                                                        >  6:236008/1‑149 (MQ=255)
                             ccGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATAAGGAAATTTTTCGTTCAAACTCCCTGACAAGCCGGGAGTTTGTAGTAGGTTATTc                                                                                                                         >  6:230176/1‑148 (MQ=21)
                             ccGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCCCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGACTAGGTTATTc                                                                                                                         >  3:158514/1‑148 (MQ=255)
                                 aTGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCAct                                                                                                                    <  7:147031/149‑1 (MQ=255)
                                  tGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCActc                                                                                                                   <  7:281858/149‑1 (MQ=255)
                                      gAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATTTTGCAGAGATTGCGTAGTCCTGATAAGCGTAGGGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTg                                                                                                                >  6:342216/1‑148 (MQ=255)
                                      gAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATACGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGc                                                                                                               <  5:236008/149‑1 (MQ=255)
                                       aGTTAACGCCTTATACTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGc                                                                                                               <  5:230176/148‑1 (MQ=255)
                                       aGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCa                                                                                                              <  3:41515/149‑1 (MQ=255)
                                       aGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCa                                                                                                              <  8:166233/149‑1 (MQ=255)
                                         tGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAg                                                                                                             >  3:207814/1‑148 (MQ=255)
                                             cGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAAc                                                                                                        >  1:51671/1‑149 (MQ=255)
                                              gCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACg                                                                                                       <  5:38011/149‑1 (MQ=255)
                                                cTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCt                                                                                                     >  6:125391/1‑149 (MQ=255)
                                                cTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCt                                                                                                     >  5:182891/1‑149 (MQ=255)
                                                       ggCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTCCAAACTCCCGGCTTGTCCGGGAGTTTgaac                                                                                                                                   <  5:211911/112‑4 (MQ=17)
                                                       ggCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTCCAAACTCCCGGCTTGTCCGGGAGTTTgaac                                                                                                                                   >  6:211911/1‑109 (MQ=17)
                                                       ggCCTACAAAAGCAAGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTa                                                                                              <  1:43931/149‑1 (MQ=255)
                                                           tACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTAt                                                                                                                           <  8:401985/116‑1 (MQ=255)
                                                           tACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTAt                                                                                                                           >  7:401985/1‑116 (MQ=255)
                                                              aaaaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGc                                                                                                               >  2:337110/1‑125 (MQ=255)
                                                              aaaaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGc                                                                                                               <  1:337110/125‑1 (MQ=255)
                                                              aaaaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCg                                                                                       >  4:237485/1‑149 (MQ=255)
                                                                aaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGc                                                                                     <  4:63099/149‑1 (MQ=255)
                                                                aaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGc                                                                                     <  8:373918/149‑1 (MQ=255)
                                                                aaGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGc                                                                                     <  8:114772/149‑1 (MQ=255)
                                                                   cACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTAc                                                                                  <  5:341488/149‑1 (MQ=255)
                                                                      gCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcag                                                                               <  8:258876/149‑1 (MQ=255)
                                                                      gCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcag                                                                               >  5:85282/1‑149 (MQ=255)
                                                                          aTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                             >  2:217852/1‑147 (MQ=255)
                                                                          aTTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGAGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACcagca                                                                             >  3:100000/1‑147 (MQ=255)
                                                                            tCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgc                                                                         >  2:35507/1‑149 (MQ=255)
                                                                              aaTATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcg                                                                        <  4:42227/148‑1 (MQ=255)
                                                                               atatatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAgcggc                                                                      <  4:158514/149‑1 (MQ=255)
                                                                                  tatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAg                                                                                                                                   <  6:1767/85‑1 (MQ=255)
                                                                                  tatTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAg                                                                                                                                   >  5:1767/1‑85 (MQ=255)
                                                                                    ttGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTaaa                                                                 <  7:305302/149‑1 (MQ=255)
                                                                                      gCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAAc                                                                <  6:182891/148‑1 (MQ=255)
                                                                                       cAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCg                                                                                       >  6:53768/1‑124 (MQ=255)
                                                                                       cAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCt                                                              <  5:383549/149‑1 (MQ=255)
                                                                                          agaTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGAt                                                           >  6:46174/1‑149 (MQ=255)
                                                                                            aTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCt                                                         <  2:390573/149‑1 (MQ=255)
                                                                                                cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGACgg                                                                                      >  6:319153/1‑116 (MQ=255)
                                                                                                cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGACgg                                                                                      <  5:319153/116‑1 (MQ=255)
                                                                                                     gCCTGATAAGCGTAGCGCATCAGGCAAGTTTGCGTTCAAACTCCCCGACAAGCCGGGAGTTTGGAGTAGGTTATGCCCACGCTTGCAGGAAACGCTGCCCGTACTGGTCGGCGACCAGCAGAGCGGCGTACACCTGATCTGGCGTcgcg                                                >  5:77235/1‑149 (MQ=255)
                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTg                                                        <  2:223904/140‑1 (MQ=255)
                                                                                                      ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTg                                                        >  1:223904/1‑140 (MQ=255)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccgc                                            >  6:173301/1‑149 (MQ=255)
                                                                                                         gATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgccgc                                            >  5:385623/1‑149 (MQ=255)
                                                                                                          aTAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTAGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCgccgcc                                           >  4:373989/1‑149 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAg                                         <  4:207814/148‑1 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAg                                         <  3:342082/148‑1 (MQ=255)
                                                                                                             aGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAg                                         <  6:385623/148‑1 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAt                                     >  5:312653/1‑149 (MQ=255)
                                                                                                                gTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAt                                     <  6:312653/149‑1 (MQ=255)
                                                                                                                   gcgcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgtt                                  <  7:403368/149‑1 (MQ=255)
                                                                                                                     gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgttgt                                <  1:69895/149‑1 (MQ=255)
                                                                                                                     gcATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgttgt                                <  2:51671/149‑1 (MQ=255)
                                                                                                                       aTCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGcgcc                                              >  8:403616/1‑133 (MQ=255)
                                                                                                                          aGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAAt                            <  1:217852/148‑1 (MQ=255)
                                                                                                                          aGGCAATTTTGCGCTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATg                           <  5:8472/149‑1 (MQ=255)
                                                                                                                            gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATg                           >  2:259953/1‑147 (MQ=255)
                                                                                                                            gCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATg                           <  1:259953/147‑1 (MQ=255)
                                                                                                                               aTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTc                      <  5:125391/149‑1 (MQ=255)
                                                                                                                               aTTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTc                      <  5:173301/149‑1 (MQ=255)
                                                                                                                                    tactccAAACTCCCGGCTTGTCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCCGCAGCGAGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCtt                 >  6:396279/6‑149 (MQ=255)
                                                                                                                                       ttCAAACTCCCGGACAAGCCGTGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTg              <  4:100000/149‑1 (MQ=255)
                                                                                                                                       ttCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTACACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGGTTCACCTTCTg              >  3:229715/1‑149 (MQ=255)
                                                                                                                                          aaaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGca            >  8:364097/1‑148 (MQ=255)
                                                                                                                                          aaaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTAATGGTCGGCTACCAGCAGATCGGCGTAAACCTAATCATGCGTCGAGCCGCCAGTCATGATGTGAATGGTATCACCTTCTGca            >  8:257587/1‑148 (MQ=255)
                                                                                                                                           aaCTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTTTTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGACGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGgaca          >  5:143103/1‑149 (MQ=255)
                                                                                                                                           aaCTCCCGGACAAGCCGGGAGGTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGcaca          <  2:290830/149‑1 (MQ=255)
                                                                                                                                              tCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGc       >  7:150448/1‑149 (MQ=255)
                                                                                                                                               cccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGcc      >  8:314605/1‑149 (MQ=255)
                                                                                                                                               cccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGcc      >  7:119574/1‑149 (MQ=255)
                                                                                                                                                ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgttgt                                >  6:257063/1‑122 (MQ=255)
                                                                                                                                                ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCAtgttgt                                <  5:257063/122‑1 (MQ=255)
                                                                                                                                                ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGtt                        <  3:139643/130‑1 (MQ=255)
                                                                                                                                                ccGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGtt                        >  4:139643/1‑130 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGGAGTTTCAGATTACTGGCACAGAGTGCCGGATGCGGAGTGAACGCCTTATCCGGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGATTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTCAAACTCCCGGACAAGCCGGGAGTTTGGAGTAGGTTATTCCCACTCTTGCAGGAAACGCTGACCGTACTGGTCGGCTACCAGCAGAGCGGCGTAAACCTGATCTGGCGTCGCGCCGCCAGGCATGTTGTGAATGGTTTCACCTTCTGCACATGCCGCTT  >  NC_000913/4137761‑4138056

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: