Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 447,950 C→T 5.8% L177L (CTG→TTG)  PP_0369 → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947447,9500CT 5.8% 203.9 / 2.2 69L177L (CTG→TTG) PP_0369diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base T (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

GTGCGGTGCTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAACATCATGTGCAGCCAG  >  NC_002947/447804‑448096
                                                                                                                                                   |                                                                                                                                                  
gtgcggtgcTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTg                                                                                                                                                  <  6:239435/149‑1 (MQ=255)
gtgcggtgcTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTg                                                                                                                                                  <  2:104023/149‑1 (MQ=255)
 tgcggtgcTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTgg                                                                                                                                                 <  1:186764/149‑1 (MQ=255)
 tgcggtgcTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTgg                                                                                                                                                 <  2:296038/149‑1 (MQ=255)
     gtgcTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGgtgt                                                                                                                                             >  6:18591/1‑149 (MQ=255)
        cTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGgtgtgt                                                                                                                                           >  6:340009/1‑148 (MQ=255)
        cTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGtt                                                                                                                                          <  5:340009/149‑1 (MQ=255)
           tttGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCAt                                                                                                                                       >  7:216886/1‑149 (MQ=255)
             tGCCGTGGGGAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCAtt                                                                                                                                      <  3:232543/148‑1 (MQ=255)
             tGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCAtt                                                                                                                                      <  6:164891/148‑1 (MQ=255)
               ccGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGgtgtgt                                                                                                                                           <  3:210483/141‑1 (MQ=255)
               ccGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGgtgtgt                                                                                                                                           >  4:210483/1‑141 (MQ=255)
                     gCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCTGGGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCTTGGGTGTGTTCATTCTGCtggt                                                                                                                              >  1:296080/1‑149 (MQ=255)
                     gCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCTGGGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCTTGGGTGTGTTCATTCTGCtggt                                                                                                                              >  1:291129/1‑149 (MQ=255)
                      cAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTCTTCCTGGGTGTGTTCATTCTGCtggtgg                                                                                                                            >  2:192491/1‑149 (MQ=255)
                        gtgtGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGgtggtg                                                                                                                          <  2:371115/149‑1 (MQ=255)
                                gTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCt                                                                                                                  <  2:319805/149‑1 (MQ=255)
                                     ccTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTg                                                                                                             >  5:42550/1‑149 (MQ=255)
                                     ccTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTg                                                                                                             >  3:98443/1‑149 (MQ=255)
                                       tCCAGCAAAGCTACCCCGAGGCGTTGATGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCAATGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCa                                                                                                           >  6:145617/1‑149 (MQ=255)
                                         cAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCAcc                                                                                                         >  7:261646/1‑149 (MQ=255)
                                            cAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCg                                                                                                      >  2:68968/1‑149 (MQ=255)
                                             aaaGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTACATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTAcc                                                                                                                   >  3:102872/1‑135 (MQ=255)
                                             aaaGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTACATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTAcc                                                                                                                   <  4:102872/135‑1 (MQ=255)
                                               aGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGcc                                                                                                   >  2:331835/1‑149 (MQ=255)
                                                   aCCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTg                                                                                               <  1:192491/149‑1 (MQ=255)
                                                    ccccGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCTTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCt                                                                                                <  2:296080/147‑1 (MQ=255)
                                                    ccccGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCTTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCt                                                                                                <  2:291129/147‑1 (MQ=255)
                                                    ccccGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCt                                                                                                <  8:252572/147‑1 (MQ=255)
                                                    ccccGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGc                                                                                              >  3:241087/1‑149 (MQ=255)
                                                    ccccGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGc                                                                                              <  7:222691/149‑1 (MQ=255)
                                                       cGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCAcg                                                                                           >  2:304374/1‑149 (MQ=255)
                                                          ggCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCAcgccg                                                                                        >  5:209242/1‑149 (MQ=255)
                                                            cGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCcgcg                                                                                      <  6:339844/149‑1 (MQ=255)
                                                                   gACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGa                                                                               >  6:337371/1‑149 (MQ=255)
                                                                    aCCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGa                                                                               <  3:34461/148‑1 (MQ=255)
                                                                      cT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCg                                                                            <  4:241087/149‑1 (MQ=255)
                                                                           ttGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCct                                                                        <  4:98443/149‑1 (MQ=255)
                                                                             ggCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCctgc                                                                      >  2:141391/1‑149 (MQ=255)
                                                                                ggAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCctgct                                                                     <  8:250684/147‑1 (MQ=255)
                                                                                ggAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCctgct                                                                     >  7:250684/1‑147 (MQ=255)
                                                                                    aTCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCAt                                                               >  5:252170/1‑149 (MQ=255)
                                                                                    aTCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTCTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCAt                                                               >  5:257690/1‑149 (MQ=255)
                                                                                       cGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCAtttt                                                            >  3:340100/1‑149 (MQ=255)
                                                                                         cTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCt                                                          <  6:42550/149‑1 (MQ=255)
                                                                                            gCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGa                                                       <  8:220925/149‑1 (MQ=255)
                                                                                            gCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGa                                                       <  6:209242/149‑1 (MQ=255)
                                                                                             ccGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGAt                                                      <  5:145617/149‑1 (MQ=255)
                                                                                              cGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTg                                                         >  2:353610/1‑145 (MQ=255)
                                                                                              cGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTg                                                         <  1:353610/145‑1 (MQ=255)
                                                                                              cGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGAtg                                                     >  4:375066/1‑149 (MQ=255)
                                                                                                      cATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCg                                               >  5:163943/1‑147 (MQ=255)
                                                                                                        tGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGa                                           <  1:68968/149‑1 (MQ=255)
                                                                                                            cACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCtt                                       >  5:346748/1‑149 (MQ=255)
                                                                                                                 gATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATc                                   <  1:108384/148‑1 (MQ=255)
                                                                                                                        cTGATCCAGATTGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTTCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGtt                           <  6:257690/149‑1 (MQ=255)
                                                                                                                        cTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGtt                           <  6:252170/149‑1 (MQ=255)
                                                                                                                         tGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCtc                                                         >  2:156921/1‑117 (MQ=255)
                                                                                                                         tGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTg                                                         <  1:156921/118‑1 (MQ=255)
                                                                                                                          gATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTa                         >  4:320525/1‑149 (MQ=255)
                                                                                                                            tCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTAcc                       <  1:304374/149‑1 (MQ=255)
                                                                                                                              cAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCtt                     >  8:313247/1‑149 (MQ=255)
                                                                                                                                gATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTg                    >  7:150097/1‑148 (MQ=255)
                                                                                                                                 aTGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGaa                  <  6:163943/149‑1 (MQ=255)
                                                                                                                                 aTGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGaa                  <  8:150097/149‑1 (MQ=255)
                                                                                                                                      cTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCTGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAAcatca             >  3:207027/1‑149 (MQ=255)
                                                                                                                                       tACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAAcatca             <  1:331835/148‑1 (MQ=255)
                                                                                                                                        aCCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAACATCAtg           >  6:155131/1‑149 (MQ=255)
                                                                                                                                         cccGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAACATCAtgt          <  3:375066/149‑1 (MQ=255)
                                                                                                                                                  ccTGGGTGTGTTCATTCTGCTGGTGCTGCAGTACCTGCGGGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGTTGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAACATCATGTGCAGCCAg  <  6:283047/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGCGGTGCTGTTTGCCGTGGGCAGTGTGGCGGTGGCCCTCCAGCAAAGCTACCCCGAGGCGTTGCTGACCT‑GGTTGGCGGAGATCCGCTGGCCGGCTTTGCATGGCCACTGGATGAGCCTGATCCAGATGGTCTACCCGCTGTTCCTGGGTGTGTTCATTCTGCTGGTGGTGCAGTACCTGCGTGCACCCCGGCCGCTGCACGCCGCGCAGTTGATCGGCCTGCTGGGCATTTTCTGGATGCTGCCGCAGACCTTCATCCTGCCGTTTACCTTGAACATCATGTGCAGCCAG  >  NC_002947/447804‑448096

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: ‑