Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,489,983 G→T 6.0% A484D (GCC→GAC)  clpV ← protein ClpV1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,489,9830GT 6.0% 209.6 / 2.5 66A484D (GCC→GAC) clpVprotein ClpV1
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base T (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATGC  >  NC_002947/3489845‑3490130
                                                                                                                                          |                                                                                                                                                    
aGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGc                                                                                                                                           >  1:348408/1‑150 (MQ=255)
  cTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGa                                                                                                                                         >  1:1025299/1‑150 (MQ=255)
   gccgCACGCGGGGTATTGTGCTCCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCAGTTCGACCAGTTCCAGCTCCTGCTGCCACTGCGCATTGAGCTTGTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGAt                                                                                                                                        <  2:596814/147‑1 (MQ=255)
      gCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCTAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCAGCCACTGGGCATTGAGCTCCTGTTCCTGCTGCTGCGCGGCCTGCAGGGCAGCCTGCAGCGATTcc                                                                                                                                     >  2:215598/1‑150 (MQ=255)
             gggTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAggcggc                                                                                                                              >  2:883928/1‑150 (MQ=255)
               gTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGcggcg‑g                                                                                                                           <  1:269724/150‑1 (MQ=255)
                 aTTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGc                                                                                                                         <  2:348408/150‑1 (MQ=255)
                              gCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGcccc                                                                                                            >  1:1064944/1‑150 (MQ=255)
                                   ggCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGcc                                                                                                       >  1:222021/1‑150 (MQ=255)
                                     cATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGATGCTGCGCCTGCTGCTGCGCGTCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCgg                                                                                                     <  1:883928/150‑1 (MQ=255)
                                      aTCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGc                                                                                                    <  2:815324/150‑1 (MQ=255)
                                      aTCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTACTGGTGTTTGTCGTCCTGCAGGGGAGCCTGAAGCGATTAGCGGCAACC‑TGCAGGGTCGTGCCCCTTGCCGGc                                                                                                    >  1:1096500/1‑150 (MQ=255)
                                         gTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTc                                                                                                 >  2:995161/1‑150 (MQ=255)
                                             gcggcgTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTgg                                                                                             <  2:900841/150‑1 (MQ=255)
                                                  gTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAg                                                                                        >  2:1072269/1‑150 (MQ=255)
                                                  gTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAg                                                                                        <  1:677125/150‑1 (MQ=255)
                                                      aGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACg                                                                                    >  1:511049/1‑150 (MQ=255)
                                                        cgcCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGtg                                                                                  <  1:760822/150‑2 (MQ=255)
                                                          ccTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGa                                                                                >  2:758321/1‑150 (MQ=255)
                                                            ttGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGAtt                                                                              >  1:765397/1‑150 (MQ=255)
                                                                    ttCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGt                                                                                               >  2:754810/1‑125 (MQ=255)
                                                                    ttCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGt                                                                                               <  1:754810/125‑1 (MQ=255)
                                                                       gACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAg                                                                   >  2:501430/1‑150 (MQ=255)
                                                                        aCCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAgg                                                                  >  1:376599/1‑150 (MQ=255)
                                                                         ccAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGTCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCt                                                                       >  1:390254/1‑144 (MQ=255)
                                                                         ccAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGTCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCt                                                                       <  2:390254/144‑1 (MQ=255)
                                                                         ccAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGt                                                                 >  2:724986/1‑150 (MQ=255)
                                                                         ccAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCAGCCTGCAGGTCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCt                                                                       <  2:390360/144‑1 (MQ=255)
                                                                         ccAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCAGCCTGCAGGTCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCt                                                                       >  1:390360/1‑144 (MQ=255)
                                                                          cAGTCCCAGCTCCTGCTGCCAATGGGCATCGAGCTGCTGTTCCTGCTGGTGCTCTTACCGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGtt                                                                <  1:14772/150‑1 (MQ=255)
                                                                           aGTTCCAGCTCCTGCTGCCACTGGGCATTTAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTg                                                               <  1:1072269/150‑1 (MQ=255)
                                                                               ccAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGa                                                           <  1:758321/150‑1 (MQ=255)
                                                                                cAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGc                                                                                                               <  1:124391/97‑1 (MQ=255)
                                                                                cAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGc                                                                                                               >  2:124391/1‑97 (MQ=255)
                                                                                cAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGAt                                                          >  1:277259/1‑150 (MQ=255)
                                                                                 aGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATc                                                         <  1:1075908/150‑1 (MQ=255)
                                                                                   cTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTg                                                       >  1:1356693/1‑150 (MQ=255)
                                                                                    tCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGc                                                      <  1:1377164/150‑1 (MQ=255)
                                                                                     ccTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGcc                                                     >  2:637123/1‑150 (MQ=255)
                                                                                       tgctgcCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGc                                                                         >  1:423633/1‑128 (MQ=255)
                                                                                       tgctgcCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGc                                                                         <  2:423633/128‑1 (MQ=255)
                                                                                           gcCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCggcggc                                               >  2:165728/1‑150 (MQ=255)
                                                                                            ccACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCa                                              <  1:1103704/150‑1 (MQ=255)
                                                                                            ccACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCC‑GCAGACGTGTCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCa                                              <  1:799240/150‑1 (MQ=255)
                                                                                              aCTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGt                                            >  2:231846/1‑150 (MQ=255)
                                                                                                  ggCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTc                                        <  2:1096500/150‑1 (MQ=255)
                                                                                                    cATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCa                                      <  2:1053383/150‑1 (MQ=255)
                                                                                                      ttGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGt                                    <  2:1356693/150‑1 (MQ=255)
                                                                                                      ttGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGt                                    <  2:222021/150‑1 (MQ=255)
                                                                                                      ttGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGt                                    <  1:464810/150‑1 (MQ=255)
                                                                                                          gctgctGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTgcgc                                <  1:135037/150‑1 (MQ=255)
                                                                                                           ctgctgTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTGGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGcc                               >  1:1260685/1‑150 (MQ=255)
                                                                                                            tgctgTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCa                                                                                                                                             <  2:36811/40‑1 (MQ=255)
                                                                                                            tgctgTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCa                                                                                                                                             >  1:36811/1‑40 (MQ=255)
                                                                                                                gTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCa                          <  2:277259/150‑1 (MQ=255)
                                                                                                                 ttCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCcgg                                                   >  1:1238576/1‑124 (MQ=255)
                                                                                                                 ttCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCcgg                                                   <  2:1238576/124‑1 (MQ=255)
                                                                                                                 ttCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAg                         >  2:1281927/1‑150 (MQ=255)
                                                                                                                   ccTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCg                       >  2:1057722/1‑150 (MQ=255)
                                                                                                                     tgctgctgcTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGaa                     <  2:893473/150‑1 (MQ=255)
                                                                                                                       ctgctgctCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCCGCAGGTTGCCGATCTGCCGGCGGCAGTCTTCCAGAGAGCCGGGCAGCGAAGa                   <  2:765397/150‑1 (MQ=255)
                                                                                                                           tgctCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGacaag               >  2:195229/1‑146 (MQ=255)
                                                                                                                            gctCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGGCGTGCCCCTTGGCGGATTCATGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGCTCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGATTgcg              >  2:309315/1‑150 (MQ=255)
                                                                                                                                  ccTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTTGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGATTGCGCCAggg        >  2:1410454/1‑150 (MQ=255)
                                                                                                                                     gCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGAGTGCGGCAAGGGaa     >  2:6892/1‑150 (MQ=255)
                                                                                                                                      cAGGGCAGCATGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATCTACGCCTGCAGGTTGTCGATCTGGCGGCGGCAGTCTGCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAAt    <  1:215598/150‑1 (MQ=255)
                                                                                                                                        gggCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATgc  <  1:995161/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
AGCTCGGCACGCAGGGTATTGAGCTGCTGGGCGTCGGCATCGTTGGCGGCGTCCAGCGCCTTGAGCTGTTCGACCAGTTCCAGCTCCTGCTGCCACTGGGCATTGAGCTGCTGTTCCTGCTGCTGCTCGGCCTGCAGGGCAGCCTGCAGCGATTCCAGGCGGCG‑GGCATGGTCGTGCCCCTTGCCGGCTTCGTGGCCGAGCACGTCGATTTCCGCCTGCAGGTTGTCGATCTGCCGGCGGCAGTCTTCCAGTGCGCCGGGCAGCGAAGATTGCGCCAGGGCAATGC  >  NC_002947/3489845‑3490130

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