Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 477,844 C→A 5.8% A258A (GCG→GCT)  ylaB ← putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,8440CA 5.8% 225.3 / 2.4 69A258A (GCG→GCT) ylaBputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCA  >  NC_000913/477699‑477988
                                                                                                                                                 |                                                                                                                                                
gcgcCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCt                                                                                                                                               <  6:38300/149‑1 (MQ=255)
  gcCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGc                                                                                                                                             >  8:131535/1‑149 (MQ=255)
     aGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGGGACTGAATACGGCGCAAAATACGCACCACAAACATCGCTGCcag                                                                                                                                          >  4:131487/1‑149 (MQ=255)
       tGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCcagca                                                                                                                                        <  4:194944/149‑1 (MQ=255)
        gCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCcagcag                                                                                                                                       <  8:37459/149‑1 (MQ=255)
          cTCCGCCCCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCAAGAAAATATCACGAGTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGCATTCGGCGCAAAATACTCAGCTCAAACATCGCTGCCAGCAGGc                                                                                                                                     >  5:40525/1‑149 (MQ=255)
           tCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCATGCAGCCGATGCTTCGGTGACTGAATTCGGCGCAAAATAAGCCGCACAAACATCGCTGCCAGCAGGc                                                                                                                                     >  4:139965/1‑148 (MQ=255)
                aCCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGc                                                                                                                                     <  4:175410/143‑1 (MQ=255)
                aCCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGc                                                                                                                                     >  3:175410/1‑143 (MQ=255)
                aCCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCa                                                                                                                               >  8:177321/1‑149 (MQ=255)
                aCCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCa                                                                                                                               <  7:177321/149‑1 (MQ=255)
                    aCAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCaa                                                                                                                           >  5:211626/1‑149 (MQ=255)
                               aTTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCCCCACCATCGCTGCCAGCAGGCCAATCACCAACCCCAGCGGt                                                                                                                 >  2:171338/1‑148 (MQ=255)
                               aTTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGt                                                                                                                 >  8:237255/1‑148 (MQ=255)
                                    cTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCaa                                                                                                           <  5:176874/149‑1 (MQ=255)
                                      aaGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCaaa                                                                                                          >  4:133624/1‑148 (MQ=255)
                                       aGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAAt                                                                                                         >  7:195091/1‑148 (MQ=255)
                                         gagaCAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATaaa                                                                                                      <  6:212096/149‑1 (MQ=255)
                                         gagaCAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATa                                                                                                        >  8:49205/1‑147 (MQ=255)
                                         gagaCAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATa                                                                                                        >  1:223148/1‑147 (MQ=255)
                                              aaTCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACtt                                                                                                 >  7:180776/1‑149 (MQ=255)
                                              aaTCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACtt                                                                                                 >  8:72491/1‑149 (MQ=255)
                                              aaTCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACtt                                                                                                 >  5:91003/1‑149 (MQ=255)
                                               aTCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTg                                                                                                <  2:223148/149‑1 (MQ=255)
                                                tCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGc                                                                                               <  8:195091/149‑1 (MQ=255)
                                                  ggCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGcc                                                                                              <  1:214511/148‑1 (MQ=255)
                                                    cTGATAGGGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGAt                                                                                           <  6:181305/149‑1 (MQ=255)
                                                      gATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGc                                                                                         >  6:24210/1‑149 (MQ=255)
                                                       aTAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGc                                                                                         >  2:85607/1‑148 (MQ=255)
                                                             cACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCtt                                                                                  >  1:27372/1‑149 (MQ=255)
                                                              aCGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCttt                                                                                 <  1:85607/149‑1 (MQ=255)
                                                                 cAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTg                                                                              <  6:91003/149‑1 (MQ=255)
                                                                           gATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                                                    <  8:180776/149‑1 (MQ=255)
                                                                           gATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTc                                                                    <  5:57517/149‑1 (MQ=255)
                                                                              tttCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                                                 >  1:13372/1‑149 (MQ=255)
                                                                               ttCGATAGTATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                                                 <  6:54508/148‑1 (MQ=255)
                                                                               ttCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                                                 <  3:86997/148‑1 (MQ=255)
                                                                                tCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTg                                                                <  4:117731/148‑1 (MQ=255)
                                                                                       cATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCATTCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAt                                                        <  1:171338/149‑1 (MQ=255)
                                                                                       cATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAt                                                        <  2:177316/149‑1 (MQ=255)
                                                                                           cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtga                                                    <  5:24210/149‑1 (MQ=255)
                                                                                            tGCAGCCGATGATGCGGTGCCTGAATACGGCGCAAAATACGCAGAACAAACATAGCTGCCAGCAGGCCAATCACAAACACGAGCGGTACACAAATAAAGACTTGACGATGCCATCCTTTCTGCAACATTTTCGtt                                                                 >  8:201661/1‑135 (MQ=255)
                                                                                            tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATAGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                                                 <  7:205182/135‑1 (MQ=255)
                                                                                            tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATAGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                                                 >  8:205182/1‑135 (MQ=255)
                                                                                            tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATAGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGtt                                                                 <  7:201661/135‑1 (MQ=255)
                                                                                             gCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                                                   >  1:98032/1‑148 (MQ=255)
                                                                                                    atgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                                           >  7:104256/1‑149 (MQ=255)
                                                                                                       atgCGGTGACTGAATACGGCGCAAAATACTCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAata                                         <  3:131487/148‑1 (MQ=255)
                                                                                                       atgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat                                        <  2:98032/149‑1 (MQ=255)
                                                                                                        tgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat                                        >  4:126631/1‑148 (MQ=255)
                                                                                                          cGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTc                                      <  7:72491/148‑1 (MQ=255)
                                                                                                               aCTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtctc                                 <  8:179814/148‑1 (MQ=255)
                                                                                                               aCTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtctc                                 <  3:133624/148‑1 (MQ=255)
                                                                                                                 tGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGAATTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCgg                              >  5:97447/1‑149 (MQ=255)
                                                                                                                   aaTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGt                             >  2:50776/1‑148 (MQ=255)
                                                                                                                         gcgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAggg                        >  6:180726/1‑147 (MQ=255)
                                                                                                                          cgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAg                     <  7:131535/149‑1 (MQ=255)
                                                                                                                           gcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAgg                    <  8:104256/149‑1 (MQ=255)
                                                                                                                               aaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGatat                >  2:172113/1‑149 (MQ=255)
                                                                                                                                 tACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCg              <  7:237255/149‑1 (MQ=255)
                                                                                                                                     cagcaCAACCCTCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACATGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTaaa          <  2:190133/149‑1 (MQ=255)
                                                                                                                                         acaAACCTCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCACATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCTTTGCAGCCCAtgt                                                      <  3:158451/101‑1 (MQ=255)
                                                                                                                                         acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTaa                           <  5:223570/128‑1 (MQ=255)
                                                                                                                                         acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTaa                           >  6:223570/1‑128 (MQ=255)
                                                                                                                                          caAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCg              >  2:69670/1‑140 (MQ=255)
                                                                                                                                          caAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCg              <  1:69670/140‑1 (MQ=255)
                                                                                                                                           aaaCATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTAATCTCCGGTAAGTGAAGGATATCGTAAATGATTc    >  6:128752/1‑149 (MQ=255)
                                                                                                                                            aaCATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTc    >  4:225007/1‑148 (MQ=255)
                                                                                                                                             aCATCGCTGCCAGCAGGCCAATCACCACCCCGAGCGGTACCCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCa  >  5:238677/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGCCAGTGCCTCAGCACCGACAATTTTGCCATTGGCTAAGGAGACAATCGGCTGATAGTGCACGCAAATATCACGATTTTCGATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCA  >  NC_000913/477699‑477988

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