Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,847,533 T→C 6.1% K59E (AAG→GAG)  PP_5126 ← HDOD domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,847,5330TC 6.1% 194.2 / 2.8 65K59E (AAG→GAG) PP_5126HDOD domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (36/25);  new base C (2/2);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AGGGTAACGATGGTCTCGTCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGCTTGCTGGA  >  NC_002947/5847385‑5847679
                                                                                                                                                    |                                                                                                                                                  
aGGGTAACGAT‑GTCTCGTCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCtt                                                                                                                                                    >  2:102010/1‑148 (MQ=255)
     aaCGATGGTCTCGTCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGc                                                                                                                                                >  7:215338/1‑148 (MQ=255)
      aCGATGGTCTCGTCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTAGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCCAACCCATCTAGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCATGAGc                                                                                                                                                >  1:365638/1‑147 (MQ=255)
        gATGGTCTCGTCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAgg                                                                                                                                            >  6:64923/1‑149 (MQ=255)
                  tCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGCCAAGCCGAAGTGCGGAGAGTGGACCAGCTTGAGCAAGGCGCCGTagag                                                                                                                                  >  8:72623/1‑149 (MQ=255)
                  tCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGagag                                                                                                                                  >  4:346266/1‑149 (MQ=255)
                    aCTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGc                                                                                                                                >  7:46713/1‑149 (MQ=255)
                      tCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGccc                                                                                                                              >  1:309656/1‑149 (MQ=255)
                      tCATCTCGAACTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACATAGCAGGTTCACAGCACGCTGGATCGAACAGATCTTGTTCGACAAGCCGAATTGCCGAGAGTTGACCAGCATGAGCAAGGCGACGGAGATGccc                                                                                                                              >  7:310763/1‑149 (MQ=255)
                        aTCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGccc                                                                                                                              <  5:311/147‑1 (MQ=255)
                        aTCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGccc                                                                                                                              <  4:357754/147‑1 (MQ=255)
                        aTCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCgg                                                                                                                            <  7:72623/149‑1 (MQ=255)
                        aTCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCgg                                                                                                                            >  4:95593/1‑149 (MQ=255)
                                 ttGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTggg                                                                                                                    >  4:262784/1‑148 (MQ=255)
                                 ttGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTggg                                                                                                                    >  3:86101/1‑148 (MQ=255)
                                       gACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCa                                                                                                              <  8:215338/148‑1 (MQ=255)
                                             gCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGc                                                                                                       <  1:109739/149‑1 (MQ=255)
                                                        gTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCtt                                                                                                          <  6:330127/135‑1 (MQ=255)
                                                        gTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCtt                                                                                                          >  5:330127/1‑135 (MQ=255)
                                                        gTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTc                                                                                              >  5:346475/1‑147 (MQ=255)
                                                              gatCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCg                                                                                      >  4:255646/1‑149 (MQ=255)
                                                               atCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCgg                                                                                     <  3:95593/149‑1 (MQ=255)
                                                                tCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTCGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCggg                                                                                    >  3:146305/1‑149 (MQ=255)
                                                                tCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTCGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCggg                                                                                    >  3:151821/1‑149 (MQ=255)
                                                                 cGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAAt                                                                                                                <  8:81633/120‑1 (MQ=255)
                                                                 cGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAAt                                                                                                                >  7:81633/1‑120 (MQ=255)
                                                                  gatcgtcTACCTTGCAGGTTCACCGCACGCTGGATAGATGCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGt                                                                                   <  7:332147/142‑1 (MQ=255)
                                                                  gAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGt                                                                                   >  2:100127/1‑148 (MQ=255)
                                                                    aCGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAg                                                                                >  6:135327/1‑149 (MQ=255)
                                                                      ggCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGc                                                                              <  8:46713/149‑1 (MQ=255)
                                                                      ggCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGc                                                                              <  5:64923/149‑1 (MQ=255)
                                                                        cTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCAt                                                                            <  4:86101/149‑1 (MQ=255)
                                                                         tACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATg                                                                           >  6:39875/1‑149 (MQ=255)
                                                                         tACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATg                                                                           >  5:47693/1‑149 (MQ=255)
                                                                           ccTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTCGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATg                                                                           <  4:146305/147‑1 (MQ=255)
                                                                           ccTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTCGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATg                                                                           <  4:151821/147‑1 (MQ=255)
                                                                             tAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACACGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGCCCCGGGTCCTGGGAAATCAGAATGGCGATGCTTTCCAGGTCCGGGTCAGGCAAGTACt                                                                       >  7:352154/1‑149 (MQ=255)
                                                                                      cACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAAc                                                               <  2:309656/148‑1 (MQ=255)
                                                                                      cACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAAc                                                               <  6:346475/148‑1 (MQ=255)
                                                                                                   aTCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCAt                                                 >  8:270098/1‑149 (MQ=255)
                                                                                                    tCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCggg                                                                                                                           <  5:147447/74‑1 (MQ=255)
                                                                                                    tCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCggg                                                                                                                           >  6:147447/1‑74 (MQ=255)
                                                                                                      gAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGa                                               <  3:346266/148‑1 (MQ=255)
                                                                                                      gAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGAt                                              >  7:347928/1‑149 (MQ=255)
                                                                                                            aTCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGTCCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCgg                                        >  5:70361/1‑149 (MQ=255)
                                                                                                            aTCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCgg                                        >  2:344921/1‑149 (MQ=255)
                                                                                                                 gTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGc                                          <  3:262784/142‑1 (MQ=255)
                                                                                                                 gTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCg                                         <  3:255646/143‑1 (MQ=255)
                                                                                                                  ttCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGACCGGGTCAGGCATGTACTGATCGAACTGCAGATCGACCATGATCTGCGGTTGAgg                                  >  5:350947/1‑149 (MQ=255)
                                                                                                                      aCAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGAt                                              >  3:21877/1‑133 (MQ=255)
                                                                                                                      aCAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGAt                                              <  4:21877/133‑1 (MQ=255)
                                                                                                                        aaGCCGAAGTGCGGAGAGTTGTCCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGAt                            <  6:70361/149‑1 (MQ=255)
                                                                                                                           ccGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGt                         <  5:135327/149‑1 (MQ=255)
                                                                                                                            cGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGt                         >  1:363372/1‑148 (MQ=255)
                                                                                                                                gTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGc                    <  5:39875/149‑1 (MQ=255)
                                                                                                                                  gCGGAGAGTTGACAAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACAATGATCTGCTGTTGAGGCGGGATGGTGATgaat                  >  8:340768/1‑146 (MQ=255)
                                                                                                                                   cGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCtt                 >  6:379240/1‑149 (MQ=255)
                                                                                                                                      agagTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGc               >  6:234563/1‑148 (MQ=255)
                                                                                                                                          ttGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGc          >  5:14525/1‑149 (MQ=255)
                                                                                                                                               cAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCa                                                  >  2:154483/1‑104 (MQ=255)
                                                                                                                                               cAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCa                                                  <  1:154483/104‑1 (MQ=255)
                                                                                                                                               cAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGCTTGCt     <  1:100127/149‑1 (MQ=255)
                                                                                                                                               cAGCTTGAGCAAGGCGCAGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGCTTGCt     >  6:135357/1‑149 (MQ=255)
                                                                                                                                                 gCTTGAGCAAGGGGCCGGAGAGGCCCGGTTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCTTGATCTGCGGTTGAGGCGGGATGTTGATGCCTTGCAAGGGTTGCTgg   <  5:271707/149‑1 (MQ=255)
                                                                                                                                                  cTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGCTTGCTGGa  >  6:89039/1‑149 (MQ=255)
                                                                                                                                                  cTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGCTTGCTGGa  >  5:117429/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGGGTAACGATGGTCTCGTCACTCATCTCGCCCTTGATCGACTGGGCGTTGATCAGGTTGATGATCGAACGGCTACCTAGCAGGTTCACCGCACGCTGGATCGAACCGATCTTGTTCGACAAGCCGAAGTGCGGAGAGTTGACCAGCTTGAGCAAGGCGCCGGAGAGGCCCGGGTCCTGGGAAATCAGCTTGGCGATGGTTTCCAGGTCCGGGTCAGGCATGTACTGCTCGAACTGCAGATCGACCATGATCTGCGGTTGAGGCGGGATGGTGATGCCTTGCAAGGCTTGCTGGA  >  NC_002947/5847385‑5847679

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