Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,981,815 T→A 5.7% R191R (CGA→CGT 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,8150TA 5.7% 373.4 / 9.8 106R191R (CGA→CGTotsBtrehalose‑6‑phosphate phosphatase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base T (46/54);  new base A (2/4);  total (48/58)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 106 total aligned reads displayed.

AAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCAT  >  NC_000913/1981668‑1981958
                                                                                                                                                    |                                                                                                                                               
aaGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTc                                                                                                                                                <  3:481164/149‑1 (MQ=255)
aaGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTc                                                                                                                                                >  2:788496/1‑149 (MQ=255)
aaGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTc                                                                                                                                                >  1:309754/1‑149 (MQ=255)
  gCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGc                                                                                                                                              >  2:288497/1‑149 (MQ=255)
    cAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGccc                                                                                                                                            <  4:475991/149‑1 (MQ=255)
    cAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGccc                                                                                                                                            >  3:97738/1‑149 (MQ=255)
      gCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGa                                                                                                                                          <  4:454500/149‑1 (MQ=255)
       cTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGAt                                                                                                                                         <  3:247306/149‑1 (MQ=255)
       cTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgta                                                                                                                                         <  4:870016/149‑3 (MQ=255)
        tCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtaa                                                                                                                                        <  2:793994/149‑4 (MQ=255)
         ccAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtatt                                                                                                                                       <  3:793868/149‑5 (MQ=255)
           aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAg                                                                                                                                     >  3:476964/1‑149 (MQ=255)
           aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAg                                                                                                                                     >  1:309157/1‑149 (MQ=255)
           aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAg                                                                                                                                     >  1:190100/1‑149 (MQ=255)
           aGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAg                                                                                                                                     <  1:936922/149‑1 (MQ=255)
            gACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGAACGCCCgtatttct                                                                                                                                    <  4:773293/149‑8 (MQ=255)
             aCATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAggg                                                                                                                                   <  4:684314/149‑1 (MQ=255)
               aTCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAg                                                                                                                                 >  1:692796/1‑149 (MQ=255)
               aTCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAg                                                                                                                                 >  1:800489/1‑149 (MQ=255)
                   ggCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTc                                                                                                                             <  1:788496/149‑1 (MQ=255)
                    gCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTcc                                                                                                                            >  1:863786/1‑149 (MQ=255)
                     cacaCCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCt                                                                                                                           <  4:887616/149‑1 (MQ=255)
                     cacaCCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCt                                                                                                                           <  2:294774/149‑1 (MQ=255)
                     cacaCCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCt                                                                                                                           <  3:349095/149‑1 (MQ=255)
                             ccAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                   <  3:273367/149‑1 (MQ=255)
                             ccAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                   <  1:647320/149‑1 (MQ=255)
                                gTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                   <  1:732303/146‑1 (MQ=255)
                                gTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                   >  2:732303/1‑146 (MQ=255)
                                gTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCt                                                                                                                <  4:249550/149‑1 (MQ=255)
                                gTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCt                                                                                                                <  2:73181/149‑1 (MQ=255)
                                  cGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGc                                                                                                              <  4:862919/149‑1 (MQ=255)
                                   gCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCa                                                                                                             <  4:270832/149‑1 (MQ=255)
                                    ccATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                   >  3:854405/1‑142 (MQ=255)
                                    ccATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATaaaa                                                                                                                   <  4:854405/142‑1 (MQ=255)
                                    ccATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCaa                                                                                                            <  3:26608/149‑1 (MQ=255)
                                     cATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGACATACGGGCGTTCGCCCGATAAAGGGAGCt                                                                                                                               <  3:793605/129‑1 (MQ=255)
                                     cATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGCGCt                                                                                                                               >  2:474568/1‑129 (MQ=255)
                                     cATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCt                                                                                                                               <  1:474568/129‑1 (MQ=255)
                                     cATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCt                                                                                                                               >  4:793605/1‑129 (MQ=255)
                                        gatgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                        <  1:787211/149‑1 (MQ=255)
                                              tGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCAcc                                                                                                  >  2:447/1‑149 (MQ=255)
                                                      aCCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTgg                                                                                          >  2:60333/1‑149 (MQ=255)
                                                       ccTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGt                                                                                         <  4:122958/149‑1 (MQ=255)
                                                       ccTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGt                                                                                         <  3:563692/149‑1 (MQ=255)
                                                        cTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTa                                                                                        >  1:24058/1‑149 (MQ=255)
                                                             ccAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACctct                                                                                   >  1:212500/1‑149 (MQ=255)
                                                              cAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCtctc                                                                                  >  1:889827/1‑149 (MQ=255)
                                                               aaTTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAc                                                                                       >  2:816716/1‑143 (MQ=255)
                                                               aaTTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTAc                                                                                       <  1:816716/143‑1 (MQ=255)
                                                               aaTTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCg                                                                                 <  3:82749/149‑1 (MQ=255)
                                                                 tttttACTGACATTCCGCCGAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGt                                                                               <  4:630615/149‑1 (MQ=255)
                                                                 tttttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGt                                                                               <  1:48143/149‑1 (MQ=255)
                                                                  ttttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGtt                                                                              <  3:335085/149‑1 (MQ=255)
                                                                   tttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGttt                                                                             <  2:591740/149‑1 (MQ=255)
                                                                   tttACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGttt                                                                             >  3:38639/1‑149 (MQ=255)
                                                                     tACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGa                                                                           <  2:622563/149‑1 (MQ=255)
                                                                      aCTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAt                                                                          >  4:62906/1‑149 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCg                                                                      >  4:438523/1‑149 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCg                                                                      >  3:839440/1‑149 (MQ=255)
                                                                          aCATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCg                                                                      >  1:309875/1‑149 (MQ=255)
                                                                           cATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGa                                                                     <  2:730334/149‑1 (MQ=255)
                                                                             ttCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGaca                                                                   >  1:725020/1‑149 (MQ=255)
                                                                                 gCCCAGTCGGTTAACGACTGCGAAGCGAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAacac                                                               >  2:134636/1‑149 (MQ=255)
                                                                                      agtcgTTAACCGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACtt                                                           >  3:886341/5‑149 (MQ=255)
                                                                                        cGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTACCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTcc                                                        <  4:531141/149‑1 (MQ=255)
                                                                                          gTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtctc                                                                       <  4:652298/132‑1 (MQ=255)
                                                                                          gTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGAtctc                                                                       >  3:652298/1‑132 (MQ=255)
                                                                                              aCGACAGCGAAGCCAGAATCATCGGTTGAATCATCGCCCAGAAATACGGGCGTTCGCCCGATGACGGGAGTTTCCTGCACAAAAGCTGCAATTGCCGCACCTTTACTCGTACCTCTCGGTTTGATCTCGACAACACACTTTACctgctg                                                  >  2:671899/1‑149 (MQ=255)
                                                                                                 gACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTa                                                >  3:210010/1‑149 (MQ=255)
                                                                                                      cGAAGCCAGATTCATCGGTTAAATCATCGCCCAGTAATACGGGCGTTCGCCCGATAAAGGGAGCTTCTTGCATAAATGCTGCAATTGCCTCACCATTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGcc                                           <  4:925298/149‑1 (MQ=255)
                                                                                                      cGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGcc                                           >  1:98425/1‑149 (MQ=255)
                                                                                                          gCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCAttt                                       <  2:725020/149‑1 (MQ=255)
                                                                                                             aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgtg                                    <  2:692796/149‑1 (MQ=255)
                                                                                                             aGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTtgtg                                    <  2:596969/149‑1 (MQ=255)
                                                                                                              gATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAGGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTgg                                   <  2:309754/149‑1 (MQ=255)
                                                                                                                 tCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCa                                <  2:98425/149‑1 (MQ=255)
                                                                                                                   aTCGGTTAAATCATCGCCCAGCAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCGCACCTTTACTGGTACCTCTCGGTTTGATCTGGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGa                              >  1:670583/1‑149 (MQ=255)
                                                                                                                   aTCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGa                              <  2:309165/149‑1 (MQ=255)
                                                                                                                      ggTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCt                           >  2:145823/1‑149 (MQ=255)
                                                                                                                       gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                          <  4:886341/149‑1 (MQ=255)
                                                                                                                       gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                          >  2:806800/1‑149 (MQ=255)
                                                                                                                       gTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTg                          >  4:394869/1‑149 (MQ=255)
                                                                                                                        ttAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGa                         <  1:288497/149‑1 (MQ=255)
                                                                                                                              catcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATa                   <  2:863786/149‑1 (MQ=255)
                                                                                                                                tcGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACg                 >  1:741487/1‑149 (MQ=255)
                                                                                                                                  gCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGtt               <  4:599901/149‑1 (MQ=255)
                                                                                                                                   cccAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTg              <  2:497724/149‑1 (MQ=255)
                                                                                                                                    ccAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                        <  1:713177/58‑1 (MQ=255)
                                                                                                                                    ccAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgc             >  4:650616/1‑149 (MQ=255)
                                                                                                                                    ccAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgc             >  3:119018/1‑149 (MQ=255)
                                                                                                                                    caaGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGc                                                                                                        >  2:713177/3‑58 (MQ=255)
                                                                                                                                     cAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgcg            >  3:219886/1‑149 (MQ=255)
                                                                                                                                      aGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTgcgc           <  2:309157/149‑1 (MQ=255)
                                                                                                                                       gAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGac                                                                    <  1:827930/91‑1 (MQ=255)
                                                                                                                                       gAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGac                                                                    >  2:827930/1‑91 (MQ=255)
                                                                                                                                       ctttatcGGGCGAACGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCt          >  1:368181/7‑149 (MQ=255)
                                                                                                                                           tACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATg      <  2:348724/149‑1 (MQ=255)
                                                                                                                                           atcGGGCGAACGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATg      >  4:200647/3‑149 (MQ=255)
                                                                                                                                              gggCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCa   >  3:317647/1‑149 (MQ=255)
                                                                                                                                               ggCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCAt  <  4:775854/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAACGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGGCGTTCGCCCGATAAAGGGAGCTTCCTGCATAAAAGCTGCAATTGCCTCACCTTTACTGGTACCTCTCGGTTTGATCTCGACAACACACTTTCCCTGCTGTAACGCCATTTGTGGCCAGATCTGAGTAATACGTTGCGCTAATGTCAT  >  NC_000913/1981668‑1981958

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