Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 782,852 C→A 5.2% R1025L (CGC→CTC)  PP_0672 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947782,8520CA 5.2% 316.3 / 5.6 97R1025L (CGC→CTC) PP_0672sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (45/47);  new base A (2/3);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

Only 100 of 104 total aligned reads displayed.

ATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCACTGCCATCCTG  >  NC_002947/782705‑783000
                                                                                                                                                    |                                                                                                                                                    
aTCCCCGCGACGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGgcg                                                                                                                                                     <  8:387916/149‑1 (MQ=255)
 tCACCGGGTTGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGgcgc                                                                                                                                                    <  1:184908/149‑1 (MQ=255)
   acagggATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTAGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAg                                                                                                                                                  >  6:121222/4‑149 (MQ=255)
    ccGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCGCCATGCCATGGAGCAga                                                                                                                                                 <  8:346500/149‑2 (MQ=255)
     cGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGt                                                                                                                                                >  7:111998/1‑149 (MQ=255)
      gggATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCGCCATGCCATGGAGCAgaac                                                                                                                                               <  2:182153/149‑4 (MQ=255)
       ggATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCTCTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCGCCATGCCATGGAGCAgaac                                                                                                                                               <  1:130642/149‑4 (MQ=255)
       ggATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCg                                                                                                                                              >  8:190679/1‑149 (MQ=255)
          tGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCtt                                                                                                                                           <  6:110544/149‑1 (MQ=255)
             aGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGGATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTc                                                                                                                                         >  5:338226/1‑148 (MQ=255)
             aGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTcc                                                                                                                                        <  2:148442/149‑1 (MQ=255)
              gTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGGTCCATGGCATGGCGCAGGTCGCTTTCCa                                                                                                                                       <  1:246451/149‑1 (MQ=255)
              gTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCAGCTCCATGGCATGGCGCAGGTCGCTTTCCa                                                                                                                                       <  1:39303/149‑1 (MQ=255)
                 gCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGc                                                                                                                                     <  6:245327/148‑1 (MQ=255)
                 gCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCt                                                                                                                                    >  8:171172/1‑149 (MQ=255)
                  cccGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTc                                                                                                                                   <  6:113025/149‑1 (MQ=255)
                      gcggcACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAgg                                                                                                                               <  6:181612/149‑1 (MQ=255)
                      gcggcACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAgg                                                                                                                               <  7:313070/149‑1 (MQ=255)
                       cggcACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGgc                                                                                                                              <  3:194303/149‑1 (MQ=255)
                             cAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTc                                                                                                                          >  1:176420/1‑147 (MQ=255)
                              aGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAg                                                                                                                       >  3:249193/1‑149 (MQ=255)
                               ggCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAgg                                                                                                                      >  2:203925/1‑149 (MQ=255)
                                gCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAggg                                                                                                                     >  7:173040/1‑149 (MQ=255)
                                      ggtcgggtGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTgg                                                                                                               >  4:161110/1‑149 (MQ=255)
                                          gggtGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGtt                                                                                                           <  1:424785/149‑1 (MQ=255)
                                          gggtGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGtt                                                                                                           <  1:424766/149‑1 (MQ=255)
                                            gtGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCa                                                                                                         >  6:29356/1‑149 (MQ=255)
                                                cGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTc                                                                                                     <  6:148581/149‑1 (MQ=255)
                                                cGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTACAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCAttt                                                                                                      <  5:162915/148‑1 (MQ=255)
                                                 gCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCg                                                                                                    >  1:281901/1‑149 (MQ=255)
                                                 gCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCg                                                                                                    <  3:279276/149‑1 (MQ=255)
                                                 gCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTACCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCg                                                                                                    <  7:385598/149‑1 (MQ=255)
                                                   cAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGc                                                                                                  <  6:273645/149‑1 (MQ=255)
                                                    aGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGcc                                                                                                 <  5:20670/149‑1 (MQ=255)
                                                    aGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGcc                                                                                                 >  7:400243/1‑149 (MQ=255)
                                                      cgcAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTg                                                                                               <  4:429817/149‑1 (MQ=255)
                                                      cgcAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTg                                                                                               <  8:173040/149‑1 (MQ=255)
                                                       gcagcaGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTATTAGACGATGAAGTCGTTCTGCTACATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTgg                                                                                              >  3:248572/1‑149 (MQ=255)
                                                       gcagcaGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTgg                                                                                              <  4:249193/149‑1 (MQ=255)
                                                        cagcagCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGt                                                                                             <  2:281901/149‑1 (MQ=255)
                                                         agcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCg                                                                                                             <  1:252512/132‑1 (MQ=255)
                                                         agcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCg                                                                                                             >  2:252512/1‑132 (MQ=255)
                                                         agcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCg                                                                                                             <  1:259125/132‑1 (MQ=255)
                                                           cagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAg                                                                                           >  7:87396/1‑148 (MQ=255)
                                                             gcgcTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCAttt                                                                                                      >  1:33172/1‑135 (MQ=255)
                                                             gcgcTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCAttt                                                                                                      <  2:33172/135‑1 (MQ=255)
                                                                cTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGtt                                                                                                           <  4:256807/127‑1 (MQ=255)
                                                                cTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGtt                                                                                                           >  3:256807/1‑127 (MQ=255)
                                                                cTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTgg                                                                                     <  1:197552/149‑1 (MQ=255)
                                                                   cGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTTGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCTCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGaa                                                                                   <  5:121222/148‑1 (MQ=255)
                                                                     gcgcCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAgtt                                                                                >  6:94884/1‑149 (MQ=255)
                                                                      cgcCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAgttg                                                                               <  3:161110/149‑1 (MQ=255)
                                                                           gTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCtt                                                                          <  2:176420/149‑1 (MQ=255)
                                                                           gTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCtt                                                                          <  5:424493/149‑1 (MQ=255)
                                                                                 cgcTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCcgcg                                                                    <  8:111998/149‑1 (MQ=255)
                                                                                    ttGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTc                                                                 >  2:371029/1‑149 (MQ=255)
                                                                                    ttGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTc                                                                 <  5:29356/149‑1 (MQ=255)
                                                                                    ttGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTc                                                                 >  4:234432/1‑149 (MQ=255)
                                                                                      gCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCt                                                               >  1:190021/1‑149 (MQ=255)
                                                                                          gTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTAGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTgg                                                            >  7:52173/1‑149 (MQ=255)
                                                                                           tCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGc                                                           <  5:173096/149‑1 (MQ=255)
                                                                                           tCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGc                                                           <  5:368158/149‑1 (MQ=255)
                                                                                           tCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGc                                                           >  6:13246/1‑149 (MQ=255)
                                                                                              ccgcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTg                                                        >  6:395080/1‑149 (MQ=255)
                                                                                                 cTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTa                                                     <  4:346245/149‑1 (MQ=255)
                                                                                                   gAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTACAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTAc                                                    <  1:371029/148‑1 (MQ=255)
                                                                                                        gCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCg                                              >  8:150949/1‑149 (MQ=255)
                                                                                                         cGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCg                                              <  5:23753/148‑1 (MQ=255)
                                                                                                             tGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTc                                         >  7:265021/1‑149 (MQ=255)
                                                                                                                tAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGc                                      <  2:358735/149‑1 (MQ=255)
                                                                                                                    agagGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGAGCTACAGGGCGCTGGGGTTAATTTCGGCCTGGTAGAACTGGATGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGttc                                  >  7:121798/1‑149 (MQ=255)
                                                                                                                       ggATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTGCTTGGCGTGGTACATGGCGGTGTCGGCGttcttc                               <  5:69275/149‑1 (MQ=255)
                                                                                                                        gATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTTGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTtca                              <  1:203925/149‑1 (MQ=255)
                                                                                                                        gATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTtca                              >  6:217606/1‑149 (MQ=255)
                                                                                                                        gATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTtca                              >  3:139119/1‑149 (MQ=255)
                                                                                                                          tGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGttc                                  >  2:289188/1‑143 (MQ=255)
                                                                                                                          tGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGttc                                  <  1:289188/143‑1 (MQ=255)
                                                                                                                             aTTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAg                                                                                  >  4:263153/1‑92 (MQ=255)
                                                                                                                             aTTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAg                                                                                  <  3:263153/92‑1 (MQ=255)
                                                                                                                             aTTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGc                         >  7:231479/1‑149 (MQ=255)
                                                                                                                             aTTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGc                         >  5:386553/1‑149 (MQ=255)
                                                                                                                             aTTCGGTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGc                         <  4:263373/149‑1 (MQ=255)
                                                                                                                              ttCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCt                        >  7:101726/1‑149 (MQ=255)
                                                                                                                                 gTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTgg                      <  5:13246/148‑1 (MQ=255)
                                                                                                                                 gTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGgc                     >  1:129188/1‑149 (MQ=255)
                                                                                                                                  ttCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCcgc                                                                     <  6:204815/100‑1 (MQ=255)
                                                                                                                                  ttCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCcgc                                                                     >  5:204815/1‑100 (MQ=255)
                                                                                                                                  accTGCGCCATGCCATGGAGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGgct                    >  2:362917/3‑149 (MQ=255)
                                                                                                                                  acatGCGCCATGCCATGGAGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGgct                    >  1:99409/4‑149 (MQ=255)
                                                                                                                                      gCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGAACATGGCGGTGTCGGCGTTCTTCATCAGCTGgctcag                 >  7:100931/1‑148 (MQ=255)
                                                                                                                                       cTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTg                                                                         <  4:423923/91‑1 (MQ=255)
                                                                                                                                       cTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTg                                                                         >  3:423923/1‑91 (MQ=255)
                                                                                                                                       cTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGgctcagct               >  8:134694/1‑149 (MQ=255)
                                                                                                                                        tCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGgctcagct               >  5:88747/1‑148 (MQ=255)
                                                                                                                                        tCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGgctcagct               >  3:119831/1‑148 (MQ=255)
                                                                                                                                          catggcatggCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCAc            >  4:188860/1‑149 (MQ=255)
                                                                                                                                               catggcGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCACTGCCa       >  5:50873/1‑149 (MQ=255)
                                                                                                                                               catggcGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCACTGCCa       <  2:129188/149‑1 (MQ=255)
                                                                                                                                                atggcGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCACTGCCAt      >  3:407834/1‑149 (MQ=255)
                                                                                                                                                    cgcAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTCGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCACTGCCATCCTg  <  2:190021/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCGCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGGCTCAGCTCACTGCCATCCTG  >  NC_002947/782705‑783000

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: