Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,981,814 T→A 6.2% T192S (ACG→TCG)  otsB ← trehalose‑6‑phosphate phosphatase, biosynthetic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,981,8140TA 6.2% 408.3 / 14.3 113T192S (ACG→TCG) otsBtrehalose‑6‑phosphate phosphatase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base T (52/54);  new base A (3/4);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-02

Only 100 of 121 total aligned reads displayed.

AAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAA  >  NC_000913/1981668‑1981961
                                                                                                                                                  |                                                                                                                                                   
aaGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGt                                                                                                                                                     >  2:222272/1‑147 (MQ=255)
aaGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGt                                                                                                                                                     <  1:222272/147‑1 (MQ=255)
aaGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTc                                                                                                                                                   <  3:745290/149‑1 (MQ=255)
 aGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGTTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCg                                                                                                                                                  >  2:255944/1‑149 (MQ=255)
     aGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCg                                                                                                                                              >  4:515212/1‑149 (MQ=255)
          cAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATaaa                                                                                                                                         >  1:576583/1‑149 (MQ=255)
          cAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATaaa                                                                                                                                         <  3:202894/149‑1 (MQ=255)
          cAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATaaa                                                                                                                                         >  4:744469/1‑149 (MQ=255)
               aTCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAg                                                                                                                                    <  1:142615/149‑1 (MQ=255)
                     cacaCCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCt                                                                                                                              <  2:576583/149‑1 (MQ=255)
                        acCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGtt                                                                                                                                                    <  2:863472/124‑1 (MQ=255)
                        acCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGtt                                                                                                                                                    >  2:554147/1‑124 (MQ=255)
                        acCCGCCAGTCGCCATGATGCATGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGtt                                                                                                                                                    <  1:554147/124‑1 (MQ=255)
                               aGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGc                                                                                                                    >  4:259182/1‑149 (MQ=255)
                                 tCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgt                                                                                                                                             >  2:622154/1‑121 (MQ=255)
                                 tCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgt                                                                                                                                             <  1:622154/122‑2 (MQ=255)
                                    ccATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCaa                                                                                                               >  4:431039/1‑149 (MQ=255)
                                     cATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAAt                                                                                                              >  4:660795/1‑149 (MQ=255)
                                      atgatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGc                                                                                                                    <  2:452597/142‑1 (MQ=255)
                                      atgatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGc                                                                                                                    >  1:452597/1‑142 (MQ=255)
                                         atgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGcc                                                                                                          <  2:1317154/149‑1 (MQ=255)
                                             cTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      <  2:556902/149‑1 (MQ=255)
                                              tGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAcc                                                                                                     <  3:259182/149‑1 (MQ=255)
                                              tGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      <  3:921969/148‑1 (MQ=255)
                                              tGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      >  4:921969/1‑148 (MQ=255)
                                              tGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      <  1:1037531/148‑1 (MQ=255)
                                              tGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      >  2:1037531/1‑148 (MQ=255)
                                              tGAGTTGCACCGGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAcc                                                                                                     >  1:1396705/1‑149 (MQ=255)
                                                  ttGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGccc                                                                                                                                               <  3:1019912/103‑1 (MQ=255)
                                                  ttGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGccc                                                                                                                                               >  4:1019912/1‑103 (MQ=255)
                                                  ttGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTa                                                                                                 <  4:16693/149‑1 (MQ=255)
                                                   tGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGc                                                                                                                            <  3:1846924/121‑1 (MQ=255)
                                                   tGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGc                                                                                                                            >  4:1846924/1‑121 (MQ=255)
                                                   tGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTAc                                                                                                <  3:1357269/149‑1 (MQ=255)
                                                    gCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      <  3:1705438/142‑1 (MQ=255)
                                                    gCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAc                                                                                                      >  4:1705438/1‑142 (MQ=255)
                                                    gCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACt                                                                                               <  4:1428745/149‑1 (MQ=255)
                                                    gCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACt                                                                                               <  3:431039/149‑1 (MQ=255)
                                                           tgCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACct                                                                                        >  1:49792/1‑149 (MQ=255)
                                                             ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGa                                                                                                                                     >  3:453260/1‑102 (MQ=255)
                                                             ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGa                                                                                                                                     <  4:453260/102‑1 (MQ=255)
                                                             ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACctct                                                                                      >  4:714111/1‑149 (MQ=255)
                                                             ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACctct                                                                                      >  2:996875/1‑149 (MQ=255)
                                                              cAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                      <  4:1605223/116‑1 (MQ=255)
                                                              cAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                      >  3:1605223/1‑116 (MQ=255)
                                                              cAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCaa                                                                                                               >  4:1463822/1‑123 (MQ=255)
                                                              cAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCaa                                                                                                               <  3:1463822/123‑1 (MQ=255)
                                                                aTTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCgg                                                                                   <  3:714111/149‑1 (MQ=255)
                                                                       cTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtc                                                                            <  3:660795/149‑1 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAg                                                                                                                                    <  2:224191/92‑1 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAg                                                                                                                                    >  1:224191/1‑92 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATa                                                                                                                         <  1:833169/103‑1 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAt                                                                             >  3:305164/1‑147 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAt                                                                             <  4:305164/147‑1 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtat                                                                                                                                           <  1:350756/85‑4 (MQ=255)
                                                                        tGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtat                                                                                                                                           >  2:350756/1‑82 (MQ=255)
                                                                         gACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtctc                                                                          <  1:996875/149‑1 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGt                                                                                  >  3:1106891/1‑140 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGt                                                                                  <  4:1106891/140‑1 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtctc                                                                          >  1:115791/1‑148 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtctc                                                                          <  2:115791/148‑1 (MQ=255)
                                                                            aTTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGac                                                                       >  1:1616621/1‑149 (MQ=255)
                                                                             ttCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgt                                                                                                                                             <  3:665605/78‑2 (MQ=255)
                                                                                              aCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTa                                                                                                 <  1:1027805/105‑1 (MQ=255)
                                                                                              aCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTa                                                                                                 >  2:1027805/1‑105 (MQ=255)
                                                                                              aCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaaca                                                                   <  2:1849045/135‑1 (MQ=255)
                                                                                              aCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaaca                                                                   >  1:1849045/1‑135 (MQ=255)
                                                                                                      gAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                           >  3:58003/1‑87 (MQ=255)
                                                                                                      gAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                           <  4:58003/87‑1 (MQ=255)
                                                                                                      gAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtct                                                                           <  3:1566788/119‑1 (MQ=255)
                                                                                                      gAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtct                                                                           >  4:1566788/1‑119 (MQ=255)
                                                                                                         gCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCAttt                                          <  2:1396705/149‑1 (MQ=255)
                                                                                                         gCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCAttt                                          >  1:688916/1‑149 (MQ=255)
                                                                                                            aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaa                                                                     >  3:1638475/1‑119 (MQ=255)
                                                                                                            aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaa                                                                     <  4:1638475/119‑1 (MQ=255)
                                                                                                            aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTaa                                                  >  1:129568/1‑138 (MQ=255)
                                                                                                            aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTaa                                                  <  2:129568/138‑1 (MQ=255)
                                                                                                            aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCAt                                            >  3:1624698/1‑144 (MQ=255)
                                                                                                               ttCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgtg                                       <  1:1877237/146‑1 (MQ=255)
                                                                                                               ttCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgtg                                       >  2:1877237/1‑146 (MQ=255)
                                                                                                                  aTCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACg                                                >  3:1782096/1‑134 (MQ=255)
                                                                                                                   tCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                           <  3:328836/74‑1 (MQ=255)
                                                                                                                   tCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                           >  4:328836/1‑74 (MQ=255)
                                                                                                                   tCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGacaa                                                                     <  2:1271108/112‑1 (MQ=255)
                                                                                                                     ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCt                              >  4:1657779/1‑149 (MQ=255)
                                                                                                                     ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCt                              >  4:709434/1‑149 (MQ=255)
                                                                                                                      gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                             >  2:304855/1‑149 (MQ=255)
                                                                                                                      gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                             >  2:1481139/1‑149 (MQ=255)
                                                                                                                        tAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAg                           <  1:304855/149‑1 (MQ=255)
                                                                                                                         aaaTCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAt                                                                             >  1:1225581/1‑98 (MQ=255)
                                                                                                                             catcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCCCCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTCCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATa                      <  4:1741726/149‑1 (MQ=255)
                                                                                                                               tcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACg                    >  4:1075703/1‑149 (MQ=255)
                                                                                                                                    cAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATc                               >  1:1062551/1‑133 (MQ=255)
                                                                                                                                    cAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATc                               <  2:1062551/133‑1 (MQ=255)
                                                                                                                                     aGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGa                                                                              <  1:251252/85‑1 (MQ=255)
                                                                                                                                         aTACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAAt          <  3:709434/149‑1 (MQ=255)
                                                                                                                                          tACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATg         >  4:256473/1‑149 (MQ=255)
                                                                                                                                               gCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCAtt    <  2:688916/149‑1 (MQ=255)
                                                                                                                                                  ttCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCctgc                                                       <  2:1876407/95‑1 (MQ=255)
                                                                                                                                                  ttCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCctgc                                                       >  1:1876407/1‑95 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCATTAA  >  NC_000913/1981668‑1981961

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: