Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,524,526 T→C 6.2% intergenic (+105/‑88) idh → / → PP_4013 isocitrate dehydrogenase/MutT/nudix family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,524,5260TC 6.2% 293.4 / 4.9 99intergenic (+105/‑88)idh/PP_4013isocitrate dehydrogenase/MutT/nudix family protein
Reads supporting (aligned to +/- strand):  ref base T (53/39);  new base C (3/3);  total (57/42)
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.37e-01

AAGCCAGACCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTCCTCT  >  NC_002947/4524385‑4524674
                                                                                                                                             |                                                                                                                                                    
tgaagacgtgctggttggctttgaattcctcgacgaagaggaacttgccttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTGAAcccgctc                                                                                                                                              >  1:128616‑M1/128‑148 (MQ=255)
   ccAGACCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAGCAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATg                                                                                                                                           >  2:417223/1‑150 (MQ=255)
   ccAGACCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATGTCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATg                                                                                                                                           >  1:103555/1‑150 (MQ=255)
       aCCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTc                                                                                                                                       >  2:824523/1‑150 (MQ=255)
       aCCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTc                                                                                                                                       >  2:1193063/1‑150 (MQ=255)
       aCCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTc                                                                                                                                       >  1:24180/1‑150 (MQ=255)
        cccTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCa                                                                                                                                      >  2:955155/1‑150 (MQ=255)
        cccTGAACAGCGCCATCGCTGCCCTGAAAGGTCCCCGGTTAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCa                                                                                                                                      <  2:824797/150‑1 (MQ=255)
           tGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCAttt                                                                                                                                   <  2:429690/150‑1 (MQ=255)
             aaCAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTTTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGc                                                                                                                                 <  1:275188/150‑1 (MQ=255)
              aCAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGcc                                                                                                                                >  2:1397478/1‑150 (MQ=255)
                aGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCACGGCCGCGAACCGGGGTTTGTGCTTGCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTc                                                                                                                              >  1:1243792/1‑150 (MQ=255)
                      aTCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTACATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATc                                                                                                                        <  2:1313155/150‑1 (MQ=255)
                          cTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCAt                                                                                                                    <  2:999680/150‑1 (MQ=255)
                           tGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCAtt                                                                                                                   >  2:1407840/1‑150 (MQ=255)
                            ctcgacgaagaggaacttgccttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTc                                                                                                                  >  2:768138‑M1/100‑150 (MQ=255)
                               cTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTcgcg                                                                                                               >  1:1376298/1‑150 (MQ=255)
                                tGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              <  1:1193063/150‑1 (MQ=255)
                                tGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              <  2:24180/150‑1 (MQ=255)
                                tGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              >  2:705205/1‑150 (MQ=255)
                                acgaagaggaacttgccttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              >  1:1164572‑M1/96‑150 (MQ=255)
                                 gTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCggg                                                                                                             <  2:766857/150‑1 (MQ=255)
                                      gTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAAc                                                                                                        <  1:941976/150‑1 (MQ=255)
                                      gTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAAc                                                                                                        <  1:1397478/150‑1 (MQ=255)
                                       ggaacttgccttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTCACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              <  1:190877‑M1/55‑1 (MQ=255)
                                       ggaacttgccttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTCACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              >  2:190877‑M1/89‑143 (MQ=255)
                                          gCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGc                                                                                                    >  2:778352/1‑150 (MQ=255)
                                           cTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCt                                                                                                   >  1:591670/1‑150 (MQ=255)
                                           cTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCt                                                                                                   >  2:701700/1‑150 (MQ=255)
                                               ccttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTa                                                                                               <  2:1117803‑M1/70‑1 (MQ=255)
                                               aaGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTa                                                                                               >  2:166102/1‑150 (MQ=255)
                                                aGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCTGCCGTTGGAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTAc                                                                                              >  2:760184/1‑150 (MQ=255)
                                                 ttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTCACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              >  2:880026‑M1/79‑133 (MQ=255)
                                                 ttcgtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTCACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              <  1:880026‑M1/55‑1 (MQ=255)
                                                    gtgctcgacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAggg                                                                                          >  1:1117803‑M1/76‑150 (MQ=255)
                                                    aCAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAggg                                                                                          >  1:1276120/1‑150 (MQ=255)
                                                    aCAGCACAAACCCCGGCCACGCACCGGGGTGTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAggg                                                                                          >  2:1155113/1‑150 (MQ=255)
                                                         acaAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCAt                                                                                     >  2:65574/1‑150 (MQ=255)
                                                         acaAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCAt                                                                                     <  1:701700/150‑1 (MQ=255)
                                                         acaAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCAt                                                                                     >  1:1058974/1‑150 (MQ=255)
                                                         acaAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTAATCGGCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCag                                                                                     >  1:1366061/1‑149 (MQ=255)
                                                          gacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTcgc                                                                                                                <  1:1007434‑M1/53‑1 (MQ=255)
                                                          gacgatggtggcgacggtgatgtggggttgccaggtcatgaacagctccttgaaggcgtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTcgc                                                                                                                >  2:1007434‑M1/70‑122 (MQ=255)
                                                           aaaCCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGc                                                                                   >  1:815344/1‑150 (MQ=255)
                                                               cccGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGc                                                                               >  2:1256245/1‑150 (MQ=255)
                                                               cccGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGc                                                                               <  2:1243792/150‑1 (MQ=255)
                                                                ccGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGcc                                                                              >  2:691901/1‑150 (MQ=255)
                                                                  ggCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATGTCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCtt                                                                            <  2:103555/150‑1 (MQ=255)
                                                                     gcgtggccGGGGTTTGTGTTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCTGGTGCTCCCGCTCCTACAGGGGTCATGCACGCCTTCaa                                                                         <  1:216182/144‑1 (MQ=38)
                                                                     cACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCaa                                                                         <  1:778352/150‑1 (MQ=255)
                                                                     cACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCaa                                                                         <  1:963833/150‑1 (MQ=255)
                                                                        tggccGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGa                                                                      <  2:153147/147‑1 (MQ=255)
                                                                        gCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGa                                                                      <  2:47472/150‑1 (MQ=255)
                                                                         cACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAg                                                                     <  1:417223/150‑1 (MQ=255)
                                                                           ccGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCt                                                                   <  1:1407840/150‑1 (MQ=255)
                                                                            cGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTg                                                                  <  1:166102/150‑1 (MQ=255)
                                                                               ggTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTc                                                               <  1:955155/150‑1 (MQ=255)
                                                                                 tttGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCAt                                                             <  2:1376298/150‑1 (MQ=255)
                                                                                 tttGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACACGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCAt                                                             <  2:1114913/150‑1 (MQ=255)
                                                                                   tgtgCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGa                                                           <  1:824523/150‑1 (MQ=255)
                                                                                     tgCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCt                                                                   >  1:564977/1‑140 (MQ=255)
                                                                                     tgCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCt                                                                   <  2:564977/140‑1 (MQ=255)
                                                                                      gCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCt                                                        >  2:914210/1‑150 (MQ=255)
                                                                                       cTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTg                                                       <  1:65574/150‑1 (MQ=255)
                                                                                            tGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAAc                                                  <  2:591670/150‑1 (MQ=255)
                                                                                            tGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAAc                                                  >  1:1188789/1‑150 (MQ=255)
                                                                                              gATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAg                                                                     >  2:921944/1‑129 (MQ=255)
                                                                                              gATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAg                                                                     <  1:921944/129‑1 (MQ=255)
                                                                                                tAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGACTACGTCATTTGCCTCCCGATCACATTCGGGGGTGAGCCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATCACCTGGCAACCCca                                              >  1:883099/1‑150 (MQ=38)
                                                                                                  aaTGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCa                                                                                                                             <  2:732053/69‑1 (MQ=255)
                                                                                                  aaTGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCa                                                                                                                             >  1:732053/1‑69 (MQ=255)
                                                                                                          cttgaaggcgtgcatgaccctTGTAGGAGCGGGTTTACCCGCGAGTGCGTCATTTGCCTCACCATCGCGCGTTCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCGACATCACCGTc                                    <  1:811073‑M1/129‑1 (MQ=255)
                                                                                                          cTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTc                                    <  1:705205/150‑1 (MQ=255)
                                                                                                           ttgaaggcgtgcatgaccccTGTAGGAGCGGGTTCACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              <  2:311103‑M1/55‑1 (MQ=255)
                                                                                                           ttgaaggcgtgcatgaccccTGTAGGAGCGGGTTCACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCgg                                                                                                              >  1:311103‑M1/21‑75 (MQ=255)
                                                                                                              cggcCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGcca                                <  1:691901/150‑1 (MQ=255)
                                                                                                                 ccTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAg                                                                     <  2:458686/110‑1 (MQ=255)
                                                                                                                 ccTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAg                                                                     >  1:458686/1‑110 (MQ=255)
                                                                                                                   gtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGc                                                                                                                                          <  1:693318‑M1/27‑1 (MQ=255)
                                                                                                                   gtgcatgaccccTGTAGGAGCGGGTTTACCCGCGAATGc                                                                                                                                          >  2:693318‑M1/13‑39 (MQ=255)
                                                                                                                     cATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCAtcgt                         >  1:1216003/1‑150 (MQ=255)
                                                                                                                       tCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCAtcgtcg                       >  1:821282/1‑150 (MQ=255)
                                                                                                                        cGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGa                      >  1:185914/1‑150 (MQ=255)
                                                                                                                          gCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGc                    >  2:1191589/1‑150 (MQ=255)
                                                                                                                            cTTTGTAGGAGCGGGTTTACCCACGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCAc                  >  1:485678/1‑150 (MQ=255)
                                                                                                                             tttGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACg                 >  2:720955/1‑150 (MQ=255)
                                                                                                                              ttGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGACTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGa                >  2:114684/1‑150 (MQ=255)
                                                                                                                                GTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAg              >  1:1256245‑M1/1‑150 (MQ=255)
                                                                                                                                GTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAg              >  2:52858‑M1/1‑150 (MQ=255)
                                                                                                                                  AGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCCACGTCGAGCACGAAGGc            >  1:760184‑M1/1‑150 (MQ=255)
                                                                                                                                     AGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAg         >  2:821282‑M1/1‑150 (MQ=255)
                                                                                                                                      GCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCACACATCACCGTCGCCACCATCGTCGAGc                    >  1:381280‑M1/1‑138 (MQ=255)
                                                                                                                                      GCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCACACATCACCGTCGCCACCATCGTCGAGc                    <  2:381280‑M1/138‑1 (MQ=255)
                                                                                                                                      GCGCGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGt        >  2:883099‑M1/1‑150 (MQ=255)
                                                                                                                                         GGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGACATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTcc     <  2:479705‑M1/150‑1 (MQ=255)
                                                                                                                                          GTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTCct    >  1:914210‑M1/1‑150 (MQ=255)
                                                                                                                                          GTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTCct    >  1:720955‑M1/1‑150 (MQ=255)
                                                                                                                                           TTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTCctc   >  2:1188789‑M1/1‑150 (MQ=255)
                                                                                                                                            TTACCCGCGAATGCGTTATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTCctct  >  2:811073‑M1/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AAGCCAGACCCTGAACAGCGCCATCGCTGCCCTGTAAGGTTCGCTTGAAGTAACAGCACAAACCCCGGCCACGCACCGGGGTTTGTGCTTTCTGGATAAATGATATCTGGCGGCCTGCATCGGCCTTTGTAGGAGCGGGTTTACCCGCGAATGCGTCATTTGCCTCACCATCGCATTCGCGGGTGAACCCGCTCCTACAGGGGTCATGCACGCCTTCAAGGAGCTGTTCATGACCTGGCAACCCCACATCACCGTCGCCACCATCGTCGAGCACGAAGGCAAGTTCCTCT  >  NC_002947/4524385‑4524674

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: