Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,749,284 T→C 5.1% T254A (ACT→GCT)  U712_19325 ← Hydroxyethylthiazole kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,749,2840TC 5.1% 252.3 / 3.1 79T254A (ACT→GCT) U712_19325Hydroxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (48/27);  new base C (2/2);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ATCGGCTTTCGTATTGTTTGACCCCATAATGAAGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCCACAGCGC  >  CP006881/3749148‑3749421
                                                                                                                                        |                                                                                                                                         
aTCGGCTTTCGTATTGTTTGACCCCATAATGAAGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCAGTCTTGGCCGTCTTGTTCAGGTACAGTTg                                                                                                                                         >  2:210190/1‑139 (MQ=255)
   ggCTTTCGTATTGTTTGACCCCATAATGAAGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGa                                                                                                                                        >  5:145668/1‑137 (MQ=255)
                          tAATGAAGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTg                                                                                                                >  6:262885/1‑138 (MQ=255)
                          tAATGAAGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCAATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCCGCAAGTCAATCTg                                                                                                                >  8:43935/1‑138 (MQ=255)
                                aGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTg                                                                                                                                         >  8:277136/1‑107 (MQ=255)
                                aGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTg                                                                                                                                         <  7:277136/107‑1 (MQ=255)
                                 gTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTc                                                                                                        >  6:240345/1‑139 (MQ=255)
                                  tataCAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTc                                                                                                        >  4:121347/1‑138 (MQ=255)
                                    taCAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAg                                                                                                     >  5:86262/1‑139 (MQ=255)
                                    taCAGATAATAACTCTTTAATCATTTCCCGAGAAATACTAGTCATGAGACAGTCACACTTTATATAGTCTATCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGATTGTTCATAAATTCAATCTGACAGCTTCCAg                                                                                                     >  1:9285/1‑139 (MQ=255)
                                     aCAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGCTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAgg                                                                                                    >  8:72881/1‑139 (MQ=255)
                                     aCAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGCTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAgg                                                                                                    >  8:68108/1‑139 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACCGTCACCGTTTCTATAGTCGCCCATTCTTGGACGTCGTGTTAAGTAACAGTTGCAAGCTTGTTCAGCAATTGAATCTGAACGCTTCCAgtgact                                                                                                >  4:290379/1‑134 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAAtt                                                                                                                       <  3:86525/116‑1 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAAtt                                                                                                                       >  4:86525/1‑116 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCt                                                                                                >  3:71951/1‑139 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCt                                                                                                >  3:37381/1‑139 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCt                                                                                                >  5:179392/1‑139 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTAAGTAACAGTTTAAATCTTTTTAAGCAATTAAATCAGAAATCATCCAGGGCCt                                                                                                >  3:273887/1‑139 (MQ=255)
                                         ataataACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCATTTCTATAGTCGCGCAGTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTAGTTCAGCAATTCAATCTGATAGAGTCCGGGGCCt                                                                                                >  4:290381/1‑139 (MQ=255)
                                          taataaCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCtt                                                                                               >  7:152841/1‑139 (MQ=255)
                                             taaCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGt                                                                                            <  5:262885/139‑1 (MQ=255)
                                             taaCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGCTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGt                                                                                            <  7:68108/139‑1 (MQ=255)
                                             taaCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGCTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGt                                                                                            <  7:72881/139‑1 (MQ=255)
                                              aaCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAAtt                                                                                                                       <  4:201075/111‑1 (MQ=255)
                                              aaCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAAtt                                                                                                                       >  3:201075/1‑111 (MQ=255)
                                               aCTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCa                                                                                          <  3:121347/139‑1 (MQ=255)
                                                   tttCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATACAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCg                                                                                      >  6:49005/1‑139 (MQ=255)
                                                     tcatcaTTTCCCGAGAAATACGAGTCATGAGACAGTCCCCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTTAAAGCTTCCAGGGCCTTTGTCAGCCGTc                                                                                    >  5:234295/1‑139 (MQ=255)
                                                     tcatcaTTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCGGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGcc                                                                                       <  4:108469/136‑1 (MQ=255)
                                                          atTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTctgct                                                                                <  1:103297/138‑1 (MQ=255)
                                                           tttCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCtgt                                                                              >  7:32810/1‑139 (MQ=255)
                                                               ccGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCg                                                                           <  1:137127/138‑1 (MQ=255)
                                                                           aGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAgcgg                                                              <  8:32810/139‑1 (MQ=255)
                                                                                tGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGAcc                                                         <  4:71951/139‑1 (MQ=255)
                                                                                      aGTCACCCTTTCTATAGTCGCCCATTCTTGGCCGTCTTGTTCAGTAACAGTTTAAAGCTTGTTCAGCAATTCAAACTGAAAGCTTCCAGGGCCGTTGTCAGCCGTCTGCTGTGCGGCCAGCTGAGCGGCGCCCCCATAc                                                   >  8:267372/1‑139 (MQ=255)
                                                                                      aGTCACCCTATCTAAAGTCCCCCAATCTTGGACGTCATGTTCAGTAACAGTTGAAAGCATGTTCAGCAATTCAAACAGAAAGCTTCCAGGGCCTATGTCAGCCGACAGCTGTCCGGCAAGCAGAGCCGCGACCCCAAAc                                                   >  7:4108/1‑139 (MQ=255)
                                                                                       gTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATAc                                                   >  5:54379/1‑138 (MQ=255)
                                                                                                tCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATgg                                         >  8:93494/1‑139 (MQ=255)
                                                                                                    ttttCGCCCATTATTGGACGTCTTGTTCAGTAACTGTTGAAAGCTTGTTCTGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGc                                        <  5:34015/133‑1 (MQ=255)
                                                                                                    taGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGc                                        <  6:54379/136‑1 (MQ=255)
                                                                                                    taGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGc                                     >  1:129836/1‑139 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAATCTTGTTCAGCAATTCAATCAGAAAGCTTACAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTTAGCGGCGCCCCCATACGAAGAAAAGGCCGCa                                    >  2:88552/1‑138 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCCGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGt                                                                                            >  6:63246/1‑82 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGt                                                                                            <  5:63246/82‑1 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGc                                                                         <  7:114989/101‑1 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTTAGCGGCGACCCCATACGAAGAAATGGCCGCa                                    >  1:152270/1‑138 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGGCTGCTGTGCGGc                                                                         >  8:114989/1‑101 (MQ=255)
                                                                                                      gTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATATGAAAGCTTCCAGGGCCTTTGTCAGCCCTCTGCTGTGCGGCCAGCTGAGCGGCGACCCCAAACGAAGAAATTGCCGCa                                    >  1:154441/1‑138 (MQ=255)
                                                                                                       tCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTTAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAAt                                  <  6:234295/139‑1 (MQ=255)
                                                                                                        cGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCa                                                                                          >  8:48376/1‑82 (MQ=255)
                                                                                                        cGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCa                                                                                          <  7:48376/82‑1 (MQ=255)
                                                                                                         gCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCAGACGAAGAAATGGCCGCAATag                                >  6:81208/1‑139 (MQ=255)
                                                                                                         gCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTAAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATag                                >  5:235087/1‑139 (MQ=255)
                                                                                                         acccATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCAAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCACTACGAAGAAATGGCCGCAATag                                >  3:162638/2‑139 (MQ=255)
                                                                                                          cccATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAAtt                                                                                                                       <  4:32698/51‑1 (MQ=255)
                                                                                                          cccATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAAtt                                                                                                                       >  3:32698/1‑51 (MQ=255)
                                                                                                          cccAATATTGGCCGTATTGTTAAGTAAAACTTGAAAGCTTGTTCAGCAATTAAATATTAAAGCTTAAAGTGCCTTTTTAAGCCGAATTCTGTGCGGCAAGATGAGCGGCGACCCCATACGAAGAAAAGGCCACAAAag                                >  6:202212/1‑138 (MQ=255)
                                                                                                            cATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAgcgg                                                              >  1:68011/1‑106 (MQ=255)
                                                                                                            cATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAgcgg                                                              <  2:68011/106‑1 (MQ=255)
                                                                                                                       gTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGa                  <  5:240345/139‑1 (MQ=255)
                                                                                                                       gTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGa                  <  2:129836/139‑1 (MQ=255)
                                                                                                                       gTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGa                  <  3:143089/139‑1 (MQ=255)
                                                                                                                         cTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGaa                                                                                                              <  6:75331/45‑1 (MQ=255)
                                                                                                                         cTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGaa                                                                                                              >  5:75331/1‑45 (MQ=255)
                                                                                                                         cTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGa                  >  7:85886/1‑137 (MQ=255)
                                                                                                                         cTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGa                  <  8:85886/137‑1 (MQ=255)
                                                                                                                         cTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGAtt                >  6:147827/1‑139 (MQ=255)
                                                                                                                           tGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGAtttt              >  4:96867/1‑139 (MQ=255)
                                                                                                                             ttCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGCCCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTttct            >  8:33558/1‑139 (MQ=255)
                                                                                                                                 gTAACAGTTTAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCca        <  7:267372/139‑1 (MQ=255)
                                                                                                                                 gTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCca        <  1:88552/139‑1 (MQ=255)
                                                                                                                                 gTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCca        <  1:15359/139‑1 (MQ=255)
                                                                                                                                  tAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGCAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGCCCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCcac       >  7:50499/1‑139 (MQ=255)
                                                                                                                                  tAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGAGTTTCTTCcac       >  4:92196/1‑139 (MQ=255)
                                                                                                                                  tAACACTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCca        >  8:58296/1‑138 (MQ=255)
                                                                                                                                       gTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCca        <  3:171213/133‑1 (MQ=255)
                                                                                                                                       gTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCca        >  4:171213/1‑133 (MQ=255)
                                                                                                                                       gTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCCACAgcg   >  4:210382/1‑138 (MQ=255)
                                                                                                                                        ttGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCCACAgcgc  >  1:185626/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                         
ATCGGCTTTCGTATTGTTTGACCCCATAATGAAGTATACAGATAATAACTCTTTCATCATTTCCCGAGAAATACGAGTCATGAGACAGTCACCCTTTCTATAGTCGCCCATTCTTGGACGTCTTGTTCAGTAACAGTTGAAAGCTTGTTCAGCAATTCAATCTGAAAGCTTCCAGGGCCTTTGTCAGCCGTCTGCTGTGCGGCAAGCTGAGCGGCGACCCCATACGAAGAAATGGCCGCAATAGCAGCAAACAATGGATTTTCTTCCACAGCGC  >  CP006881/3749148‑3749421

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