Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,535,691 G→A 6.6% S56S (AGC→AGT)  yddE ← PhzC‑PhzF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,535,6910GA 6.6% 281.8 / 5.9 91S56S (AGC→AGT) yddEPhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base G (46/38);  new base A (3/3);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCAT  >  NC_000913/1535548‑1535833
                                                                                                                                               |                                                                                                                                              
cgATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCACTgtc                                                                                                                                           >  4:414877/1‑149 (MQ=255)
cgATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCACTgtc                                                                                                                                           >  4:421562/1‑149 (MQ=255)
     ttttCCTGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGtcgtcg                                                                                                                                       <  7:361398/148‑1 (MQ=255)
           tGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCtgtg                                                                                                                                  >  7:127623/1‑147 (MQ=255)
           tGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAAACTAACACCTTCGCACGTACATAGAGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACATCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTgca                                                                                                                                >  8:192058/1‑149 (MQ=255)
              cAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTgcag                                                                                                                               >  1:268297/1‑147 (MQ=255)
                   aTGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGaaaa                                                                                                                        <  7:392568/149‑1 (MQ=255)
                             gATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTcc                                                                                                               >  3:415090/1‑148 (MQ=255)
                                  tGCAATTTCCTAAACCTAACACCTTCTCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCa                                                                                                                     >  8:449553/1‑137 (MQ=255)
                                  tGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTgcagca                                                                                                                             >  1:309542/1‑129 (MQ=255)
                                  tGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTgcagca                                                                                                                             <  2:309542/129‑1 (MQ=255)
                                  tGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCa                                                                                                                     <  7:449553/137‑1 (MQ=255)
                                  tGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATgg                                                                                                         >  8:141663/1‑149 (MQ=255)
                                    cAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGcc                                                                                                       >  7:92193/1‑149 (MQ=255)
                                     aaTTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGcc                                                                                                       >  4:272595/1‑148 (MQ=255)
                                     aaTTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCACTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGcc                                                                                                       >  5:345175/1‑148 (MQ=255)
                                     aaTTTCCTAAACCTAACACCTTCGCACGTACATAGTGAGCAGCTACTGTCTCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGcc                                                                                                       >  5:81370/1‑148 (MQ=255)
                                      aTTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCt                                                                                                      >  2:348931/1‑148 (MQ=255)
                                            tAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGaaaa                                                                                                                        <  7:197430/124‑1 (MQ=255)
                                            tAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTcgc                                                                                               >  6:278972/1‑149 (MQ=255)
                                            tAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTGGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTcgc                                                                                               >  1:319123/1‑149 (MQ=255)
                                             aaaCCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTcgcg                                                                                              >  5:394605/1‑149 (MQ=255)
                                              aaCCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTcgcg                                                                                              <  6:394605/148‑1 (MQ=255)
                                              aaCCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCgg                                                                                             <  8:92193/149‑1 (MQ=255)
                                                     cacCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGaa                                                                                                            <  3:314881/127‑1 (MQ=255)
                                                     cacCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGaa                                                                                                            >  4:314881/1‑127 (MQ=255)
                                                       ccTTCGCACGTACATAGGGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGc                                                                                    <  1:144197/149‑1 (MQ=255)
                                                        cTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCt                                                                                   >  8:361702/1‑149 (MQ=255)
                                                           cGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGc                                                                                 <  6:199904/148‑1 (MQ=255)
                                                              aCGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCAttt                                                                             <  7:141663/149‑1 (MQ=255)
                                                              aCGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTATCCGAATGGCCTAACTCGCGGGCGATAAGCTGCAttt                                                                             <  8:127623/149‑1 (MQ=255)
                                                                 tACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTgcg                                                                          <  5:86645/149‑1 (MQ=255)
                                                                     tAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTc                                                                      >  6:371814/1‑149 (MQ=255)
                                                                     tAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTc                                                                      >  5:403531/1‑149 (MQ=255)
                                                                     tAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTc                                                                      >  4:95171/1‑149 (MQ=255)
                                                                          cgcAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTaa                                                                 <  4:415090/149‑1 (MQ=255)
                                                                           gcAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGaaaa                                                                                                                        >  5:316481/1‑93 (MQ=255)
                                                                           gcAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGaaaa                                                                                                                        <  6:316481/93‑1 (MQ=255)
                                                                           gcAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCAt                                                                               <  7:154259/134‑1 (MQ=255)
                                                                           gcAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCAt                                                                               >  8:154259/1‑134 (MQ=255)
                                                                            cgatctgtGTCGCGTGTCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCATCAGAAATGCAGTTTCCGAATTGACTAACTCGCGGGCGATAAGCTGGATTTGCGCTTCGCTAAGa                                                               <  3:27766/142‑1 (MQ=255)
                                                                                  cTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCACTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTAt                                                           <  3:421562/147‑1 (MQ=255)
                                                                                  cTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCACTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTAt                                                           <  3:414877/147‑1 (MQ=255)
                                                                                   tGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGa                                                               >  2:38532/1‑142 (MQ=255)
                                                                                   tGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGa                                                               <  1:38532/142‑1 (MQ=255)
                                                                                   tGTCGCGTCACCGCAAATCGGCACTTCCACCGTTGGCGTAAAGTAGCGGATGCGCCCCTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCCCGGGCCATAAGCTGCATTTGCGCTTCGCTAAGa                                                               >  4:27766/1‑142 (MQ=255)
                                                                                     tCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATcggcg                                                      >  8:271000/1‑149 (MQ=255)
                                                                                     tCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATcggcg                                                      <  5:278972/149‑1 (MQ=255)
                                                                                     tCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATcggcg                                                      >  6:202633/1‑149 (MQ=255)
                                                                                       gcgTGACCGCAAATCGGCACTTCTACCGTTGGCTTAAAGTAGCGGATGCGCACATCGCTGTCGTCTCTGTGAAGCAGAACAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCggg                                                    <  6:81370/149‑1 (MQ=255)
                                                                                       gcgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATcggcgg                                                     <  7:361702/148‑1 (MQ=255)
                                                                                       gcgTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATcggcgg                                                     <  2:49266/148‑1 (MQ=255)
                                                                                           gACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATc                                                          >  8:329676/1‑139 (MQ=255)
                                                                                           gACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATc                                                          <  7:329676/139‑1 (MQ=255)
                                                                                            aCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGaaaa                                               <  5:371814/149‑1 (MQ=255)
                                                                                            aCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGaaaa                                               <  2:136247/149‑1 (MQ=255)
                                                                                              cGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaac                                              >  7:305109/1‑148 (MQ=255)
                                                                                               gCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGGTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaacaa                                            >  8:2820/1‑149 (MQ=255)
                                                                                               gCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaacaa                                            >  3:137793/1‑149 (MQ=255)
                                                                                                 aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAacaaca                                          >  8:147439/1‑149 (MQ=255)
                                                                                                 aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAacaaca                                          >  5:295479/1‑149 (MQ=255)
                                                                                                   aTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACAcc                                        >  4:185307/1‑149 (MQ=255)
                                                                                                           ttCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAAt                                >  2:166162/1‑149 (MQ=255)
                                                                                                           ttCAACAGTTGTCGTAAAGTAGAGTATGCGCACATCTCTGTAGTAGCTGTGCAGCAGAAAAGCAGTTTCCTAATGGACTACCTCGCGGTCGATAAGCTGCATTTGCGCTTCGCTAAGAGTATCGGCTGGGAAAACACAACCGGCAGAAt                                >  8:295129/1‑149 (MQ=255)
                                                                                                            tCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAAtt                               >  1:105901/1‑149 (MQ=255)
                                                                                                                      gCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGaaa                      <  3:185307/148‑1 (MQ=255)
                                                                                                                      gCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAAc                     >  8:107786/1‑149 (MQ=255)
                                                                                                                         tAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCt                                                                                   <  6:149279/84‑1 (MQ=255)
                                                                                                                           aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGtt                 <  4:137793/148‑1 (MQ=255)
                                                                                                                           aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGtt                 <  1:166162/148‑1 (MQ=255)
                                                                                                                           aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtg                <  7:192058/149‑1 (MQ=255)
                                                                                                                           aaGTAGCGGATGCGCACATCACTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGtt                 <  6:345175/148‑1 (MQ=255)
                                                                                                                             gTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgtg              <  5:202633/149‑1 (MQ=255)
                                                                                                                             gTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt               >  7:337631/1‑148 (MQ=255)
                                                                                                                              tAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt               >  7:153225/1‑147 (MQ=255)
                                                                                                                              tAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt               >  8:170614/1‑147 (MQ=255)
                                                                                                                              tAGCGGATGCGCACATCGCTGGCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGa             <  2:149717/149‑1 (MQ=255)
                                                                                                                                gCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCgg                                      <  6:189481/122‑1 (MQ=255)
                                                                                                                                gCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCgg                                      >  5:189481/1‑122 (MQ=255)
                                                                                                                                gCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGtt                 >  7:304534/1‑143 (MQ=255)
                                                                                                                                gCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGtt                 <  8:304534/143‑1 (MQ=255)
                                                                                                                                     tGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTaa        >  1:245983/1‑147 (MQ=255)
                                                                                                                                        gcACATCGCTGTCGTCGCTGTGCAGCGGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGa             >  6:166242/1‑139 (MQ=255)
                                                                                                                                        gcACATCGCTGTCGTCGCTGTGCAGCGGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGa             <  5:166242/139‑1 (MQ=255)
                                                                                                                                        gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt               <  8:434459/137‑1 (MQ=255)
                                                                                                                                        gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTtgt               >  7:434459/1‑137 (MQ=255)
                                                                                                                                        gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCa   >  1:257673/1‑149 (MQ=255)
                                                                                                                                        gcACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCa   >  7:385778/1‑149 (MQ=255)
                                                                                                                                         cacaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCAt  <  8:337631/149‑1 (MQ=255)
                                                                                                                                          acaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTg                                                                            <  4:91513/74‑1 (MQ=255)
                                                                                                                                          acaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTg                                                                            >  3:91513/1‑74 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CGATGTTTTCCTGCCAGCGATGTTTGCCAGATCGTGCAATTTCCTAAACCTAACACCTTCGCACGTACATAGTGCGCAGCTACTGTCGCGTGACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGTGAGGTAAAAGCAT  >  NC_000913/1535548‑1535833

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