Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,126,856 G→C 8.4% A190G (GCC→GGC)  PP_2744 ← ribose‑phosphate pyrophosphokinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,126,8560GC 8.4% 213.4 / 6.9 68A190G (GCC→GGC) PP_2744ribose‑phosphate pyrophosphokinase family protein
Reads supporting (aligned to +/- strand):  ref base G (34/28);  new base C (4/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAGGCGTTGTCG  >  NC_002947/3126708‑3127002
                                                                                                                                                       |                                                                                                                                                   
cGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAggg                                                                                                                                                     >  3:44324/1‑149 (MQ=255)
      ccccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGc                                                                                                                                               >  4:136622/1‑149 (MQ=255)
      ccccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCAGCTG‑‑GCGGTGTTTCTCCAGACTGGCCACCGAGCCCGCTCGCGCCAGCAGGTGCTCCAGGCCCTGGc                                                                                                                                               >  4:410220/1‑149 (MQ=255)
                 tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGCCCGGTCGCGCCCGCAGGTGCGCCAGGGCCTGGCGAAAATGCTcc                                                                                                                                   >  2:217579/1‑148 (MQ=255)
                 tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGAACAGGTGTTCGCCGCTAAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCTCGCCAGAAGGTGCGCAAGTGCCTGGCGAAACTGCTCg                                                                                                                                   >  2:349774/1‑149 (MQ=255)
                    ggTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTTCGCCAGGCCCTGGCGCACCTGCTGgcgc                                                                                                                               <  5:81634/149‑1 (MQ=255)
                         tcgAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGttt                                                                                                                           <  3:136622/149‑1 (MQ=255)
                            aaTACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGAc                                                                                                                        <  2:59966/149‑1 (MQ=255)
                            aaTACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGAc                                                                                                                        >  1:54767/1‑149 (MQ=255)
                              tACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGc                                                                                                                      >  3:307383/1‑149 (MQ=255)
                                cGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCa                                                                                                                    >  1:154045/1‑149 (MQ=255)
                                 gATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGAGAGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGc                                                                                                                      >  1:348563/1‑148 (MQ=255)
                                 gATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCAc                                                                                                                   >  4:323299/1‑149 (MQ=255)
                                    gACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCg                                                                                                                 <  6:52936/148‑1 (MQ=255)
                                       ggTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGatat                                                                                                             >  2:274635/1‑149 (MQ=255)
                                          gctgcCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCg                                                                                                          <  4:290790/149‑1 (MQ=255)
                                           ctgcctgcCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCgg                                                                                                         <  4:22838/149‑1 (MQ=255)
                                           ctgcctgcCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCgg                                                                                                         <  4:310198/149‑1 (MQ=255)
                                            tgcctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGTGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGt                                                                                                        >  1:164919/1‑149 (MQ=255)
                                            tgcctgccACGTTACCGACCAGGTGTTAGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTAGGCGCGTTTGACGCCACCGATATCCGGt                                                                                                        >  7:64194/1‑149 (MQ=255)
                                              cctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCtggc                                                                                                                                  <  6:176166/122‑4 (MQ=37)
                                              cctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCGAGCAGTTTCGCCAGGCCCTGGCGCACCTGCtggc                                                                                                                                  >  5:176166/1‑119 (MQ=37)
                                              cctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGa                                                                                                      >  2:96312/1‑149 (MQ=255)
                                              cctgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGa                                                                                                      >  2:304184/1‑149 (MQ=255)
                                                tgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGaca                                                                                                    >  6:430956/1‑149 (MQ=255)
                                                     cGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGAcc                                                                                               <  5:430956/149‑1 (MQ=255)
                                                       ttACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGt                                                                                             <  3:323299/149‑1 (MQ=255)
                                                         aCCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGACAACTGCTCGGCTCGTTTGACGCCACCGATATCCGGTGACACGACCGt                                                                                             <  1:96312/148‑1 (MQ=255)
                                                         aCCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTg                                                                                            <  2:154045/148‑1 (MQ=255)
                                                               cAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCt                                                                                     <  2:348563/149‑1 (MQ=255)
                                                                aGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTg                                                                                    <  4:44324/149‑1 (MQ=255)
                                                                     ttCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCAc                                                                               <  2:141944/149‑1 (MQ=255)
                                                                     ntCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACTGTGACCGCCTGCTCAc                                                                               >  4:191327/2‑149 (MQ=255)
                                                                       cgccgcTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCAcc                                                                              >  6:207747/1‑148 (MQ=255)
                                                                          cgcTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAg                                                                           >  3:423815/1‑148 (MQ=255)
                                                                                ggCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCggg                                                                    <  1:274635/149‑1 (MQ=255)
                                                                                ggCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCggg                                                                    >  7:320151/1‑149 (MQ=255)
                                                                                ggCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCggg                                                                    >  3:162365/1‑149 (MQ=255)
                                                                                  cctgcctgCTG‑‑GCGGTGTTTCTCCCTCATGGCCACCGAGACCGGTCGCGCCATCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCggg                                                                    <  7:149343/147‑1 (MQ=255)
                                                                                     gcctgcTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCgggg                                                                   <  8:256915/145‑1 (MQ=255)
                                                                                     gcctgcTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCgggg                                                                   >  7:256915/1‑145 (MQ=255)
                                                                                      cctgcTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAgg                                                              >  3:49574/1‑149 (MQ=255)
                                                                                       ctgctg‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGc                                                             >  1:417954/1‑149 (MQ=255)
                                                                                        tgctg‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCg                                                            >  3:10316/1‑149 (MQ=255)
                                                                                               gCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTggg                                                       >  6:216656/1‑149 (MQ=255)
                                                                                               gCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTggg                                                       <  1:304184/149‑1 (MQ=255)
                                                                                                   tgtTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGTGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                                                   <  2:164919/149‑1 (MQ=255)
                                                                                                   tgtTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCAAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGaa                                                   <  8:155178/149‑1 (MQ=255)
                                                                                                          ccATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTgg                                            <  4:10316/149‑1 (MQ=255)
                                                                                                           catcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCt                                              <  4:49574/146‑1 (MQ=255)
                                                                                                            atcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGt                                          <  5:216656/149‑1 (MQ=255)
                                                                                                            atcatGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGt                                          >  2:345361/1‑149 (MQ=255)
                                                                                                              catGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCgg                                                                     >  3:154616/1‑120 (MQ=255)
                                                                                                              catGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCgg                                                                     <  4:154616/120‑1 (MQ=255)
                                                                                                              catGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTg                                         >  8:94834/1‑148 (MQ=255)
                                                                                                               atGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCCCCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCa                                       >  7:299706/1‑149 (MQ=255)
                                                                                                                   ccACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGc                                   <  4:423815/149‑1 (MQ=255)
                                                                                                                    cACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGGTGGGCGAATAGCTGGGTGCATTg                                    >  7:345765/1‑147 (MQ=255)
                                                                                                                       cGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGt                               >  2:181836/1‑149 (MQ=255)
                                                                                                                          gaCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTc                             >  8:28718/1‑148 (MQ=255)
                                                                                                                                  gcgcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCgg                     >  3:32598/1‑148 (MQ=255)
                                                                                                                                   cgcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGa                   <  5:278716/149‑1 (MQ=255)
                                                                                                                                    gccgagcagTTTCGCCAGGCCCTGGCGCACCTGCTGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACGGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGAt                  >  1:139866/5‑149 (MQ=38)
                                                                                                                                    gcCAGCAGGTGCGCCAGGGCGTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGAt                  <  1:217579/149‑1 (MQ=255)
                                                                                                                                    gcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCa                           <  4:344930/140‑1 (MQ=255)
                                                                                                                                    gcCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCa                           >  3:344930/1‑140 (MQ=255)
                                                                                                                                           agGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACTGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTa            <  3:191327/149‑1 (MQ=255)
                                                                                                                                             tttcgcCAGGCCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAgg          >  3:429882/4‑149 (MQ=255)
                                                                                                                                                 gcCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAGGCGtt      >  6:270379/1‑149 (MQ=255)
                                                                                                                                                       gCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGTGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAGGCGTTGTCg  >  4:74911/1‑147 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTG‑‑GCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAGGCGTTGTCG  >  NC_002947/3126708‑3127002

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: