Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 404,513 C→A 5.7% L271I (CTT→ATT)  yaiC → diguanylate cyclase, cellulose regualtor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913404,5130CA 5.7% 258.7 / 3.7 79L271I (CTT→ATT) yaiCdiguanylate cyclase, cellulose regualtor
Reads supporting (aligned to +/- strand):  ref base C (28/46);  new base A (3/2);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTAC  >  NC_000913/404369‑404660
                                                                                                                                                |                                                                                                                                                   
cATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTAc                                                                                                                                                 >  1:255600/1‑149 (MQ=255)
 aTTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTAcc                                                                                                                                                >  5:308412/1‑149 (MQ=255)
     ggAAACTATGTTAAGCAATGAATTTCAAAACTGTCCGCGGCATAATCGCGATGCAACTTTACAGATAATAGATAACGACCATTTCAAGAGCATCAACGATAAATGGGGCCAAGATGTGGGCGATGAAGCGATTGTGCCGCTTACCCGa                                                                                                                                             >  6:126794/1‑148 (MQ=255)
        aaCTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCTATGCCACGTTAATGATTATCGATATCGCCCATTTCAAGAGAATCAACGATCCCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCCCTTACCCGCCAGt                                                                                                                                         >  1:210477/1‑149 (MQ=255)
        aaCTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGt                                                                                                                                         >  5:249591/1‑149 (MQ=255)
        aaCTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGt                                                                                                                                         >  6:187689/1‑149 (MQ=255)
          cTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTa                                                                                                                                       >  5:74055/1‑149 (MQ=255)
             tGTTACGAAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTAc                                                                                                                                      <  2:243406/147‑1 (MQ=255)
               tatcGCATTGAATTTGATATCTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTCGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACaaa                                                                                                                                   <  5:126794/145‑1 (MQ=255)
                  cGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTa                                                                                                                                <  5:243336/148‑1 (MQ=255)
                  cGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTAc                                                                                                                               <  6:212581/149‑1 (MQ=255)
                   gCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTAcc                                                                                                                              <  7:292204/149‑1 (MQ=255)
                    cAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTAccc                                                                                                                             <  2:210477/149‑1 (MQ=255)
                     aaTGAATTTGATAACTGTCGGCGGCATAATCGCGATGCCACGTTACTGATTTTCGCGATCGACCATTTCAAGAGCATCAACGATCCCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCt                                                                                                                            >  1:134568/1‑149 (MQ=255)
                     aaTGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGTATAAACGATACCTGGGGCCATGATTTGGGCGATGAAGCGATTGTGGCGCTTACAATAAAGTTACAAATTAACCt                                                                                                                            >  2:54905/1‑149 (MQ=255)
                     aaTGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCt                                                                                                                            <  1:240022/149‑1 (MQ=255)
                     aaTGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCt                                                                                                                            >  1:287738/1‑149 (MQ=255)
                     aaTGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGATTACCCGACAGTTACAAATTACCCt                                                                                                                            >  2:240341/1‑149 (MQ=255)
                     aaTGAAATTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGATTACCCGACAGTTACAAATTACCCt                                                                                                                            >  2:235752/1‑149 (MQ=255)
                       tGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACa                                                                                                                                           <  8:186464/132‑1 (MQ=255)
                              gATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAg                                                                                                                   >  1:180463/1‑149 (MQ=255)
                                 aaCTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGAATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGa                                                                                                                >  5:293236/1‑149 (MQ=255)
                                  aCTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACTATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGAGAt                                                                                                               <  2:254391/149‑1 (MQ=255)
                                  aCTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAt                                                                                                               <  6:249591/149‑1 (MQ=255)
                                  aCTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAt                                                                                                               <  3:303669/149‑1 (MQ=255)
                                      tCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAtgtg                                                                                                            <  7:156721/148‑1 (MQ=255)
                                      tCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGAtgt                                                                                                             <  3:297141/147‑1 (MQ=255)
                                        ggcggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGATTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGAt                                                                                                          <  1:240341/148‑1 (MQ=255)
                                        ggcggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGATTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGAt                                                                                                          <  1:235752/148‑1 (MQ=255)
                                         gcggcATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTg                                                                                                        >  1:238833/1‑149 (MQ=255)
                                             cATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCg                                                                                                    <  6:288116/149‑1 (MQ=255)
                                                 aTCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGttt                                                                                                <  8:125685/149‑1 (MQ=255)
                                                    gcgATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggc                                                                                             >  7:107385/1‑149 (MQ=255)
                                                       aTGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTggcggc                                                                                          <  2:287738/149‑1 (MQ=255)
                                                         gCAACGTTACTGATTATCGATATCGACCATTTAAAGAGCATCAAAGATACCTGGGGCCATGAAGTGGGCGATGAAGCGATTGTGGCGATTACACGACAGTTACAAATTACCCTGCGCGGTAGCGAGGTGATTGGTCGGTTTtgcggcta                                                                                        >  6:72273/1‑147 (MQ=255)
                                                           aaCGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGAt                                                                                       >  8:267677/1‑148 (MQ=255)
                                                            aCGTTAGTGATTATCGATATCGACCCTTTCAAGAGTATCAACGATACCTGGGGCAATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGa                                                                                     <  1:54905/149‑1 (MQ=255)
                                                            aCGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGa                                                                                     <  4:290437/149‑1 (MQ=255)
                                                            aCGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGa                                                                                     <  7:222682/149‑1 (MQ=255)
                                                               ttACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGtt                                                                                  <  8:79027/149‑1 (MQ=255)
                                                                  cTGATTATCGATATCGACCATTTCAAGAGAATCAACGTTAGCTGGGGCCTTGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTg                                                                                <  8:107385/148‑1 (MQ=255)
                                                                   tGATTATCGATATCTACCATTTCTAGAGCATCAACGATTGCTGGGGCCATGATGTGGGCGATGAAGCGATTGTTGCGCTTACCCGACAGTTACAAATTTCCCTTCGCGGTAGCGTTGTGATTGGTCGGTTTGGCGGCGATGAGTTTTCa                                                                              <  6:293236/149‑1 (MQ=255)
                                                                     agtatcgatatcgaCCATGTCAAGAGCAGCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACACGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                                                                              <  8:304734/145‑1 (MQ=255)
                                                                     aTTATCGATATCGACCATTTCAAGAGCTTCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGGGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                                                                              <  6:19854/147‑1 (MQ=255)
                                                                     aTTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                                                                              <  2:238833/147‑1 (MQ=255)
                                                                     aTTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                                                                              <  4:164429/147‑1 (MQ=255)
                                                                         tcgatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAAt                                                                         >  6:325802/1‑148 (MQ=255)
                                                                         tcgatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATc                                                                        <  2:292631/149‑1 (MQ=255)
                                                                         tcgatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGAGGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAAt                                                                         <  6:303087/148‑1 (MQ=255)
                                                                           gatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCa                                                                       <  5:187689/148‑1 (MQ=255)
                                                                           gatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCAt                                                                      <  7:264696/149‑1 (MQ=255)
                                                                           gatatcgaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCAt                                                                      <  6:74055/149‑1 (MQ=255)
                                                                            atatcgcCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATg                                                                     >  4:117118/1‑149 (MQ=255)
                                                                                 gaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCg                                                                                                    >  1:112359/1‑113 (MQ=255)
                                                                                 gaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCg                                                                                                    <  2:112359/113‑1 (MQ=255)
                                                                                 gaCCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCgg                                                                >  3:143175/1‑149 (MQ=255)
                                                                                  aacATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGGCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGt                                                               <  7:81668/147‑1 (MQ=255)
                                                                                    cATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTAc                                                             >  3:82469/1‑149 (MQ=255)
                                                                                     aTTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                                                                              <  4:8037/131‑1 (MQ=255)
                                                                                     aTTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCa                                                                              >  3:8037/1‑131 (MQ=255)
                                                                                         cAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCa                                                         <  2:255600/148‑1 (MQ=255)
                                                                                         cAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCa                                                         <  2:134568/148‑1 (MQ=255)
                                                                                                aTCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGagag                                                 >  7:301977/1‑149 (MQ=255)
                                                                                                 tCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAgc                                                >  6:115725/1‑149 (MQ=255)
                                                                                                    aCGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTa                                                                                                                    <  3:253459/78‑1 (MQ=255)
                                                                                                      gATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCAt                                                                      >  6:33544/1‑122 (MQ=255)
                                                                                                      gATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCAt                                                                      <  5:33544/122‑1 (MQ=255)
                                                                                                      gATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCAt                                           <  5:243840/149‑1 (MQ=255)
                                                                                                         aCCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTAc                                        >  4:73314/1‑149 (MQ=255)
                                                                                                          ccTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTAc                                        <  2:180463/148‑1 (MQ=255)
                                                                                                             ggggCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGcc                                    <  2:156922/149‑1 (MQ=255)
                                                                                                               ggCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCAt                                  <  3:117118/149‑1 (MQ=255)
                                                                                                                gCCATTATGTGGGCGATGAAGCGTTTGTGGCGCTTACCCGACATTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGGCGGGTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATg                                 <  5:72273/149‑1 (MQ=255)
                                                                                                                   atgatgTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTa                              >  3:297010/1‑149 (MQ=255)
                                                                                                                      atgTGGGCGAGGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCCGTAGCGATGTGATTGGTCGGATTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACgg                           <  6:308412/149‑1 (MQ=255)
                                                                                                                           ggCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCa                      >  4:77439/1‑149 (MQ=255)
                                                                                                                                    gCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCt             <  4:82469/149‑1 (MQ=255)
                                                                                                                                      gATTGTGGCTCTTACCCGACAGTTACAAATTCCCCTGAGCGGTAGCGATGTTATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTActgttga                       <  3:30382/137‑7 (MQ=255)
                                                                                                                                      gATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGAGGCGATGAGTTTTCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGc                       >  4:30382/1‑137 (MQ=255)
                                                                                                                                             gcgcTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACAt     >  8:67199/1‑148 (MQ=255)
                                                                                                                                               gcTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTAc  <  3:73314/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CATTGGGAAACTATGTTACGCAATGAATTTGATAACTGTCGGCGGCATAATCGCGATGCAACGTTACTGATTATCGATATCGACCATTTCAAGAGCATCAACGATACCTGGGGCCATGATGTGGGCGATGAAGCGATTGTGGCGCTTACCCGACAGTTACAAATTACCCTGCGCGGTAGCGATGTGATTGGTCGGTTTGGCGGCGATGAGTTTGCAGTAATCATGTCCGGTACGCCAGCTGAGAGCGCCATTACCGCCATGTTACGGGTGCATGAAGGGCTAAATACATTAC  >  NC_000913/404369‑404660

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