Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,254,842 T→A 13.8% K299M (AAG→ATG) ‡ PP_2852 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,254,8420TA13.8% 167.7 / 16.2 65K299M (AAG→ATG) ‡PP_2852sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (3/6);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ATCGAAACCGATGTGCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCTTCGT  >  NC_002947/3254694‑3254986
                                                                                                                                                    |                                                                                                                                                 
aTCGAAACCGATGTGCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCt                                                                                                                                                   >  5:358190/1‑149 (MQ=255)
         gATGTGCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTgg                                                                                                                                          <  1:35295/149‑2 (MQ=255)
           tgtgCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCg                                                                                                                                        <  6:4967/149‑1 (MQ=255)
           tgtgCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTggaa                                                                                                                                        <  2:84054/149‑4 (MQ=255)
             tgCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGtt                                                                                                                                      <  2:300068/149‑1 (MQ=255)
             tgCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTggaac                                                                                                                                       <  7:166408/148‑5 (MQ=255)
             tgCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTggaac                                                                                                                                       <  6:295249/148‑5 (MQ=255)
                 gCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGc                                                                                                                                   <  7:209547/148‑1 (MQ=255)
                 gCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGc                                                                                                                                   >  8:209547/1‑148 (MQ=255)
                     gatgCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCt                                                                                                                              >  7:63829/1‑149 (MQ=255)
                       tgCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTgg                                                                                                                                          <  5:301743/135‑2 (MQ=255)
                       tgCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTgg                                                                                                                                          >  6:301743/1‑134 (MQ=255)
                        gCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGt                                                                                                                           >  4:291942/1‑149 (MQ=255)
                         cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  6:358190/149‑1 (MQ=255)
                         cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  8:280460/149‑1 (MQ=255)
                         cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  4:241121/149‑1 (MQ=255)
                         cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  6:195402/149‑1 (MQ=255)
                             ccTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTaa                                                                                                                      >  5:67140/1‑149 (MQ=255)
                               ttGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTggaa                                                                                                                                        >  2:46268/1‑126 (MQ=255)
                               ttGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTggaa                                                                                                                                        <  1:46268/129‑4 (MQ=255)
                                tGGCCAGGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATcg                                                                                                                   >  4:47581/1‑149 (MQ=255)
                                  gCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTTCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATcgcc                                                                                                                 <  3:164642/149‑1 (MQ=255)
                                   ccATGAAGCGCAGGCCTGTGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATcgccg                                                                                                                <  4:322711/149‑1 (MQ=255)
                                    cATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGc                                                                                                                               >  1:55308/1‑133 (MQ=255)
                                    cATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGc                                                                                                                               <  2:55308/133‑1 (MQ=255)
                                        aaGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGcaca                                                                                                            <  8:72347/148‑1 (MQ=255)
                                        aaGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGcaca                                                                                                            <  8:63829/148‑1 (MQ=255)
                                              aGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGcc                                                                                                                                  <  3:85577/120‑1 (MQ=255)
                                              aGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGcc                                                                                                                                  >  4:85577/1‑120 (MQ=255)
                                              aGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCg                                                                                                      >  8:54724/1‑148 (MQ=255)
                                                               gTAGTGGGTGGTAAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATgg                                                                                   >  4:109254/1‑149 (MQ=255)
                                                                tAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAACATCGCATACAGGCCACG‑GATGGt                                                                                  <  7:54724/149‑1 (MQ=255)
                                                                    gggtggTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCtgg                                                                              >  4:331211/1‑149 (MQ=255)
                                                                      gtggtgAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCtggtg                                                                            <  3:331211/149‑1 (MQ=255)
                                                                       tggtgAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCtggtgg                                                                           <  6:67140/149‑1 (MQ=255)
                                                                       tggtgAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCtggtgg                                                                           <  3:161303/149‑1 (MQ=255)
                                                                       tggtgAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAACTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCtggtgg                                                                           >  1:371317/1‑149 (MQ=255)
                                                                            aaaCCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGtgtg                                                                      <  3:47581/149‑1 (MQ=255)
                                                                            aaaCCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATCGTCTGGTGGGtgtg                                                                      >  4:229648/1‑149 (MQ=255)
                                                                              aCCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGtgtg                                                                      >  2:85283/1‑147 (MQ=255)
                                                                                   ccTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GTTGGTCTGGTGGGTGTGCAGCACg                                                               >  2:216395/1‑149 (MQ=255)
                                                                                       gCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAg                                                            >  8:142930/1‑148 (MQ=255)
                                                                                       gCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGt                                                           <  3:229648/149‑1 (MQ=255)
                                                                                         gTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCt                                                         >  4:32130/1‑149 (MQ=255)
                                                                                          tAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTg                                                        <  3:109254/149‑1 (MQ=255)
                                                                                            aTTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTggg                                                      <  7:142930/149‑1 (MQ=255)
                                                                                            aTTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTggg                                                      <  2:198195/149‑1 (MQ=255)
                                                                                                   cGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GTTGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCa                                               <  1:216395/149‑1 (MQ=255)
                                                                                                         aCAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGc                                                                                                          >  2:205064/1‑85 (MQ=255)
                                                                                                                       gTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCCCG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTg                           >  7:257872/1‑149 (MQ=255)
                                                                                                                       gTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACGCGATGGTCTGGTGGGTGTGCAGCACGTAGTCt                                                         <  1:38977/120‑1 (MQ=255)
                                                                                                                       gTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCt                                                         >  2:38977/1‑119 (MQ=255)
                                                                                                                       gTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTg                           >  8:81296/1‑149 (MQ=255)
                                                                                                                           gTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGa                       <  7:81296/149‑1 (MQ=255)
                                                                                                                            ttCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGTAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAg                      <  3:291942/149‑1 (MQ=255)
                                                                                                                            ttCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAg                      >  2:140768/1‑149 (MQ=255)
                                                                                                                                    cagcagTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGc               <  2:285781/148‑1 (MQ=255)
                                                                                                                                    cagcagTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGc               >  1:285781/1‑148 (MQ=255)
                                                                                                                                       cagTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATc            >  2:282459/1‑148 (MQ=255)
                                                                                                                                        aacggcACGCCCAAGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATc            >  2:1017/6‑147 (MQ=255)
                                                                                                                                             cACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGc                                  >  7:74128/1‑120 (MQ=255)
                                                                                                                                             cACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGc                                  <  8:74128/120‑1 (MQ=255)
                                                                                                                                             cACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCtt     >  6:329983/1‑149 (MQ=255)
                                                                                                                                             cACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCtt     >  2:81611/1‑149 (MQ=255)
                                                                                                                                             cACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCtt     >  1:336220/1‑149 (MQ=255)
                                                                                                                                                gCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCTTCGt  >  3:328521/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCGAAACCGATGTGCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCTTCGT  >  NC_002947/3254694‑3254986

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