Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 676,098 T→C 8.2% T105A (ACC→GCC)  PP_0578 ← putative Methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947676,0980TC 8.2% 172.3 / 5.0 61T105A (ACC→GCC) PP_0578putative Methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (3/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

CGATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGC  >  NC_002947/675953‑676245
                                                                                                                                                 |                                                                                                                                                   
cGATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAg                                                                                                                                                  <  7:191762/149‑1 (MQ=255)
 gATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAgg                                                                                                                                                 >  2:103757/1‑148 (MQ=255)
   tAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTg                                                                                                                                               <  3:201290/149‑1 (MQ=255)
       ccAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCg                                                                                                                                           <  2:368635/149‑1 (MQ=255)
          gggCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCgagg                                                                                                                                        <  3:268901/149‑1 (MQ=255)
          gggCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGCCAACTGACCgagc                                                                                                                                        <  5:57689/149‑2 (MQ=255)
                 gcgGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCtt                                                                                                                                 <  4:114610/149‑1 (MQ=255)
                     ttGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGc                                                                                                                               >  4:273933/1‑147 (MQ=255)
                       gTAGAAGGACAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAg                                                                                                                            <  8:154776/148‑1 (MQ=255)
                            aGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGcc                                                                                                                      <  1:150120/149‑1 (MQ=255)
                             ggCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCt                                                                                                                     <  3:140320/149‑1 (MQ=255)
                               ccAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGCCAACTGACCgagc                                                                                                                                        >  2:110546/1‑127 (MQ=255)
                               ccAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGCCAACTGACCgagc                                                                                                                                        <  1:110546/128‑2 (MQ=255)
                                cAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCATTTGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGa                                                                                                                  <  1:311547/149‑1 (MQ=255)
                                cAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGGGCCTGGa                                                                                                                  >  8:95951/1‑149 (MQ=255)
                                cAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGa                                                                                                                  >  7:1587/1‑149 (MQ=255)
                                     aTAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGaaa                                                                                                             <  1:197500/149‑1 (MQ=255)
                                     aTAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGaaa                                                                                                             <  1:244755/149‑1 (MQ=255)
                                          caaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgc                                                                                                        >  7:286959/1‑149 (MQ=255)
                                           aaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgcg                                                                                                       <  3:359395/149‑1 (MQ=255)
                                           aaccGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTcgcg                                                                                                       >  4:178423/1‑149 (MQ=255)
                                            accGCCGGGGGCCGTGAGCTGGCGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCgg                                                                                                      >  2:209186/1‑149 (MQ=255)
                                             ccgccgGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCgga                                                                                                     >  1:319441/1‑148 (MQ=255)
                                                        cGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCtt                                                                                          >  2:358549/1‑149 (MQ=255)
                                                        cGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCtg                                                                                          >  4:281519/1‑148 (MQ=255)
                                                         gTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTc                                                                                         <  6:5981/149‑1 (MQ=255)
                                                            aGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCgg                                                                                       <  7:112020/148‑1 (MQ=255)
                                                              cTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCAGCAGGTCGTAGCGCTCGGTCAGTTGGCCGCGGAGGTCTTGCCAGGGCGCCTGGACGAAAGTCGCGGTCTGGCCGGCTTCGGCaa                                                                                    >  7:336854/1‑149 (MQ=255)
                                                                 gTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCaaa                                                                                   <  1:156392/147‑1 (MQ=255)
                                                                 gTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCaaa                                                                                   <  8:198488/147‑1 (MQ=255)
                                                                 gTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGc                                                                                 <  2:100452/149‑1 (MQ=255)
                                                                      gcaCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGcc                                                                            >  3:125531/1‑149 (MQ=255)
                                                                       caCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGccc                                                                           >  7:217130/1‑149 (MQ=255)
                                                                        aCCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCt                                                                          >  6:209239/1‑149 (MQ=255)
                                                                                cctGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACGAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGTGCCGt                                                                  >  3:291492/3‑149 (MQ=255)
                                                                                 aTGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTc                                                                 <  8:217130/149‑1 (MQ=255)
                                                                                  tGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTc                                                                 >  5:1664/1‑148 (MQ=255)
                                                                                  tGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTc                                                                 <  6:1664/148‑1 (MQ=255)
                                                                                     tctcGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCGGAGGAGGTCTTGCCAGGGCGCCTGGAGGAAAGTCGCGGTCTGGCCGGATTCGGCAACGCGTGCCCTGGCGCCGTCGAGc                                                             >  5:206296/1‑149 (MQ=255)
                                                                                         gggTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCAt                                                           <  4:125531/147‑1 (MQ=255)
                                                                                         gggTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATg                                                          <  3:178423/148‑1 (MQ=255)
                                                                                               ggCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGcgccg                                                   >  4:301646/1‑149 (MQ=255)
                                                                                                  cagccaTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAg                                                <  3:31046/149‑1 (MQ=255)
                                                                                                    gccaTTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGt                                              <  1:358549/149‑1 (MQ=255)
                                                                                                           gAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGTTTCGGc                                         >  7:950/1‑147 (MQ=255)
                                                                                                            aGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAg                                      <  8:1587/149‑1 (MQ=255)
                                                                                                             gCACGGCGTGGCAGAGAACCACCAACAACCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCATGGATGAAAGTCGCGGTCTGGCCGGCATCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAgg                                     >  6:175805/1‑149 (MQ=255)
                                                                                                                  gCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAg                                >  4:112748/1‑149 (MQ=255)
                                                                                                                     tGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGCGCGCCGCAGGTTCGGCCAGGGTCAGCtg                             <  1:209186/149‑1 (MQ=255)
                                                                                                                       gCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCtggt                           >  3:354507/1‑149 (MQ=255)
                                                                                                                       gCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCtggt                           >  1:182551/1‑149 (MQ=255)
                                                                                                                         agagCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCtggtgg                         >  7:59440/1‑149 (MQ=255)
                                                                                                                            gCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCcg                      >  5:263677/1‑149 (MQ=255)
                                                                                                                                cAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCt                  >  6:77196/1‑149 (MQ=255)
                                                                                                                                         ctccTCGGCCAACTGACCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCT‑GCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGgccagcca        >  2:314739/4‑149 (MQ=255)
                                                                                                                                          tccTCGGCCAACTGACCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGgccagcca        >  5:129711/3‑149 (MQ=255)
                                                                                                                                              cGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAg    <  3:112748/149‑1 (MQ=255)
                                                                                                                                              cGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAg    <  6:263677/149‑1 (MQ=255)
                                                                                                                                                gTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAgg                                     >  1:252230/1‑114 (MQ=255)
                                                                                                                                                gTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAgg                                     <  2:252230/114‑1 (MQ=255)
                                                                                                                                                gTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGc  >  5:311482/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGATAGACCAGGGCGTCGCGGTTGTAGAAGGCCAGCGATAACCAACCGCCGGGGGCCGTGAGCTGGTGCAGCACCGGCAGGATGCTCTCGGGTTCGGCCAGCCATTCGAGCACGGCGTGGCAGAGCACCAGGTCGTAGCGCTCGGTCAGTTGGCCGAGGAGGTCTTGCCAGGGCGCCTGGATGAAAGTCGCGGTCTGGCCGGCTTCGGCAAAGCGTGCCCTGGCGCCGTCGAGCATGGGCGCCGCAGGTTCGGCCAGGGTCAGCTGGTGGCCGCGCTCGGCCAGCCACAAGGC  >  NC_002947/675953‑676245

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