Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,813 T→A 6.3% F145I (TTC→ATC)  PP_0891 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,8130TA 6.3% 237.7 / 4.5 79F145I (TTC→ATC) PP_0891hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base A (2/3);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

TGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCT  >  NC_002947/1033667‑1033960
                                                                                                                                                  |                                                                                                                                                   
tGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCACCTGGAAGCGCCGCGCATc                                                                                                                                                   >  4:124036/1‑149 (MQ=255)
 gAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTc                                                                                                                                                   >  5:291652/1‑148 (MQ=255)
   aGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCa                                                                                                                                                 <  5:245259/148‑1 (MQ=255)
      cacaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCTAggtg                                                                                                                                             <  5:108169/149‑1 (MQ=255)
       acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCACCTGGAAGCGCCGCGCATCCAggcgg                                                                                                                                            <  3:124036/149‑1 (MQ=255)
           aTCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAc                                                                                                                                        <  3:93876/149‑1 (MQ=255)
                 cTGCATCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGa                                                                                                                                  <  4:207876/149‑1 (MQ=255)
                     aGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCct                                                                                                                              >  3:254083/1‑149 (MQ=255)
                      gCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGTGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCctc                                                                                                                             <  8:199393/149‑1 (MQ=255)
                       ccGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCctct                                                                                                                            >  7:185490/1‑149 (MQ=255)
                             aaTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGtgct                                                                                                                      <  4:254083/149‑1 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtggt                                                                                                                                           <  3:271716/127‑1 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtggt                                                                                                                                           >  4:271716/1‑127 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtg                                                                                                                                             >  5:13600/1‑125 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAccc                                                                                                                                      >  4:242930/1‑132 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAccc                                                                                                                                      <  3:242930/132‑1 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGtgctg                                                                                                                     >  6:8875/1‑149 (MQ=255)
                              aTGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGAGGCGCTTCCAggtg                                                                                                                                             <  6:13600/125‑1 (MQ=255)
                                 aCCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTgcg                                                                                                                   <  8:185490/148‑1 (MQ=255)
                                   cTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTgcgcgc                                                                                                                <  2:115828/149‑1 (MQ=255)
                                                 acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCTCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCt                                                                                                  <  3:54873/149‑1 (MQ=255)
                                                   aGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCAggcgg                                                                                                                                            >  6:163820/1‑105 (MQ=255)
                                                   aGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGAAGCGCCGCGCATCCAggcgg                                                                                                                                            <  5:163820/105‑1 (MQ=255)
                                                     aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                     <  5:61457/142‑1 (MQ=255)
                                                     aTCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTg                                                                                                     >  6:61457/1‑142 (MQ=255)
                                                         aaTTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCg                                                                                          <  2:108792/149‑1 (MQ=255)
                                                          aTTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgg                                                                                         >  8:279125/1‑149 (MQ=255)
                                                           ttGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCggg                                                                                        >  3:120784/1‑149 (MQ=255)
                                                              cgcCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCt                                                                                     <  4:120784/149‑1 (MQ=255)
                                                                      gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                             >  3:201021/1‑149 (MQ=255)
                                                                      gCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGcc                                                                             <  1:129426/149‑1 (MQ=255)
                                                                          aCGGTCACCCTGGATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCGTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGTGCTGGAAAGCCTgc                                                                          <  6:291652/148‑1 (MQ=255)
                                                                             gTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCa                                                                       >  8:275523/1‑148 (MQ=255)
                                                                              acACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCAAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTAGCCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAg                                                                     >  2:307442/2‑149 (MQ=255)
                                                                               cACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGt                                                                    >  5:159629/1‑149 (MQ=255)
                                                                                aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgg                                                                                         >  6:227765/1‑127 (MQ=255)
                                                                                aCCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCgg                                                                                         <  5:227765/127‑1 (MQ=255)
                                                                                  ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                  >  2:311841/1‑148 (MQ=255)
                                                                                  ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                  <  1:311841/148‑1 (MQ=255)
                                                                                  ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCg                                                                 >  4:41018/1‑149 (MQ=255)
                                                                                   cTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGc                                                                <  7:279125/149‑1 (MQ=255)
                                                                                      cATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCa                                                              <  3:295161/148‑1 (MQ=255)
                                                                                          gCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGtt                                                          >  5:64638/1‑148 (MQ=255)
                                                                                                  cacaGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                                 >  2:159248/1‑149 (MQ=255)
                                                                                                     aGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCaa                                              >  2:333746/1‑149 (MQ=255)
                                                                                                     aGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCaa                                              >  7:166352/1‑149 (MQ=255)
                                                                                                       cTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                      <  2:129628/123‑1 (MQ=255)
                                                                                                       cTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                      >  1:129628/1‑123 (MQ=255)
                                                                                                          cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCt                                                                            >  2:196915/1‑114 (MQ=255)
                                                                                                          cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                      <  4:162623/120‑1 (MQ=255)
                                                                                                          cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                                      >  3:162623/1‑120 (MQ=255)
                                                                                                          cAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATATTGGCTTGCTGCCGGGGCTGGAAAGCCt                                                                            <  1:196915/114‑1 (MQ=255)
                                                                                                           aaTGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTc                                        <  6:159629/149‑1 (MQ=255)
                                                                                                               tacTGTTGCTCGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTa                                    >  3:219542/3‑149 (MQ=255)
                                                                                                                 cTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAgg                                                            >  6:57269/1‑123 (MQ=255)
                                                                                                                 cTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAgg                                                            <  5:57269/123‑1 (MQ=255)
                                                                                                                    ttGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCg                                 <  7:275523/147‑1 (MQ=255)
                                                                                                                    ttGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGt                               >  1:183647/1‑149 (MQ=255)
                                                                                                                      gCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGt                               >  8:182622/1‑147 (MQ=255)
                                                                                                                       cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                  >  2:319362/1‑111 (MQ=255)
                                                                                                                       cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTc                                                                  <  1:319362/111‑1 (MQ=255)
                                                                                                                       cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAg                            >  4:66425/1‑149 (MQ=255)
                                                                                                                       cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTAGGTACCGGTTa                             <  1:159248/148‑1 (MQ=255)
                                                                                                                        tGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGGCGGTACCGGTTAgc                           >  7:296125/1‑149 (MQ=255)
                                                                                                                        tGGCCACCCCCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAgc                           >  3:259078/1‑149 (MQ=255)
                                                                                                                          gCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAgcgc                         >  8:30589/1‑149 (MQ=255)
                                                                                                                          gCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAgcgc                         >  8:74786/1‑149 (MQ=255)
                                                                                                                           ccACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGcgcg                        <  4:201021/149‑1 (MQ=255)
                                                                                                                            cACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCa                                                  >  1:80368/1‑122 (MQ=255)
                                                                                                                            cACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCa                                                  <  2:80368/122‑1 (MQ=255)
                                                                                                                             aCCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGt                      >  2:16536/1‑149 (MQ=255)
                                                                                                                               ccGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGc                    <  3:41018/149‑1 (MQ=255)
                                                                                                                               ccGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGc                    >  4:113003/1‑149 (MQ=255)
                                                                                                                                cGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCt                   >  3:205306/1‑149 (MQ=255)
                                                                                                                                   cTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGAt                <  3:66425/149‑1 (MQ=255)
                                                                                                                                        tGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCa           <  1:333746/149‑1 (MQ=255)
                                                                                                                                        tGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCa           <  6:64638/149‑1 (MQ=255)
                                                                                                                                           gcggcgCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTg        >  5:223862/1‑149 (MQ=255)
                                                                                                                                                 cTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCt  <  1:16536/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATCTTCACTGCGCGCT  >  NC_002947/1033667‑1033960

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: