Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,222,096 A→G 5.9% Y200H (TAC→CAC)  PP_2823 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,222,0960AG 5.9% 199.1 / 2.4 68Y200H (TAC→CAC) PP_2823methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base G (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CTGCACGCGGCGGGTCAGGTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATGGGGGCGT  >  NC_002947/3221953‑3222238
                                                                                                                                               |                                                                                                                                              
cTGCACGCGGCGGGTCAGGTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGa                                                                                                                                           >  1:51149/1‑149 (MQ=255)
      gcggcgGGTCAGGTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCTCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGAAGGTAGCCGAcgcgcg                                                                                                                                     >  6:231182/1‑149 (MQ=255)
      gcggcgGGTCAGGTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGAcgcgcg                                                                                                                                     >  7:159445/1‑149 (MQ=255)
      gcggcgGGTCAGGTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGAcgcgcg                                                                                                                                     <  8:288210/149‑1 (MQ=255)
                 ggTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCg                                                                                                                          >  5:301000/1‑149 (MQ=255)
                   tCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACTAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGtt                                                                                                                        <  1:325595/149‑1 (MQ=255)
                   tCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGtt                                                                                                                        <  6:58430/149‑1 (MQ=255)
                   tCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGtt                                                                                                                        <  6:301802/149‑1 (MQ=255)
                   tCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGtt                                                                                                                        <  6:296925/149‑1 (MQ=255)
                     gCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGt                                                                                                                      >  5:179114/1‑149 (MQ=255)
                      ccATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGtt                                                                                                                        >  1:35402/1‑146 (MQ=255)
                      ccATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGtt                                                                                                                        <  2:35402/146‑1 (MQ=255)
                        aTCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTGGCCGACGCGCgtggtg                                                                                                                                >  6:232012/1‑136 (MQ=255)
                        aTCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTGGCCGACGCGCgtggtg                                                                                                                                <  5:232012/136‑1 (MQ=255)
                          cgccgcTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCg                                                                                                                          >  8:301117/1‑140 (MQ=255)
                          cgccgcTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCg                                                                                                                          <  7:301117/140‑1 (MQ=255)
                           gccgcTGGCAATGGTCTTGAGCGATCAGATGACGCCGCGGATATTGCTCATGACCTGACGACTCACCCAGGCCGCGCTGGCGCTGAACACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCTTGGTGGTGGCGTTGTCGgcgc                                                                                                                <  5:231182/149‑1 (MQ=255)
                              gcTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGa                                                                                                             <  4:276085/149‑1 (MQ=255)
                                        gTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTc                                                                                                     >  4:184109/1‑147 (MQ=255)
                                          cTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAgc                                                                                                 >  4:224559/1‑149 (MQ=255)
                                               gctaTCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGt                                                                                                      >  4:306106/4‑139 (MQ=255)
                                               gCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGt                                                                                                      <  3:306106/139‑1 (MQ=255)
                                                cGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGctct                                                                                           <  4:106210/149‑1 (MQ=255)
                                                          aCGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGc                                                                                 <  6:179114/149‑1 (MQ=255)
                                                              cgcgGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCgg                                                                                  >  7:28527/1‑144 (MQ=255)
                                                              cgcgGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCgg                                                                                  <  8:28527/144‑1 (MQ=255)
                                                                  gACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGAcgcgc                                                                                                                                      >  4:233481/1‑88 (MQ=255)
                                                                     aTTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGc                                                                      >  1:6509/1‑149 (MQ=255)
                                                                         cTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGc                                                                   <  7:70708/148‑1 (MQ=255)
                                                                            aTGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGaaa                                                               <  6:301000/149‑1 (MQ=255)
                                                                                 cTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTgg                                                                                                        <  7:159726/103‑1 (MQ=255)
                                                                                 cTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTgg                                                                                                        >  8:159726/1‑103 (MQ=255)
                                                                                 cTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAg                                                                             >  6:334056/1‑130 (MQ=255)
                                                                                 cTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAg                                                                             <  5:334056/130‑1 (MQ=255)
                                                                                  tGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGc                                                                   <  5:252404/139‑1 (MQ=255)
                                                                                  tGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGc                                                                   >  6:252404/1‑139 (MQ=255)
                                                                                     cGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGcc                                                      <  3:224559/149‑1 (MQ=255)
                                                                                         tCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTGGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAgc                                                  >  8:212095/1‑149 (MQ=255)
                                                                                               aGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTgg                                            <  2:6509/149‑1 (MQ=255)
                                                                                                ggCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTggg                                           <  8:159445/149‑1 (MQ=255)
                                                                                                    gcgctggcgctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGtt                                       <  4:309088/149‑1 (MQ=255)
                                                                                                    gcgctggcgctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCACTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGtt                                       >  7:50596/1‑149 (MQ=255)
                                                                                                       ctggcgctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTGGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCtt                                    <  7:212095/149‑1 (MQ=255)
                                                                                                         ggcgctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGa                                  <  3:184109/149‑1 (MQ=255)
                                                                                                          gcgctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGa                                  <  2:51149/148‑1 (MQ=255)
                                                                                                          gcgctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAg                                 >  4:129954/1‑149 (MQ=255)
                                                                                                            gctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCtt                                    >  4:68780/1‑144 (MQ=255)
                                                                                                            gctgACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCtt                                    <  3:68780/144‑1 (MQ=255)
                                                                                                                accaccAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCt                           >  5:205617/1‑149 (MQ=255)
                                                                                                                 cccccaGCCGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTg                          <  2:308936/146‑1 (MQ=255)
                                                                                                                       gcagcagGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGa                                  >  3:122763/1‑135 (MQ=255)
                                                                                                                       gcagcagGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGa                                  <  4:122763/135‑1 (MQ=255)
                                                                                                                         agcagGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCACTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCg                  <  8:50596/149‑1 (MQ=255)
                                                                                                                              gCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATg                    <  6:347444/142‑1 (MQ=255)
                                                                                                                              gCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATg                    >  5:347444/1‑142 (MQ=255)
                                                                                                                              gCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCa              >  6:57550/1‑148 (MQ=255)
                                                                                                                              gCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCa              >  8:12330/1‑148 (MQ=255)
                                                                                                                              gCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGAGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCa              >  5:257885/1‑148 (MQ=255)
                                                                                                                              gCAAAAACTGAGCAGGTAGCCGACGCGCGTGGAGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGt                                                                         >  4:104452/1‑89 (MQ=255)
                                                                                                                                aaaaaCTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCaaaa           >  7:253346/1‑149 (MQ=255)
                                                                                                                                   aaCTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCt                                     >  4:80776/1‑120 (MQ=255)
                                                                                                                                   aaCTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCt                                     <  3:80776/120‑1 (MQ=255)
                                                                                                                                    aCTGAGCAGGTAGCCGACTCTCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATggg       <  7:145039/149‑1 (MQ=255)
                                                                                                                                     cTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATgggg      <  6:257885/149‑1 (MQ=255)
                                                                                                                                        aGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATggg       <  7:247197/145‑1 (MQ=255)
                                                                                                                                        aGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATggg       >  8:247197/1‑145 (MQ=255)
                                                                                                                                         gCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATGGGGGCGt  <  6:205617/149‑1 (MQ=255)
                                                                                                                                         gCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATGGGGGCGt  <  8:253346/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CTGCACGCGGCGGGTCAGGTCGCCATCGCCGCTGGCAATGGTCTTGAGCGATCCGATGACGCCGCGGACATTGCTCATGACCTGACGAATCACCCAGGCCGCGCTGGCGCTGACCACCAGCAGCAGGCAAAAACTGAGCAGGTAGCCGACGCGCGTGGTGGTGGCGTTGTCGGCGCGGACTTGGGTCAAGCTGCTCTGGAAGTCGGCGTAGGCGTGGCTGCGAAACTGCTCGCCGAGCGCCTGGGCGTTCTTGAGGTCGCTGGCCATGCGATCAAAATGGGGGCGT  >  NC_002947/3221953‑3222238

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