Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,869,622 G→C 6.8% A44P (GCC→CCC)  PP_5144 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,869,6220GC 6.8% 188.2 / 4.5 59A44P (GCC→CCC) PP_5144hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/22);  new base C (2/2);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCAGGG  >  NC_002947/5869481‑5869764
                                                                                                                                             |                                                                                                                                               
gagGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGgcg                                                                                                                                          <  8:238710/149‑1 (MQ=255)
   gCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGa                                                                                                                                       >  3:279375/1‑149 (MQ=255)
     ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                                     >  3:114500/1‑149 (MQ=255)
         ccTATGTGCCGGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGa                                                                                                                                 <  3:35975/149‑1 (MQ=255)
            aTGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGagg                                                                                                                               >  6:9477/1‑148 (MQ=255)
                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGaggcaggc                                                                                                                          <  5:9477/148‑1 (MQ=255)
                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                         >  5:110937/1‑149 (MQ=255)
                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                         <  4:279375/149‑1 (MQ=255)
                   tGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCgg                                                                                                                        >  5:186746/1‑148 (MQ=255)
                                cTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                          >  8:214417/1‑149 (MQ=255)
                                   gcggcCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAg                                                                                                       >  2:213170/1‑149 (MQ=255)
                                         ggggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCcac                                                                                                 <  3:194427/149‑1 (MQ=255)
                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggc                                                                                                                <  1:146596/132‑1 (MQ=255)
                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggc                                                                                                                >  2:146596/1‑132 (MQ=255)
                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                          <  1:336541/138‑1 (MQ=255)
                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                          >  2:336541/1‑138 (MQ=255)
                                              aCGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATTCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAgggg                                                                                            >  1:29029/1‑149 (MQ=255)
                                                gCCCTGCCGCTGATCGTGGCGGCCAACCGCGATGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCCGGAGGTACCTGGCTGGGCGTAGGCCCACAGGGgcg                                                                                          >  4:65760/1‑149 (MQ=255)
                                                gCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGgcg                                                                                          >  2:87501/1‑149 (MQ=255)
                                                  ccggCCGCTGATCGT‑GCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCtt                                                                                       <  4:124974/146‑1 (MQ=255)
                                                  ccTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCt                                                                                        >  4:281268/1‑149 (MQ=255)
                                                    tGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCT‑GCTGGGCGTAGGCCCACAGGGGCGCTTTg                                                                                     >  4:259521/1‑149 (MQ=255)
                                                            aTCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCt                                                                              <  2:238575/149‑1 (MQ=255)
                                                             tCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGTGGCGCTTTGCCGCGCt                                                                              <  6:110937/148‑1 (MQ=255)
                                                                 ggcggcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCa                                                                         <  1:87501/149‑1 (MQ=255)
                                                                    ggcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACa                                                                      <  6:186746/149‑1 (MQ=255)
                                                                    ggcCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACa                                                                      <  1:245731/149‑1 (MQ=255)
                                                                        aaCCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCAGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCAACAGGGGCGCTTTGCCGCGCTGACCAACATCcg                                                                  >  1:207203/1‑149 (MQ=255)
                                                                            gcgACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGAc                                                              >  2:82554/1‑149 (MQ=255)
                                                                            gcgACGAATTCTAAGCCCGCCCCACTCAGCCACTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTCCCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCAGCGCTGACCAACAGCcgcgcc                                                              >  6:324584/1‑147 (MQ=255)
                                                                              gacgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGAccc                                                            <  4:114500/149‑1 (MQ=255)
                                                                               acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGcc                                                                                   >  1:81936/1‑125 (MQ=255)
                                                                               acgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGcc                                                                                   <  2:81936/125‑1 (MQ=255)
                                                                                cgaATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑gg                                                         >  2:243245/1‑149 (MQ=255)
                                                                                    ttCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGc                                                       >  7:43303/1‑147 (MQ=255)
                                                                                    ttCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGc                                                       >  7:38224/1‑147 (MQ=255)
                                                                                     tCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAg                                                    >  2:98700/1‑149 (MQ=255)
                                                                                     tCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCCGGGGTGCATGCAGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCAAg                                                    >  2:94399/1‑149 (MQ=255)
                                                                                          gcccgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGgcg                                                 >  1:196221/1‑147 (MQ=255)
                                                                                             cgcccCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAggg                                            <  3:66421/149‑1 (MQ=255)
                                                                                                   aCTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCctcg                                      >  5:267167/1‑149 (MQ=255)
                                                                                                       aGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCctcgctcg                                  >  2:314978/1‑149 (MQ=255)
                                                                                                        gCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCtcgctcgc                                 >  2:155328/1‑149 (MQ=255)
                                                                                                          cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgg                               >  8:34041/1‑149 (MQ=255)
                                                                                                          cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgg                               <  7:214417/149‑1 (MQ=255)
                                                                                                          cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgg                               >  2:236435/1‑149 (MQ=255)
                                                                                                            cTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCgggg                             >  8:223377/1‑149 (MQ=255)
                                                                                                               gcggcCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCAGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTg                           <  1:213170/149‑1 (MQ=255)
                                                                                                                       ggaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggt                    <  1:314978/147‑1 (MQ=255)
                                                                                                                       ggaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCT‑GCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtg                   <  3:259521/147‑1 (MQ=255)
                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                            <  6:306120/134‑1 (MQ=255)
                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                            <  6:311168/134‑1 (MQ=255)
                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                            >  5:311168/4‑137 (MQ=255)
                                                                                                                         ggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                            >  5:306120/4‑137 (MQ=255)
                                                                                                                           gaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGc              <  1:243245/149‑1 (MQ=255)
                                                                                                                                   tGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTGCCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg       >  6:127827/1‑148 (MQ=255)
                                                                                                                                    ggggTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCAACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg       >  5:336976/1‑147 (MQ=255)
                                                                                                                                     gggTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTGCCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCa     <  5:127827/148‑1 (MQ=255)
                                                                                                                                       gTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCAggg  >  2:58453/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCC‑GGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTGCAGGG  >  NC_002947/5869481‑5869764

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