Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,745,479 C→A 20.5% intergenic (+39/+32) hutF → / ← PP_5037 formiminoglutamate deiminase/lipocalin family lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,745,4790CA20.5% 204.9 / 50.8 95intergenic (+39/+32)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (17/58);  new base A (9/10);  total (26/69)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

TGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGT  >  NC_002947/5745331‑5745628
                                                                                                                                                      |                                                                                                                                                      
tGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGccc                                                                                                                                                       <  1:3907668/150‑1 (MQ=255)
  cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcgc                                                                                                                                                     <  2:5547446/150‑1 (MQ=255)
  cgcgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcgc                                                                                                                                                     <  2:4826837/150‑1 (MQ=255)
  cgcgACGTGATGGTGCCTGGGCCAGGGGAGGTGCGCCAGGGGCGGCATGCGCAGGACGATGACAGCGGGCGGGCGTTTGCCAAGGTGGTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcgc                                                                                                                                                     <  2:3147691/150‑1 (MQ=255)
    cgATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGa                                                                                                                                                   <  2:3873864/150‑1 (MQ=255)
     gATGTGATGGTGCATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                                  <  1:3331599/150‑1 (MQ=255)
     gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                                  <  1:946519/150‑1 (MQ=255)
     gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                                  <  1:946650/150‑1 (MQ=255)
     gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGaa                                                                                                                                                  <  1:3023962/150‑1 (MQ=39)
     gATGTGATGGTGAATGGGCAATGGGTGGTGCGACAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                                  <  2:41527/150‑1 (MQ=255)
      aTGTTATGGTGAATGGGGTATGGGTGGTGCGCCATTGTGGTCATGCCCAGGAAGATGAAAGAGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGTGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCTCGCGatg                                                                                                                                                 <  2:1897354/150‑3 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                                 <  1:2355199/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                                 <  2:267106/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                                 <  2:4690180/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                                 <  2:1197485/150‑1 (MQ=39)
        gtgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAggg                                                                                                                                               <  2:1029026/150‑1 (MQ=39)
         tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                              <  1:1460843/150‑1 (MQ=255)
          gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                             <  1:806313/150‑1 (MQ=255)
          gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                             <  2:2537380/150‑1 (MQ=255)
          gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                             <  2:2262162/150‑1 (MQ=255)
          gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                             <  1:482778/150‑3 (MQ=38)
          gATGGTGAATGGGCAAGGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                             <  2:1021591/150‑3 (MQ=38)
          gATGGTGAAGGGGCCATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGACAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                             <  2:3337813/150‑1 (MQ=255)
           aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                            <  2:4163191/150‑1 (MQ=255)
           aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                            <  2:4163190/150‑1 (MQ=255)
           aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                            <  2:2639253/150‑1 (MQ=39)
             ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgca                                                                                                                                          <  1:3588279/150‑1 (MQ=255)
              gTGAATGGGCCATGGGTGGTGCGGCAGGGGCGGCATGCCCAGGCCGATGACCGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCTGCTTTTGGGCTGAGCGTTGCGGGGTTGCGTTGGAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcat                                                                                                                                         <  2:4244154/150‑2 (MQ=255)
              gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaa                                                                                                                                         <  1:380543/150‑1 (MQ=38)
              gTGAATGGGCAATGGGTGGAGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa                                                                                                                                         <  1:4334855/150‑1 (MQ=255)
               tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                        <  1:1506068/150‑1 (MQ=255)
               tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                        <  2:4094007/150‑1 (MQ=255)
               tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                        <  2:4094055/150‑1 (MQ=255)
               tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                        <  1:4341615/150‑1 (MQ=255)
                    gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGGAGCCCGTAGCGGGCAAGACCCCGCAGGGTTgaaacgcaac                                                                                                                                  >  2:5117264/1‑150 (MQ=255)
                              tggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGAGCCCAGCGCGGGGCAGCCCCCCCACGGGTTGAACGCCCACCCCaaaaaa                                                                                                                         >  2:4457637/1‑150 (MQ=255)
                                 tgCGCCTGGGGCGGCATCCCCCGGAAGAGGAAAGGGGGGGGGCGTTGGCCATGGTGTTGCGCCAGCTTTTGGGCTGCGCGTCCGGGGGTTGCTTTGCAGCCCTTCGCGTGAACGCCCGCTATGGTTTGCAACGCAACCCCAAAAAACTCa                                                                                                                     <  2:2536136/150‑1 (MQ=255)
                                           gcggcATGCCCCGGATTATGACCGCGGGGGGGGGTTTGCATCGGTTTTGCGCGAGCTTTTGGGGTGAGGTTACCGGGGTTGCGTTGCTGCCCTTGGGGTGCACGCCCGCGACGGTTTGAAACGCAACCCCAAAAAACTCAGTGCAttttt                                                                                                           <  2:3976542/150‑1 (MQ=255)
                                                    ccAGGAAGATGAACGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGTTTGCGTTGCAGCCCTCCGCGTGCACGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                                  <  2:5276798/150‑1 (MQ=255)
                                                         aaGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGTTGCAGCCCTTCGGGTGCAAGCCCGCTAAGGTTTGCAACGCATCCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATc                                                                                             <  2:3705035/150‑1 (MQ=255)
                                                         aaGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGGGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATc                                                                                             <  2:3188860/150‑1 (MQ=255)
                                                             tgaccgcgggcggGCGTTTGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                         <  2:3742233/145‑1 (MQ=255)
                                                                       ggggctTTTGCCATGGTGTTGCGCCAGGTTTTGGGCTGTGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCGGCGACGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGa                                                                               <  2:1563984/149‑1 (MQ=255)
                                                                         ggctTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTGGCGTTGCAGCCCTGCGCTTGCAAGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATc                                                                             <  2:3823927/150‑1 (MQ=255)
                                                                           ggTTTGCCATGGTGTTGCGCATTCTTTTGGGGTGCGCGTACGGGGTTTGCGTTGCAGCCCTTGGGGTGCAAGCCCGCTTTGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                           <  1:3675670/149‑1 (MQ=255)
                                                                            gTTTGCCTCGGTGTGGCGCCAGCTTTTGGGCTGAGCGTACGGGGTTTGCTTTGCAGCCCTTCGCGTGAAAGCCGGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAg                                                                          <  2:68901/150‑1 (MQ=255)
                                                                                  cacggTGTTGCGCCAGCTTTTGGGCGGAGCGTCCCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                    <  2:4843025/147‑1 (MQ=255)
                                                                                      gtgtGGCGCCAGATTTTGGGCTGTGGTTACGGGGGTTGGTTTGCAGCCCTTCGGGTGCAAGCCCGCTATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtgt                                                                <  1:4148304/150‑1 (MQ=255)
                                                                                           gcgcCCGCTTTGGGGCTGCGCGTCCGGGGGTTGGGTTGCAGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc                                                           <  2:112173/150‑1 (MQ=255)
                                                                                               ccgCTTTTGGGCTGAGCGTACCGGGGTGGCGTGGCAGCCCTTCGCGTGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCAGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                       <  2:4431453/148‑1 (MQ=255)
                                                                                               cAGCTTTGGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                       <  2:1846702/150‑1 (MQ=255)
                                                                                                   ttttGGGCTGCGGGTCGGGGGGTGGGGTGGCCGCCCTCCGCGGCAACGCCCGCGAGGGTTTCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgctg                                                   <  1:1677566/150‑1 (MQ=255)
                                                                                                       gggCTGAGCGTACGGGGGTTGCGTTCCAGCCCTTCGCGTGAAAGCCCGCTAAGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                               <  2:2100582/150‑1 (MQ=255)
                                                                                                       gggCTGAGCGTACCGGGGTTGCGTTCCACCCCTACGCGTCCAACCCCGCTAAGGTTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                               <  2:2100490/150‑1 (MQ=255)
                                                                                                                 tACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgca                                                                                                                                          <  2:4551026/50‑1 (MQ=255)
                                                                                                                 tACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgca                                                                                                                                          >  1:4551026/1‑50 (MQ=255)
                                                                                                                    gggggTTGGGTTGCAGCCCTTCGCGGGCATGCCCGCGACGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                  <  2:1012404/149‑1 (MQ=255)
                                                                                                                    gggggTTGGGTTGCAGCCCTGCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTAAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                  <  2:4390373/149‑1 (MQ=255)
                                                                                                                    cGGGTTTGCGTTGCAGCCCTCCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                  <  2:4950710/150‑1 (MQ=255)
                                                                                                                    cGGGGTTGCGGTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                  <  2:3691764/150‑1 (MQ=255)
                                                                                                                      gggTTGCGTGGCACCCGTTCGGGTGCACGCCGGGGTTGGTTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg                                <  1:4403749/150‑1 (MQ=255)
                                                                                                                       ggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                               <  2:4508567/150‑1 (MQ=255)
                                                                                                                         ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                             >  1:870741/1‑150 (MQ=255)
                                                                                                                         ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                             >  1:871393/1‑150 (MQ=255)
                                                                                                                         ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                             <  2:3640908/150‑1 (MQ=255)
                                                                                                                          tgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtct                            <  2:893792/150‑1 (MQ=255)
                                                                                                                          tgcgttccTGCCCTTGGCTTGCACGCCCGCTAGGGTTGGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtct                            <  1:2186680/150‑1 (MQ=255)
                                                                                                                           gcgttgcTGCCCTTCGCGTGCAAGCCGGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                           >  2:715858/1‑150 (MQ=255)
                                                                                                                           gcgttgcAGCGCTTCGCTTGCACGCCCGGTATGGGTTGAACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                           <  2:1858768/150‑1 (MQ=255)
                                                                                                                           gcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                           >  2:2719147/1‑150 (MQ=255)
                                                                                                                            ggttgcCGCCCTTCGCGTGCACGCCCGCTCCGGGTCCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                          <  1:3595353/149‑1 (MQ=255)
                                                                                                                            cgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                          <  2:1847497/150‑1 (MQ=255)
                                                                                                                             gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                    <  2:2414207/107‑1 (MQ=255)
                                                                                                                             gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                    >  1:2414207/1‑107 (MQ=255)
                                                                                                                             gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                         >  2:108775/1‑150 (MQ=255)
                                                                                                                              gtgcCGCCCTTCGCGGGCACGCCCGCGACGGGTTGCACCGCAACCCCAAAACACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                        <  2:2111105/149‑1 (MQ=255)
                                                                                                                                 caacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                  <  2:560837/134‑1 (MQ=38)
                                                                                                                                 caacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                  >  1:560837/4‑137 (MQ=38)
                                                                                                                                 cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc                     >  1:70054/1‑150 (MQ=255)
                                                                                                                                  aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  2:221151/3‑150 (MQ=38)
                                                                                                                                  aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  1:5445229/1‑150 (MQ=255)
                                                                                                                                  aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  1:2623044/1‑150 (MQ=255)
                                                                                                                                   acccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                   >  2:1590953/2‑150 (MQ=255)
                                                                                                                                     ccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                 >  2:4867387/1‑150 (MQ=39)
                                                                                                                                      cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  1:44375/1‑150 (MQ=255)
                                                                                                                                      cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCGGTACGCTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  1:4498479/1‑150 (MQ=37)
                                                                                                                                      cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  1:4519765/1‑150 (MQ=39)
                                                                                                                                      cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  1:5086871/1‑150 (MQ=39)
                                                                                                                                        tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc              >  2:429711/1‑150 (MQ=255)
                                                                                                                                        tCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc              >  2:3969098/1‑150 (MQ=255)
                                                                                                                                        tCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc              >  2:3968891/1‑150 (MQ=255)
                                                                                                                                          gcgTGCAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg            <  2:353458/150‑1 (MQ=255)
                                                                                                                                               gCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt        <  2:1743684/150‑1 (MQ=255)
                                                                                                                                                 aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt      >  1:1927219/1‑150 (MQ=255)
                                                                                                                                                  aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc     >  1:4266147/1‑150 (MQ=255)
                                                                                                                                                  aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCACCCAGGCCGGCAGTGTc     >  1:2785305/1‑150 (MQ=255)
                                                                                                                                                      ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  >  2:5229887/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGT  >  NC_002947/5745331‑5745628

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