Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,254,843 T→A 13.8% K299* (AAG→TAG) ‡ 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,8430TA13.8% 167.4 / 16.5 65K299* (AAG→TAG) ‡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 = 7.29e-01

GATGTGCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCTTCGT  >  NC_002947/3254703‑3254986
                                                                                                                                            |                                                                                                                                                
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                                                                                                                                          >  6:301743/1‑134 (MQ=255)
              tgCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCAAGGGCGTgg                                                                                                                                          <  5:301743/135‑2 (MQ=255)
               gCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGt                                                                                                                           >  4:291942/1‑149 (MQ=255)
                cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  6:358190/149‑1 (MQ=255)
                cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  4:241121/149‑1 (MQ=255)
                cGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTc                                                                                                                          <  8:280460/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                                                                                                                               <  2:55308/133‑1 (MQ=255)
                           cATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGc                                                                                                                               >  1:55308/1‑133 (MQ=255)
                               aaGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGcaca                                                                                                            <  8:63829/148‑1 (MQ=255)
                               aaGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGcaca                                                                                                            <  8:72347/148‑1 (MQ=255)
                                     aGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGcc                                                                                                                                  >  4:85577/1‑120 (MQ=255)
                                     aGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGcc                                                                                                                                  <  3:85577/120‑1 (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)
                                                                                                                                            |                                                                                                                                                
GATGTGCGGCATGATGCGGTCCTTGGCCATGAAGCGCAGGCCTGCGCCTTGCAGGTAGTGGGTGGTGAAACCGGCCTGGCGTAATTGTGCCGGCAGACAGGCCTGGTTGCGTTGGTTCTGGGTCAGCAGTTCCACGCCCTTGGGCGTGCCGTTGGCCAGCTTGTCGTAATCGCCGCACAGCATCGCATACAGGCCACG‑GATGGTCTGGTGGGTGTGCAGCACGTAGTCTGGGGTGTTCATGCCGCGCTCGGCCCATTTGCTGAGCTTGGGCATCAAGTCTTCGT  >  NC_002947/3254703‑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: ‑