Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,033,782 G→C 7.5% Q35H (CAG→CAC)  dfp → bifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,033,7820GC 7.5% 204.4 / 6.4 67Q35H (CAG→CAC) dfpbifunctional 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base C (3/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

ATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCC  >  NC_002947/6033634‑6033917
                                                                                                                                                     |                                                                                                                                       
aTTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGCCAAGCCTATGCAGCGGCTGTATCGAAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAg                                                                                                                                         >  4:305524/1‑149 (MQ=255)
     aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgcg                                                                                                                                    >  4:270700/1‑149 (MQ=255)
     aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgc                                                                                                                                     <  3:82929/148‑1 (MQ=255)
     aaaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgc                                                                                                                                     <  4:190931/148‑1 (MQ=255)
      aaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTgcgc                                                                                                                                   >  4:12860/1‑149 (MQ=255)
      aaCGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTgcgc                                                                                                                                   >  2:309896/1‑149 (MQ=255)
         gtcgtcGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCGCCCGTCTCCTGGAGCACGGCGCGCAGGTGCGcgt                                                                                                                                 <  1:47920/149‑1 (MQ=255)
           cgtcgCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTGCGcgccgt                                                                                                                              <  1:309896/149‑1 (MQ=255)
                 ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                                             >  4:58424/1‑144 (MQ=255)
                 ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCgtcgtc                                                                                                                             <  3:58424/144‑1 (MQ=255)
                 ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGAc                                                                                                                        >  5:57272/1‑149 (MQ=255)
                 ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTGCGcgccgtgc                                                                                                                            >  6:138448/1‑143 (MQ=255)
                 ccGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGACGCGCACCTGCGcgccgtgc                                                                                                                            <  5:138448/145‑3 (MQ=255)
                     tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGcgtc                                                                                                                                <  4:141464/137‑1 (MQ=255)
                     tctcACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGcgtc                                                                                                                                >  3:141464/1‑137 (MQ=255)
                         aCTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtg                                                                                                                >  1:299987/1‑149 (MQ=255)
                           tGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtggt                                                                                                              >  8:340005/1‑149 (MQ=255)
                             gAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACGAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGgtggtg                                                                                                             >  8:117467/1‑148 (MQ=255)
                             gAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGggtggt                                                                                                              <  3:12860/147‑1 (MQ=255)
                              aaaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCa                                                                                                           <  5:341138/149‑1 (MQ=255)
                               aaaaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCGTGACCCGGGGTGGTGCa                                                                                                           <  4:343676/148‑1 (MQ=255)
                                 aaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGaa                                                                                                        <  3:305524/149‑1 (MQ=255)
                                 aaGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGaa                                                                                                        >  2:145218/1‑149 (MQ=255)
                                      aaGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATtcat                                                                                                   <  3:270700/149‑1 (MQ=255)
                                      aaGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATtcat                                                                                                   <  1:166490/149‑1 (MQ=255)
                                            aTGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCAcccc                                                                                             >  5:240827/1‑149 (MQ=255)
                                            aTGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CAGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACACGGGGTGGTGCAGAATTCATCAcccc                                                                                             >  5:246061/1‑149 (MQ=255)
                                                aGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCt                                                                                          >  2:113085/1‑148 (MQ=255)
                                                 gCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGa                                                                                        >  2:247614/1‑149 (MQ=255)
                                                   ggCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGAcc                                                                                      <  4:125224/149‑1 (MQ=255)
                                                        tATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGc                                                                                  <  1:145218/148‑1 (MQ=255)
                                                          tCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAgg                                                                               <  6:3784/149‑1 (MQ=255)
                                                              aaGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCt                                                                           <  2:129743/149‑1 (MQ=255)
                                                                gcgcATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGt                                                                         <  6:57272/149‑1 (MQ=255)
                                                                    aTCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                   <  4:198417/103‑1 (MQ=255)
                                                                    aTCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCggg                                                                                                                   >  3:198417/1‑103 (MQ=255)
                                                                             ggCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGt                                                            <  6:246061/149‑1 (MQ=255)
                                                                             ggCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGt                                                            <  6:240827/149‑1 (MQ=255)
                                                                                    gtggCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATg                                                     <  7:340005/149‑1 (MQ=255)
                                                                                       gcggcATTGCCGCCTACGAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  <  7:117467/149‑1 (MQ=255)
                                                                                       gcggcATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCCTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  <  1:113085/149‑1 (MQ=255)
                                                                                         ggcAATGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  <  1:176359/147‑1 (MQ=255)
                                                                                         ggcAATGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGAc                                                  >  2:176359/1‑147 (MQ=255)
                                                                                            aTTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCt                                             <  3:13383/149‑1 (MQ=255)
                                                                                            aTTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCt                                             <  1:247614/149‑1 (MQ=255)
                                                                                                  gccTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGAccc                                       >  3:17539/1‑149 (MQ=255)
                                                                                                      aCAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGgcc                                   >  8:128966/1‑149 (MQ=255)
                                                                                                       cAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGgccg                                  >  7:118134/1‑149 (MQ=255)
                                                                                                         aaaGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGgccgcc                                >  3:31141/1‑149 (MQ=255)
                                                                                                         aaaGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGAGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGCCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCCgccgcc                                >  5:206963/1‑149 (MQ=255)
                                                                                                          aaGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAAGTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGGACACATGGACCTGCTAGACCCGgccg                                  >  8:263070/1‑146 (MQ=255)
                                                                                                                  aGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACAt                                                      <  8:181143/118‑1 (MQ=255)
                                                                                                                  aGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACAt                                                      >  7:181143/1‑118 (MQ=255)
                                                                                                                    cTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCAt                     <  7:128966/149‑1 (MQ=255)
                                                                                                                        ttCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGc                 <  4:17539/149‑1 (MQ=255)
                                                                                                                           g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGc                 <  6:121012/146‑1 (MQ=255)
                                                                                                                           g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGc                 >  5:121012/1‑146 (MQ=255)
                                                                                                                           g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCa               >  1:65922/1‑148 (MQ=255)
                                                                                                                           g‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCa               >  3:258372/1‑148 (MQ=255)
                                                                                                                             ccGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCata             <  8:118134/149‑1 (MQ=255)
                                                                                                                             cactctcCTGGAGCACGGCGCGCAGCTGCGCGTCGTCAAGACCCGGGGTGGTGCAGAATTAAACACCCCGCTGACCCTGCAGGCGCTCTCCCGCCACCCCGAACACATGGACCTGCTCGACCCGGCAGCCGAACCTGCCATGGGCCata             >  6:150123/4‑149 (MQ=255)
                                                                                                                                    cTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGa         <  3:154783/146‑1 (MQ=255)
                                                                                                                                    cTGGAGCACGGCGCGCAGGTGCGCGGCGTCATGACCCGGGGTTGTGCAGATTTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCt      <  6:206963/149‑1 (MQ=255)
                                                                                                                                       gAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc   >  6:174357/1‑149 (MQ=255)
                                                                                                                                       gAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc   <  2:299987/149‑1 (MQ=255)
                                                                                                                                        atgacgacgCGCACCTGCGCGCCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGc   >  3:112219/4‑148 (MQ=255)
                                                                                                                                          cACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGCCCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGcc  >  8:308303/1‑147 (MQ=255)
                                                                                                                                          cACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGcc  >  4:219504/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                       
ATTTCAAACGTCGTCGCCCGCTCTCACTGGAAAAAGGCAAGCCTATGCAGCGGCTGTATCGCAAGCGCATCGTTCTCGGCGTGGGTGGCGGCATTGCCGCCTACAAAAGCGCCGAGCTGATTCG‑CCGTCTCCTGGAGCACGGCGCGCAGGTGCGCGTCGTCATGACCCGGGGTGGTGCAGAATTCATCACCCCGCTGACCCTGCAGGCGCTGTCGGGCCACCCCGTACACATGGACCTGCTCGACCCGGCCGCCGAAGCTGCCATGGGCCATATCGAGCTGGCC  >  NC_002947/6033634‑6033917

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