Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,549,921 A→T 6.2% H418L (CAC→CTC) ‡ PP_3133 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,549,9210AT 6.2% 193.4 / 2.8 64H418L (CAC→CTC) ‡PP_3133oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (2/2);  total (32/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 = 9.81e-01

GCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTG  >  NC_002947/3549775‑3550066
                                                                                                                                                  |                                                                                                                                                 
gCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATcaca                                                                                                                                                 >  6:234362/1‑149 (MQ=255)
 cTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTAAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATcaca                                                                                                                                                 <  7:169289/148‑1 (MQ=255)
  ttAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATg                                                                                                                                               >  5:200164/1‑149 (MQ=255)
    aGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATggg                                                                                                                                             >  6:269774/1‑149 (MQ=255)
    aGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATg                                                                                                                                               <  5:183844/147‑1 (MQ=255)
        ccGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcatg                                                                                                                                         <  6:157240/149‑3 (MQ=255)
        ccGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCCTCAAGCTCAACGATTTCGTTTACGACACCACGATACCGCTCAGAATGGCACCCCATGCCCATCACATGGGaacc                                                                                                                                         >  8:145878/1‑149 (MQ=255)
                cGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCAt                                                                                                                                 <  6:200164/149‑1 (MQ=255)
                    acacCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCAt                                                                                                                                 <  5:58917/145‑1 (MQ=255)
                    acacCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCAt                                                                                                                                 >  6:58917/1‑145 (MQ=255)
                    acacCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCAt                                                                                                                                 >  2:109939/1‑145 (MQ=255)
                    acacCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAACCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCAt                                                                                                                                 <  1:109939/145‑1 (MQ=255)
                      acCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGgccgc                                                                                                                           <  5:269774/149‑1 (MQ=255)
                          tCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcat                                                                                                                                          <  3:254896/130‑2 (MQ=255)
                          tCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATGTGATGGGcat                                                                                                                                          >  4:254896/1‑129 (MQ=255)
                            aaGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTcgcc                                                                                                                     >  8:22082/1‑149 (MQ=255)
                            aaGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTcgcc                                                                                                                     >  1:196593/1‑149 (MQ=255)
                            aaGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTcgc                                                                                                                      <  8:262354/148‑1 (MQ=255)
                              gTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTcgccgc                                                                                                                   >  2:94792/1‑149 (MQ=255)
                                 cATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGgcc                                                                                                                             <  7:90912/136‑1 (MQ=255)
                                 cATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGgcc                                                                                                                             >  8:90912/1‑136 (MQ=255)
                                   ttGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTAAGACACATCGATACAGCTGAGAATGGCACCCCATTCCAATCACATGGTCACAAAGCGAATGGCACCCTCGCCGCAGTTc                                                                                                              >  5:23455/1‑149 (MQ=255)
                                      gCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGgt                                                                                                           <  7:145878/149‑1 (MQ=255)
                                       cAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGgtg                                                                                                          <  3:104576/149‑1 (MQ=255)
                                         gCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGtgtg                                                                                                        <  6:218553/149‑1 (MQ=255)
                                         gCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGtgtg                                                                                                        <  7:284793/149‑1 (MQ=255)
                                             aCTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTc                                                                                                     >  6:57575/1‑148 (MQ=255)
                                                 gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGAt                                                                                                  <  5:57575/147‑1 (MQ=255)
                                                 gCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                >  6:259084/1‑149 (MQ=255)
                                                          ttCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATca                                                                                                                                                   >  5:113512/1‑89 (MQ=255)
                                                          ttCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATca                                                                                                                                                   <  6:113512/89‑1 (MQ=255)
                                                          ttCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGc                                                                                         >  6:2280/1‑147 (MQ=255)
                                                          ttCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCaa                                                                                       <  5:234362/149‑1 (MQ=255)
                                                           tCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTTAGAATGGCACCCCATGCCCATAACATGGGCACCACGCGCATGGCCGCCTCTCNGCAGTTCGGTGTGTTCGATGCCAATTGCaa                                                                                       <  1:31803/148‑1 (MQ=255)
                                                            cGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAg                                                                                      <  7:22082/148‑1 (MQ=255)
                                                            cGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAg                                                                                      <  1:202539/148‑1 (MQ=255)
                                                             gCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGgt                                                                                    >  6:152615/1‑149 (MQ=255)
                                                                   aTGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                >  5:83437/1‑131 (MQ=255)
                                                                   aTGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                <  6:83437/131‑1 (MQ=255)
                                                                   aTGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCGAGGTGTTCgg                                                                              <  2:196593/149‑1 (MQ=255)
                                                                       aCCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTAcc                                                                          >  2:69392/1‑149 (MQ=255)
                                                                            gggTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGaaaa                                                                     >  6:103375/1‑149 (MQ=255)
                                                                             ggTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAAc                                                                    <  5:103375/149‑1 (MQ=255)
                                                                                  gTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTa                                                               <  5:2280/149‑1 (MQ=255)
                                                                                   tCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGc                                                                                                >  6:330003/1‑115 (MQ=255)
                                                                                         tCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATAACATGGGCACCACGCGCATGGCCGCCTCGACGCAGTTCGGTTTGTTCGATGCCAATTGCAAGGTGTACTGTACCTAAAACATTTATGTCGcc                                                        >  4:176504/1‑149 (MQ=255)
                                                                                           aaCGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTa                                                                            <  8:51568/127‑1 (MQ=255)
                                                                                           aaCGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCgg                                                      >  1:342360/1‑149 (MQ=255)
                                                                                               aTTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATg                                                                                                                                >  2:102297/1‑71 (MQ=255)
                                                                                               aTTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATg                                                                                                                                <  1:102297/71‑1 (MQ=255)
                                                                                                 ttCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGcc                                                        >  4:272457/1‑141 (MQ=255)
                                                                                                 ttCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGcc                                                        <  3:272457/141‑1 (MQ=255)
                                                                                                     tttACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAtt                                            <  5:259084/149‑1 (MQ=255)
                                                                                                         cGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGGCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCAttttt                                         <  6:23455/148‑1 (MQ=255)
                                                                                                               cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCaa                                                                                             >  5:47244/1‑90 (MQ=255)
                                                                                                               cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCaa                                                                                             <  6:47244/90‑1 (MQ=255)
                                                                                                               cTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAgg                                  >  1:254945/1‑149 (MQ=255)
                                                                                                               cTCGATAACGCTGAGAATGGCACCCCATGCACATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGATCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGACGCCCGCAGCAGCATTATTGCCAAgg                                  >  3:36810/1‑149 (MQ=255)
                                                                                                                tCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAggg                                 >  8:106567/1‑149 (MQ=255)
                                                                                                                       cGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGcc                          >  7:71676/1‑149 (MQ=255)
                                                                                                                                  gCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGAc               <  5:152615/149‑1 (MQ=255)
                                                                                                                                         gtGCCCATGTGATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCg        >  2:40347/2‑149 (MQ=255)
                                                                                                                                          tGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGc       <  1:94792/149‑1 (MQ=255)
                                                                                                                                             ccATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGtt                                                                                 >  8:138368/1‑72 (MQ=255)
                                                                                                                                             ccATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGt    >  8:32118/1‑149 (MQ=255)
                                                                                                                                             ccATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGt    <  8:71676/149‑1 (MQ=255)
                                                                                                                                              cATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGtt   >  5:239810/1‑149 (MQ=255)
                                                                                                                                               aTCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTg  <  2:254945/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTTAGCGCCGAAGACCGGCACACCATCAAGTGCATTGGCAGCGAACTGGCCAAGCAGTTCGCCGACATGGACCTGGGGTTCGTCAAGCTCAACGATTTCGTTTACGACACCTCGATACCGCTGAGAATGGCACCCCATGCCCATCACATGGGCACCACGCGCATGGCCGCCTCGCCGCAGTTCGGTGTGGTCGATGCCAATTGCAAGGTGTTCGGTACCGAAAACCTGTATGTCGCCGGCAGCAGCATTTTTGCCAAGGGCGGCGCCTCCAACCCGACCATGCCGCTGTTG  >  NC_002947/3549775‑3550066

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