Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,437,241 T→C 5.5% V134A (GTG→GCG)  PP_3932 → diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,437,2410TC 5.5% 215.4 / 2.5 74V134A (GTG→GCG) PP_3932diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (39/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.53e-01

TACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCTTG  >  NC_002947/4437094‑4437386
                                                                                                                                                   |                                                                                                                                                  
tACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggtg                                                                                                                                                   <  6:36955/149‑1 (MQ=255)
tACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggtg                                                                                                                                                   >  5:245000/1‑149 (MQ=255)
  ccGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAtggtgg                                                                                                                                                  <  7:48231/148‑1 (MQ=255)
     gCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTggcggc                                                                                                                                              >  7:178828/1‑149 (MQ=255)
             gTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATgg                                                                                                                                      <  8:321222/149‑1 (MQ=255)
             gTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATgg                                                                                                                                      >  3:16955/1‑149 (MQ=255)
              ttAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATgg                                                                                                                                      >  3:191048/1‑148 (MQ=255)
                    cTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGt                                                                                                                               <  8:178828/149‑1 (MQ=255)
                     tgtgCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGt                                                                                                                               >  7:43311/1‑148 (MQ=255)
                        gCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATg                                                                                                                           >  2:366228/1‑149 (MQ=255)
                         cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATgg                                                                                                                          <  4:351407/149‑1 (MQ=255)
                         cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATgg                                                                                                                          >  2:233215/1‑149 (MQ=255)
                         cAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATgg                                                                                                                          >  7:180712/1‑149 (MQ=255)
                          aaTCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCAAGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGc                                                                                                                         >  2:365452/1‑149 (MQ=255)
                              ttGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCa                                                                                                                     >  8:161241/1‑149 (MQ=255)
                                       tggtCCCCCGCTGGCAGCGGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGGTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTc                                                                                                            <  8:155357/149‑1 (MQ=255)
                                       tggtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTc                                                                                                            <  1:366228/149‑1 (MQ=255)
                                       tggtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTc                                                                                                            <  3:57534/149‑1 (MQ=255)
                                         gtCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGa                                                                                                          <  6:55718/149‑1 (MQ=255)
                                          tCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGaa                                                                                                         >  7:261541/1‑149 (MQ=255)
                                            aCCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACa                                                                                                       <  1:233215/149‑1 (MQ=255)
                                             cccTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAg                                                                                                      >  5:20890/1‑149 (MQ=255)
                                                    gCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGt                                                                                                  <  3:311239/146‑1 (MQ=255)
                                                    gCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTg                                                                                                <  4:275247/148‑1 (MQ=255)
                                                      aGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGGGTACTCAGCGTGTTCGGCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTgcgc                                                                                             <  6:245000/149‑1 (MQ=255)
                                                      aGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGAAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGAGGTAATGCCCGCACCCTGGATTGAACAGTCGTTgcgc                                                                                             >  2:27898/1‑149 (MQ=255)
                                                           ctgTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCAc                                                                                        >  7:138765/1‑149 (MQ=255)
                                                              tCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTg                                                                                     >  2:174805/1‑149 (MQ=255)
                                                                  gcTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTa                                                                                 <  1:174805/149‑1 (MQ=255)
                                                                    ttGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                                                                            >  2:338278/1‑104 (MQ=255)
                                                                    ttGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                                                                            <  1:338278/104‑1 (MQ=255)
                                                                       ccACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTg                                                                            >  5:21805/1‑149 (MQ=255)
                                                                        cACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTgg                                                                           >  7:219536/1‑149 (MQ=255)
                                                                          cGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCAGGCTAGCCAGGtg                                                                         <  4:191048/149‑2 (MQ=255)
                                                                             attcACCGTCGCATTCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCg                                                                                                           <  8:213083/109‑1 (MQ=255)
                                                                             aTGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCAAGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTAGATGCGCGCCACCTGGCTAGCCTGGTTCAg                                                                      >  3:116674/1‑149 (MQ=255)
                                                                             aTGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAACGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTgcgcgc                                                                                           >  7:329715/1‑128 (MQ=255)
                                                                             aTGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAACGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTgcgcgc                                                                                           <  8:329715/128‑1 (MQ=255)
                                                                             aTGCACCGGCGCATGCGCAGCCCGCAAGGAGTACTCAGCGTGATCCCCGCCTGAGACGCCCCTGCCATGCAGGCGGCGAGCATGGCAGCGGTAATGGCCGCACACTGGTTCg                                                                                                           >  7:213083/1‑112 (MQ=255)
                                                                              tGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                     <  8:138765/149‑1 (MQ=255)
                                                                              tGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                     <  6:20890/149‑1 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCAt                                                                                                                                        <  7:338501/79‑1 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCAt                                                                                                                                        <  5:131954/79‑1 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCAt                                                                                                                                        >  8:338501/1‑79 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCAt                                                                                                                                        >  6:131954/1‑79 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                                                                            <  4:312961/91‑1 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAAt                                                                                                                            >  3:312961/1‑91 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                     >  1:360451/1‑146 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGc                                                                     <  2:360451/146‑1 (MQ=255)
                                                                                 aCCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGaa                                                                  <  8:43311/149‑1 (MQ=255)
                                                                                   cGGTGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGaa                                                                  >  3:282222/1‑147 (MQ=255)
                                                                                   cGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGaa                                                                  >  4:88717/1‑147 (MQ=255)
                                                                                       gcATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTtc                                                           <  3:88717/149‑1 (MQ=255)
                                                                                             gcAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGAcc                                                     >  8:351164/1‑149 (MQ=255)
                                                                                                gCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGc                                                  <  6:66890/149‑1 (MQ=255)
                                                                                                    gCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCgtgc                                              >  3:182028/1‑149 (MQ=255)
                                                                                                             aCTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGtgctgc                                     >  1:228313/1‑149 (MQ=255)
                                                                                                               tCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGc                                                                               >  7:58187/1‑106 (MQ=255)
                                                                                                               tCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGc                                                                               <  8:58187/106‑1 (MQ=255)
                                                                                                                  gCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtg                                >  8:101355/1‑149 (MQ=255)
                                                                                                                   cGTGTTGGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGGGGTAATGGCCGCACCCTTGTTTGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACATTATCTCGACCGGCGTGCTGGTGCTGCCGGtg                                <  1:27898/149‑1 (MQ=255)
                                                                                                                    gtgtTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtgct                              >  2:259473/1‑149 (MQ=255)
                                                                                                                        gagccgccTTCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTgctgct                           <  4:16955/146‑1 (MQ=255)
                                                                                                                          gccgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGc                                                                               >  8:260402/1‑95 (MQ=255)
                                                                                                                          gccgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGc                                                                               <  7:260402/95‑1 (MQ=255)
                                                                                                                          gccgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTa                         >  8:272667/1‑148 (MQ=255)
                                                                                                                          gccgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTa                         <  6:21805/148‑1 (MQ=255)
                                                                                                                           ccgccTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTa                         >  8:23370/1‑147 (MQ=255)
                                                                                                                                    tcgccgcCACCATGGCAGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtgc                               <  4:358624/131‑1 (MQ=255)
                                                                                                                                    tcgccgcCACCATGGCAGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGtgc                               >  3:358624/2‑132 (MQ=255)
                                                                                                                                    tcgccgcCACCATGGCAGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTc               >  5:65396/2‑148 (MQ=255)
                                                                                                                                      gccgccACCATGGCAGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTc               <  6:65396/146‑1 (MQ=255)
                                                                                                                                        cgcTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGc          <  7:101355/149‑1 (MQ=255)
                                                                                                                                          cTGCCATGGTGTCGGCGAGCATGGCAGCGGTAATGGCCGCACGCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGTACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACATCGCCGTCTGTGCgg        <  4:116674/149‑1 (MQ=255)
                                                                                                                                          cTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCgg        <  8:261541/149‑1 (MQ=255)
                                                                                                                                           tGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCggg       >  3:114616/1‑149 (MQ=255)
                                                                                                                                               aTGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCtt   >  2:149962/1‑149 (MQ=255)
                                                                                                                                                tggtggCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCTTg  <  7:23370/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACCGGCTCTGTGGTTAAACCTGTGCAATCTTGGCGTGGTGGTCACCCTCTGGCAGCTGCTGTCGCGCTTGCCACGCATGCACCGGCGCATGCGCAGCCCGCACGGCGTACTCAGCGTGTTCGCCGCCTGCGCCGCCGCTGCCATGGTGGCGGCGAGCATGGCAGCGGTAATGGCCGCACCCTGGTTCGAACAGTCGTTGCGCGCCACCTGGCTAGCCTGGTTCAGCGAACAGTTCTCGACCGGCGTGCTGGTGCTGCCGGTGCTGCTTACAGCGCCCTCGGTGCGGGCCTTG  >  NC_002947/4437094‑4437386

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: