Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,972,673 A→C 8.3% E51A (GAG→GCG) ‡ rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6730AC 8.3% 326.6 / 20.5 96E51A (GAG→GCG) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base A (40/48);  new base C (3/5);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 106 total aligned reads displayed.

AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAC  >  NC_000913/3972525‑3972820
                                                                                                                                                    |                                                                                                                                                   
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcga                                                                                                                                                     <  3:164205/149‑1 (MQ=255)
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcga                                                                                                                                                     >  3:529496/1‑149 (MQ=255)
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcga                                                                                                                                                     <  4:501061/149‑1 (MQ=255)
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGCGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcga                                                                                                                                                     >  4:495241/1‑149 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgag                                                                                                                                                    <  2:470010/149‑1 (MQ=255)
  aaaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  4:437726/1‑149 (MQ=255)
  aaaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  2:612740/1‑149 (MQ=255)
  aaaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  1:564842/1‑149 (MQ=255)
  aaaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  4:249310/1‑149 (MQ=255)
       ttCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCttt                                                                                                                                              >  3:62656/1‑149 (MQ=255)
         cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGc                                                                                                                                            <  1:286161/149‑1 (MQ=255)
           gATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                          >  3:228157/1‑149 (MQ=255)
           gATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                          >  2:340838/1‑149 (MQ=255)
            aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCtt                                                                                                                                         <  1:59033/149‑1 (MQ=255)
             tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                                                                                                        <  3:409845/149‑1 (MQ=255)
              aaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       >  2:588944/1‑149 (MQ=255)
              aaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       <  1:156311/149‑1 (MQ=255)
              aaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       >  3:269487/1‑149 (MQ=255)
               aCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                      <  4:62656/149‑1 (MQ=255)
               aCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                      <  2:604935/149‑1 (MQ=255)
                  ggtggtGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGaa                                                                                                                                   >  3:200972/1‑149 (MQ=255)
                   gtggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGaaa                                                                                                                                  <  2:564842/149‑1 (MQ=255)
                           gggTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGatat                                                                                                                          >  4:281213/1‑149 (MQ=255)
                             gTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCttt                                                                                                                                              <  3:345363/127‑1 (MQ=255)
                             gTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCttt                                                                                                                                              >  4:345363/1‑127 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAGCGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCCTTGAGAAAGTTGATATCTGCGATCg                                                                                                                 >  1:509646/1‑149 (MQ=255)
                                     gCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCgg                                                                                                                <  1:395642/149‑1 (MQ=255)
                                         ggCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAg                                                                                                            >  4:390599/1‑149 (MQ=255)
                                           cgcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGaa                                                                                                          <  4:247024/149‑1 (MQ=255)
                                           cgcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGaa                                                                                                          <  1:376750/149‑1 (MQ=255)
                                              tGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTg                                                                                                       <  3:245268/149‑1 (MQ=255)
                                                     ttATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                >  2:38500/1‑149 (MQ=255)
                                                       atatCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCtt                                                                                                                                               >  3:9369/1‑100 (MQ=255)
                                                        tatCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTAt                                                                                             >  1:521186/1‑149 (MQ=255)
                                                            atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGATCTGGCAcgc                                                                                                 <  1:536097/141‑1 (MQ=255)
                                                            atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTgg                                                                                                      >  3:470217/1‑136 (MQ=255)
                                                            atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTgg                                                                                                      <  4:470217/136‑1 (MQ=255)
                                                            atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgc                                                                                                 >  2:536097/1‑141 (MQ=255)
                                                            atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc                                                                                         >  3:456018/1‑149 (MQ=255)
                                                             tcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGa                                                                                                                             <  1:609889/112‑1 (MQ=255)
                                                             tcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGa                                                                                                                             >  2:609889/1‑112 (MQ=255)
                                                             tcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACt                                                                                        >  1:680/1‑149 (MQ=255)
                                                              caACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTg                                                                                       <  3:425754/149‑1 (MQ=255)
                                                               aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGa                                                                                      >  4:500776/1‑149 (MQ=255)
                                                               aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGa                                                                                      <  2:538134/149‑1 (MQ=255)
                                                               aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGa                                                                                      <  1:245662/149‑1 (MQ=255)
                                                                acgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTGCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCggg                                                                                                               >  3:379717/1‑123 (MQ=255)
                                                                acgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTCGCCTTTGAGAAAGTTGATATCTGCGATCggg                                                                                                               <  4:379717/123‑1 (MQ=255)
                                                                  gaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCa                                                                                   >  1:318922/1‑149 (MQ=255)
                                                                     acgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa                                                                                <  4:529496/149‑1 (MQ=255)
                                                                      cgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAg                                                                               <  1:283788/149‑1 (MQ=255)
                                                                       gaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGc                                                                              >  2:483467/1‑149 (MQ=255)
                                                                         gCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCa                                                                            <  4:200972/149‑1 (MQ=255)
                                                                            aCGCGGTGGTAGTGGTTGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGAc                                                                         <  3:33664/149‑1 (MQ=255)
                                                                            aCGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGAc                                                                         >  4:569698/1‑149 (MQ=255)
                                                                               cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtgt                                                                      <  1:14476/149‑1 (MQ=255)
                                                                               cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtgt                                                                      <  3:390599/149‑1 (MQ=255)
                                                                               cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtgt                                                                      <  2:318922/149‑1 (MQ=255)
                                                                                 gtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGa                                                                                      >  4:454830/1‑131 (MQ=255)
                                                                                 gtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGa                                                                                      <  3:454830/131‑1 (MQ=255)
                                                                                 gtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTgtgt                                                                    <  3:397818/149‑1 (MQ=255)
                                                                                 gtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTgtgt                                                                    >  4:100424/1‑149 (MQ=255)
                                                                                  tggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTc                                                                   <  2:230589/149‑1 (MQ=255)
                                                                                    gtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAt                                                                 <  1:588944/149‑1 (MQ=255)
                                                                                          gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCt                                                           >  3:334658/1‑149 (MQ=255)
                                                                                          gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCt                                                           >  3:360690/1‑149 (MQ=255)
                                                                                                  gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                   >  3:344834/1‑149 (MQ=255)
                                                                                                     gACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaaa                                                <  3:495241/149‑1 (MQ=255)
                                                                                                        cTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGa                                                                                                                             >  1:368120/1‑69 (MQ=255)
                                                                                                        cTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGa                                                                                                                             <  2:368120/69‑1 (MQ=255)
                                                                                                          aCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                           <  3:249310/149‑1 (MQ=255)
                                                                                                                   aCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGt                                  <  4:228157/149‑1 (MQ=255)
                                                                                                                     cTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                <  4:76892/149‑1 (MQ=255)
                                                                                                                     cTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                >  3:434347/1‑149 (MQ=255)
                                                                                                                     cTGATGTCGCTGGCACCGGTCGCGCAAAGAGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCCGACTGTGTCATGCATCTGGCAGCAGCAAGCCATGTTGACCGTTc                                >  2:482896/1‑149 (MQ=255)
                                                                                                                      tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCt                               <  3:100424/149‑1 (MQ=255)
                                                                                                                      tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCt                               >  4:584343/1‑149 (MQ=255)
                                                                                                                         tGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAtt                            <  3:281213/149‑1 (MQ=255)
                                                                                                                             gCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAATGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACg                        >  2:594972/1‑149 (MQ=255)
                                                                                                                               tGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGc                      <  3:500776/149‑1 (MQ=255)
                                                                                                                                 gCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGccc                    <  2:521186/149‑1 (MQ=255)
                                                                                                                                   aCCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGAc                         <  3:534368/142‑1 (MQ=255)
                                                                                                                                   aCCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGAc                         >  4:534368/1‑142 (MQ=255)
                                                                                                                                      tcaaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGa                                      >  4:130282/7‑126 (MQ=255)
                                                                                                                                      tcaaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGa                                      <  3:130282/120‑1 (MQ=255)
                                                                                                                                       gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGc              >  1:6239/1‑149 (MQ=255)
                                                                                                                                          gcgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt           >  1:478669/1‑149 (MQ=255)
                                                                                                                                          aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaaa                                                <  1:111548/110‑1 (MQ=255)
                                                                                                                                          aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaaa                                                >  2:111548/3‑112 (MQ=255)
                                                                                                                                          aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaaa                                                <  1:528451/110‑1 (MQ=255)
                                                                                                                                          aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCAGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaaa                                                >  2:528451/3‑112 (MQ=255)
                                                                                                                                             cAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAt        <  1:612740/149‑1 (MQ=255)
                                                                                                                                             cAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAt        <  2:117946/149‑1 (MQ=255)
                                                                                                                                             cAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAt        <  4:269487/149‑1 (MQ=255)
                                                                                                                                              aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAtt       >  1:456839/1‑149 (MQ=255)
                                                                                                                                              aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAtt       >  1:158418/1‑149 (MQ=255)
                                                                                                                                              aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAtt       <  2:6239/149‑1 (MQ=255)
                                                                                                                                              aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAtt       >  3:348360/1‑149 (MQ=255)
                                                                                                                                               aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTg      <  3:274228/149‑1 (MQ=255)
                                                                                                                                                   gagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAc  >  4:360925/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAC  >  NC_000913/3972525‑3972820

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