Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,437,180 C→T 8.0% R114C (CGC→TGC)  PP_3932 → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,437,1800CT 8.0% 179.7 / 4.5 62R114C (CGC→TGC) PP_3932diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (2/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

GGGCTTTAGCCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTT  >  NC_002947/4437032‑4437328
                                                                                                                                                    |                                                                                                                                                    
gggCTTTAGCCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGgc                                                                                                                                                      >  4:4291/1‑149 (MQ=255)
gggCTTTAGCCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGgc                                                                                                                                                      <  8:362398/149‑1 (MQ=255)
  gCTTTAGCCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGgcg                                                                                                                                                     >  3:226841/1‑148 (MQ=255)
     ttAGCCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATg                                                                                                                                                 >  4:28835/1‑149 (MQ=255)
          cTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCGCCGGTGCATgcgtgg                                                                                                                                            <  8:61468/149‑4 (MQ=255)
          cTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAg                                                                                                                                            >  7:322587/1‑149 (MQ=255)
             gtgtGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGc                                                                                                                                           <  4:226841/147‑1 (MQ=255)
             gtgtGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGc                                                                                                                                           <  6:375324/147‑1 (MQ=255)
                tggctggcCATGCTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACATGTGCAATCTTGGCGTGGTGGTCACCCTCTAGCAGCTGCTGTCTCGCTTGCCACGCATGAACAGACGCATGCGCAGCCCgaa                                                                                                                                      >  8:116141/1‑147 (MQ=255)
                    tggcCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGc                                                                                                                                  >  4:365936/1‑149 (MQ=255)
                         aTGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCCTGCGCAGCCCGCACGGCGTACt                                                                                                                             >  5:764/1‑149 (MQ=255)
                                 ggCCGACCTGCTGTGTGGCAGCGCCGGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCgtgt                                                                                                                     <  8:2683/149‑1 (MQ=255)
                                   ccGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTc                                                                                                                   <  8:322587/149‑1 (MQ=255)
                                    cGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTcg                                                                                                                  >  6:19622/1‑149 (MQ=255)
                                       ccTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTcgccg                                                                                                               >  1:212888/1‑149 (MQ=255)
                                        ctgctgTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCGCCGGTGCATgcgt                                                                                                                                              >  6:173738/1‑116 (MQ=255)
                                        ctgctgTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCGCCGGTGCATgcgt                                                                                                                                              <  5:173738/117‑2 (MQ=255)
                                         tgctgTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCgccgcc                                                                                                             >  2:210272/1‑149 (MQ=255)
                                          gctgTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCt                                                                                                            >  3:106028/1‑149 (MQ=255)
                                                  gCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGcgccgcc                                                                                                    <  7:281936/149‑1 (MQ=255)
                                                          tGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAt                                                                                            >  3:78863/1‑149 (MQ=255)
                                                             gTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggt                                                                                         <  2:212888/149‑1 (MQ=255)
                                                             gTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTCTTCGCCGCCTGCGCCGCCGCTGACAtgt                                                                                          >  4:136340/1‑147 (MQ=255)
                                                             gTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGGCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggt                                                                                         <  6:359118/149‑1 (MQ=255)
                                                                 cGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAg                                                                                                                          <  7:86649/112‑1 (MQ=255)
                                                                 cGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAg                                                                                                                          >  8:86649/1‑112 (MQ=255)
                                                                     tctGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTggcggc                                                                                   >  6:97872/1‑147 (MQ=255)
                                                                       tgtgGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGa                                                                                 <  6:6155/147‑1 (MQ=255)
                                                                       tgtgGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGa                                                                                 <  5:19622/147‑1 (MQ=255)
                                                                        gtgGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCa                                                                              >  1:114160/1‑149 (MQ=255)
                                                                         tgGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCCACCATGGCAgcggcggc                                                                                <  3:28835/146‑3 (MQ=255)
                                                                            ttAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCGCCGGTGCAt                                                                                                                                                  >  1:38455/1‑77 (MQ=255)
                                                                            ttAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCGCCGGTGCAt                                                                                                                                                  <  2:38455/77‑1 (MQ=255)
                                                                                aCCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCg                                                                      <  2:114160/149‑1 (MQ=255)
                                                                                aCCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCg                                                                      >  3:118212/1‑149 (MQ=255)
                                                                                  cTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGt                                                                    <  7:116141/149‑1 (MQ=255)
                                                                                     tgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                 >  3:248717/1‑149 (MQ=255)
                                                                                     tgCAATCTTGGCGTGGTGGTCACACTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                 >  5:162343/1‑149 (MQ=255)
                                                                                       cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATgg                                                               <  4:118212/149‑1 (MQ=255)
                                                                                        aaTCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGc                                                              <  5:129935/149‑1 (MQ=255)
                                                                                          tCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCg                                                            >  6:111448/1‑149 (MQ=255)
                                                                                               gCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCAcc                                                        <  3:4291/148‑1 (MQ=255)
                                                                                                            ccTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                 >  3:75233/1‑126 (MQ=255)
                                                                                                            ccTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                 <  4:75233/126‑1 (MQ=255)
                                                                                                                tGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGcgccgccgc                                                                                                  >  4:256747/1‑89 (MQ=255)
                                                                                                                tGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGcgccgccgc                                                                                                  <  3:256747/89‑1 (MQ=255)
                                                                                                                tGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGtt                                      >  7:127604/1‑149 (MQ=255)
                                                                                                                  gCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAg                                                                                >  8:227777/1‑105 (MQ=255)
                                                                                                                  gCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAg                                                                                <  7:227777/105‑1 (MQ=255)
                                                                                                                  gCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTgc                                    >  3:57802/1‑149 (MQ=255)
                                                                                                                      ctgctgTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTgcgcgc                                <  5:97872/149‑1 (MQ=255)
                                                                                                                      ctgctgTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGGGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTgcgcgc                                <  3:365936/149‑1 (MQ=255)
                                                                                                                          tgTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCAcc                            >  3:85695/1‑149 (MQ=255)
                                                                                                                             cgcgCTTGCCACGCATGCACCGCCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCCCCTgg                         >  7:58221/1‑149 (MQ=255)
                                                                                                                                    gCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGcc                   <  4:85695/148‑1 (MQ=255)
                                                                                                                                     ccACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTg                 >  5:148540/1‑149 (MQ=255)
                                                                                                                                      cACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTgg                >  4:331529/1‑149 (MQ=255)
                                                                                                                                        cGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGtt              <  4:78863/149‑1 (MQ=255)
                                                                                                                                            tGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc          <  6:162343/149‑1 (MQ=255)
                                                                                                                                             gCACCGGCGCATGAGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc          <  6:247740/148‑1 (MQ=255)
                                                                                                                                              cACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCg         >  3:263997/1‑148 (MQ=255)
                                                                                                                                               aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGt                                                   <  2:145541/105‑1 (MQ=255)
                                                                                                                                               aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGt                                                   >  1:145541/1‑105 (MQ=255)
                                                                                                                                               aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGaa       >  6:202244/1‑149 (MQ=255)
                                                                                                                                                ccGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGaa       >  3:10609/1‑148 (MQ=255)
                                                                                                                                                    cgcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGtt  <  4:57802/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGCTTTAGCCTGGTGTGGCTGGCCATGGTCGGGGCCGACCTGCTGTGTGGCAGCGCCTGGGTACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTT  >  NC_002947/4437032‑4437328

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: