Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,307,299 G→A 7.8% P168S (CCC→TCC)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,307,2990GA 7.8% 156.9 / 6.1 65P168S (CCC→TCC) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (3/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-02

GGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGC  >  NC_000913/4307151‑4307447
                                                                                                                                                    |                                                                                                                                                    
ggggggTTTTCCGCATGTACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTCGTTCGTACCAACGTCTTAGCGCTATTCCTGATCAATTCGTTTGCAGGTAGCCAGGAGTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                      <  1:231355/149‑1 (MQ=255)
ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTCCCAGCGTCTTAGCGCCATGCCAGGGCAATTGTTTGCCAGGCACGCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                      <  1:231360/149‑1 (MQ=255)
ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                      <  1:482144/149‑1 (MQ=255)
 gggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAGTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGGGTTGGGTCATATGCCCCCggggg                                                                                                                                                     >  2:115045/1‑149 (MQ=255)
  ggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGa                                                                                                                                                    >  1:866922/1‑149 (MQ=255)
   gCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGTCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCAAGGATTGCCACACGCGCAGGGTTGGGTCATATCTCCCAGGGGGat                                                                                                                                                   >  1:353198/1‑149 (MQ=255)
    cGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggggtg                                                                                                                                                  >  1:911593/1‑146 (MQ=255)
     gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATCCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGatat                                                                                                                                                 <  2:347063/149‑1 (MQ=255)
       tttCCGCATGGACGAAAAGATCGCGTAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                               <  2:495356/149‑1 (MQ=255)
       tttCCGCATGGACGAAAAGATCGCGCAATGGGTACTTTCGAGGCAGTTATAAACCCAGGTAGCTTATCACCATTCCAGAGCATTTAAGTTCCAGCCACCCATGATTGCCACACGCGCAGTGCTGGGTCATATCCCCCTGGGAGATATGa                                                                                                                                               <  2:1104036/149‑1 (MQ=255)
         tCCGCTTGGACGCAACGTCCGCGCAATGTGTAATTTCTGGGCTGTTGCTCGGACCAGCTTCTAAGCGCCAGGCCAAACCTATTGGTATCCATCCACCCAAGATTGCCACTCGCGCGGTGTTTGGTCATATCTCCCAGGGGGATATGccc                                                                                                                                             <  2:184358/149‑1 (MQ=255)
         cccGCATCGACGAAAAGTCCCCGCAATGGGTAATTTCGTGGTCGTTGTTCGTCCCGGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATGGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAcc                                                                                                                                             <  2:787015/148‑1 (MQ=255)
               tGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTc                                                                                                                                       <  2:461014/149‑1 (MQ=255)
                   cGAAAAGATCGCGCAAGGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCt                                                                                                                                   <  2:500255/149‑1 (MQ=255)
                    gAAAAGAGCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATGGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                  >  2:798269/1‑149 (MQ=255)
                      aaaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGGGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGcc                                                                                                                                >  1:863374/1‑149 (MQ=255)
                       aaGCTCGCGCAGTGGGTAATTTCGGGGCTGTCGTTCGGACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCca                                                                                                                               <  1:888231/149‑1 (MQ=255)
                         gATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCcaca                                                                                                                             >  1:437857/1‑149 (MQ=255)
                         gATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCcaca                                                                                                                             <  2:194552/149‑1 (MQ=255)
                          aTCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGGTATGACCCAGTTCGCCTGCCacac                                                                                                                            >  1:651578/1‑149 (MQ=255)
                            cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACCCCg                                                                                                                          >  1:727528/1‑149 (MQ=255)
                                   tGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAc                                                                                                                   <  2:147743/149‑1 (MQ=255)
                                     ggTAATTCCGGGACTGTTGTCCGTACCAGCGTCTTGGAGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAGCCAGTTCGCCTGCCACACCGTGTGCAccc                                                                                                                 <  2:40671/149‑1 (MQ=255)
                                      gTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACCCGCGCAGTGTTGGGTCGTATCCCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCg                                                                                                                >  2:261811/1‑149 (MQ=255)
                                           ttCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGGCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCAc                                                                                                           >  1:953704/1‑149 (MQ=255)
                                                    gttgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGACACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                                  >  2:782420/1‑149 (MQ=255)
                                                    ggtgttCGTACCAGCGTCTCACGGCCATTCCATAGCAATTGGTTTCCAGGCTCCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGAGATGACCCAGTTCGCCTGCAACACCGTGTGCACCCGTCCACGGCGCACCg                                                                                                  <  2:252593/147‑1 (MQ=255)
                                                          cGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGAGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCGGTCCACGGCGCACCGTGCATc                                                                                            >  1:768986/1‑149 (MQ=255)
                                                          cGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGCTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGTACCCGTCCACGGCGCACCGTTCATc                                                                                            <  1:643276/149‑1 (MQ=255)
                                                                      ttAGCGCCATTCCAGGGCAATTGGTTTCCAGGCCCCCAGGCTTGCCAGACGCGCAGTGTTGTGTTATACCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGGCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                <  2:500905/149‑1 (MQ=255)
                                                                      tgaGCGCCATTCCACCGCAATTGGTTTCCAGGCACCGAGGAGTGCCCCACGCGCAGTGTTGGGTCATACCTCCCAGTGGGACATGACCCAGTCCCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAc                                                                                <  2:467271/147‑1 (MQ=255)
                                                                       caGCGCCAGTCCAGAGCAATTGGTTTCCAGGCCCCCAGCATTGCCACACGCGCAGTGTTGGCTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCCTCCACGGCGCACCGTCCACCCACACTGCGAAcc                                                                               <  2:483879/148‑1 (MQ=255)
                                                                        agggcCATTCCACGGCAATGGGTGTCCAGGCACCCAGGTTTGCCACACGCGCAGCGTTGGGTCATATCTCCCAGGGTGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                                              <  1:626570/146‑1 (MQ=255)
                                                                        aGCGCCATTCCCGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                                              <  1:429489/149‑1 (MQ=255)
                                                                                gccggggCAATTGGTTTCCAGGCCCCCAGGATTGCCACAGGGGCAGTGTTGGGTCATATCTCCCAGGGGCATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCg                                                                      <  2:542536/143‑1 (MQ=255)
                                                                                 ccAGAGCAATTGGTTTCCAGGCACCGAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGt                                                                     >  1:687326/1‑149 (MQ=255)
                                                                                       cgaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgag                                                               <  1:630495/147‑1 (MQ=255)
                                                                                       cAATTGGTTTCCAGGCACTCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCAAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCCCGGCGCACCGTTCATCCACACTCCGAACCCCATCCCCGTACCgag                                                               >  2:741534/1‑149 (MQ=255)
                                                                                        aaTTGGTTTCCAGGCACCCAGGCTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgaga                                                              <  1:1005346/149‑1 (MQ=255)
                                                                                         aTTCGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGAGATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTCCCGAGAt                                                             >  2:104445/1‑149 (MQ=255)
                                                                                          ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATa                                                            <  1:841302/149‑1 (MQ=255)
                                                                                                 ccAGGCACCCAGGATTGCCACACACGCAATGTTGCGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccgcc                                                     >  1:634130/1‑149 (MQ=255)
                                                                                                  cAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGACAGGCCGCCa                                                    >  1:454375/1‑149 (MQ=255)
                                                                                                   aGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAg                                                   >  1:938254/1‑149 (MQ=255)
                                                                                                      cccctatgATCGCCACACGCGCAGTTTGGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCAGCCCCGTACCGAGATAGGCCGCCAGAAc                                                <  2:547896/142‑1 (MQ=255)
                                                                                                           aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGGTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCGGAACCAgtt                                           >  1:865611/1‑149 (MQ=255)
                                                                                                              atggCCACACGCGCAGGTTTGGGTCATCTCTCCCGGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTGCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                        <  2:1116790/146‑1 (MQ=255)
                                                                                                              aTTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgttgtt                                        <  1:105268/149‑1 (MQ=255)
                                                                                                                  ccACACGCTCAGGGTGGGGTCATATCGCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                    <  1:245375/149‑1 (MQ=255)
                                                                                                                  acacgCGCGCAGTGTTGGTTCATACCTCCCAGGGGGATATGTCCCAGTTCGCGTGCCACATCGTGTGCACCCCTCCACGGCGGACCGTTCATGCACACAGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                    <  2:520977/148‑1 (MQ=255)
                                                                                                                   cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCGGTACCGGGATAGGCGGCCAGAACCAGTTGTTGCgga                                   >  1:265038/1‑147 (MQ=255)
                                                                                                                        gcgcAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACAGTGTGCACCCGTCCACGGCACACCGTTCATCCCCACTGCGAACCCCATCCCCGTACCGAGCTAGACACCCAGAACCAGTTGTGGCGTATGGCg                              >  2:881632/1‑149 (MQ=255)
                                                                                                                         cgcAGGTTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGGGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCgg                             <  2:1047258/149‑1 (MQ=255)
                                                                                                                           cAGTGGTGGGTCATATCTCCCAGGGGGATATCGCCCAGTTCGCCTGCCACACCGTGTGCCCCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtt                           <  1:255738/149‑1 (MQ=255)
                                                                                                                             gtgtTGGGTGATCCTCCCCAGGGGGATATGCCCCAGTTCGCCTGCCACCCCGTCTGCACCCGTCCACGGCGCACCGTTCATCCACACTGAGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTTAGGCGGtttt                         <  2:852319/149‑1 (MQ=255)
                                                                                                                              tgtTGGGTCATATCGCCCAGGGGGATATGCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc                        <  1:862001/149‑1 (MQ=255)
                                                                                                                              tgtTGGATCATCTCTCCCAGGGGGGTATGACCCAGTTGGCCGGGCGCACCGTGAGCACCCATACACGGCGCCTCGTTCAGCCACACGGCGAAGCCCATCACCGTACCGACGTAGGCCGCCAGAACCAGTTGTTAGGGAAGGCggatatc                        >  2:183170/1‑144 (MQ=255)
                                                                                                                                   ggTCATATCTCCCAGGGGGATGTGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGGACCAGTTGTTGCGTAAGGCGGTTTTctacta                   >  1:601279/1‑149 (MQ=255)
                                                                                                                                     tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                 >  1:881366/1‑149 (MQ=255)
                                                                                                                                     tCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTCCg                 >  2:845889/1‑149 (MQ=255)
                                                                                                                                     tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                 >  1:725968/1‑149 (MQ=255)
                                                                                                                                     tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                 >  1:725952/1‑149 (MQ=255)
                                                                                                                                         atCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc             >  2:1144193/1‑149 (MQ=255)
                                                                                                                                          tctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa            >  1:818366/1‑149 (MQ=255)
                                                                                                                                          tccccccGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCa            <  1:405841/146‑1 (MQ=255)
                                                                                                                                           ctcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATGGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCGCAg           >  2:470585/1‑149 (MQ=255)
                                                                                                                                             cccAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGga         <  2:938254/149‑1 (MQ=255)
                                                                                                                                               cAGGGGGATTTGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTCCCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgaga       <  2:866922/149‑1 (MQ=255)
                                                                                                                                               cAGGGGGATATGTCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGgaga       <  1:122626/149‑1 (MQ=255)
                                                                                                                                                aGGGGGGTATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag      <  2:437857/149‑1 (MQ=255)
                                                                                                                                                   gggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTg   >  1:1080973/1‑149 (MQ=255)
                                                                                                                                                    ggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCGGCGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGCGGTGc  >  1:561239/1‑149 (MQ=255)
                                                                                                                                                    ggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGc  >  1:992836/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGC  >  NC_000913/4307151‑4307447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑