Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,745,489 T→C 9.4% intergenic (+49/+22) 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,4890TC 9.4% 360.6 / 21.4 128intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (44/71);  new base C (10/2);  total (55/73)
Fisher's exact test for biased strand distribution p-value = 4.12e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

Only 100 of 136 total aligned reads displayed.

TGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACA  >  NC_002947/5745340‑5745638
                                                                                                                                                        |                                                                                                                                                      
tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                                        <  1:829443/150‑1 (MQ=255)
tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                                        <  1:1929516/150‑1 (MQ=255)
 gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                       <  2:4149516/150‑1 (MQ=255)
 gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                       <  1:4076772/150‑1 (MQ=255)
 gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                       <  1:3234471/150‑1 (MQ=255)
  aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                      <  1:2016480/150‑1 (MQ=255)
    ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGGCCGTCGCGGGCCAGCCCGCGAAGGGGTgca                                                                                                                                                    >  2:6184048/1‑150 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  2:3522027/150‑1 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  1:4610857/150‑1 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                                 <  1:164880/150‑1 (MQ=38)
                       gggcgcCTGGGGGGTCAGGCCCAGTACTATGAACGCGGGCGGGCTTTTGCCACGGTGTTGCGCATGCTTTGGGGCTTCGCGTACGGGGGTTGCGTCGACGCCGTTCGCGTGCACGCCCGCGTCGGGTTGAAACGCAACCCCAAAAAACTc                                                                                                                                <  2:1921778/148‑1 (MQ=255)
                          cgcCTGGGGCGGATTGCCCTGGAAGATGAACGCGGGCGGGGGTTGGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGGGAAGGTTTGAAACGCAACCCCAAAAAACTCAGt                                                                                                                             <  2:456389/150‑1 (MQ=255)
                           gcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGAGCCCCAAGCCGGGCAACCCCCCCAAGGGTGGCAACCAAACCCCaaaaaaaaaagg                                                                                                                             >  2:6031636/1‑144 (MQ=255)
                                           ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTGGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                                            <  2:5119977/150‑1 (MQ=255)
                                                aaGATGAACGCGGGGGGGCGTTTGCCACGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACCGGGGTGGCGTGGCCGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATc                                                                                                       <  2:4737678/150‑1 (MQ=255)
                                                  gATGAACGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCgg                                                                                                     <  2:5494885/150‑1 (MQ=255)
                                                  gATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                                 <  2:4905750/106‑1 (MQ=25)
                                                            ggggggcgTTTGCCACGGTGTTGCGCCAGCTTTGGGGCTGAGCGTACCGGGGTGGCGTTGCAGCCCTCCGCGTGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                           <  2:4875822/147‑1 (MQ=255)
                                                               gggcgCTCGCCACGGTGGTCCGCCAGCTTTTGGGCGGTGCGTCCCGGGGTGGCGTGCCAGCCCTCCGCGTGCACGCCCGCGAAGGTTCCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGAt                                                                                        <  1:6705665/150‑1 (MQ=255)
                                                                    ttggCCAAGTTGTTTCGCCCGTTTTTGGGCTTAGCGTACCGGTGTTGCTTTGCAGCCATGCGCGTGAAACCCCGCGATGTTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc                                                                                   <  2:4041669/147‑1 (MQ=255)
                                                                        ccACGGTGTTGCGCCCGCTTTGGGGCTGCGCGCCCGGGGGTTGCGTTGCAGCCCTGCGCGGGCACGCCCGCGACGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCt                                                                               <  2:3067956/150‑1 (MQ=255)
                                                                         catggTGTGGCGCCAGCTTTGGGGCTGAGCGTACGGGGGTTGCGTGGCAGCCCTTCGGTTGCAAGCCCGCGACGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                              <  2:5813369/147‑1 (MQ=255)
                                                                           tggTGTGGCGCCAGCTTTTGGGCTGCGCGTACGGGGGTTGCTTGGCAGCCCTTCGCTTGCACGCCCGCTAGGGTTCGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  1:518598/149‑1 (MQ=255)
                                                                           aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  2:5062013/150‑1 (MQ=255)
                                                                           aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  2:6644608/150‑1 (MQ=255)
                                                                           aGGTGTTGCGCCAGCTTTTGGGCTGAGCGGACCGGGGTTGCGTTGCAGCCGTTCGCGTGCAAGCCCGCGAAGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt                                                                            <  2:3879032/150‑1 (MQ=255)
                                                                                   cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAcc                                                                    <  2:2187860/150‑1 (MQ=255)
                                                                                     cctgCTTTTGGGGTGAGCGTACGGGGGTTGCGTTGCAGCCCTGCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAcccc                                                                  <  2:3055532/147‑1 (MQ=255)
                                                                                      cAGCTTTTGGGCTGAGCGTCCCGGGGTTGCGTTGCAGCCCTTCGCGTGCACCCCCGCGACGGGTTGAACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                                 <  2:3142643/150‑1 (MQ=255)
                                                                                      cAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                                 <  2:3364619/150‑1 (MQ=255)
                                                                                              gggCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCcgc                                                         <  2:4312844/150‑1 (MQ=255)
                                                                                                 cTGAGCTTACGGGGTTTGCGTGCCCCCCCTCCGCGTCCACCCCCGCAACGGGTACAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                                      <  1:3702124/150‑1 (MQ=255)
                                                                                                 cTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCATCCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                                      <  2:1220782/150‑1 (MQ=255)
                                                                                                   gAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                                    <  1:512179/150‑1 (MQ=255)
                                                                                                    aGCGTACCGGGGTTGCGTTGCAGCCCTTGGCGTGCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                                   <  2:3088112/150‑1 (MQ=255)
                                                                                                     tggttccGGGGTTGGGTTTGCAGCCTTTCGGTTCAACGCCCGGTTTGGTTTCAACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                                   <  1:1877688/145‑1 (MQ=255)
                                                                                                      ggTACGGGGTTGGGGTTGCAGCCCTTGGGTTGAAAGCCCGGTATGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                                 <  1:5462751/149‑1 (MQ=255)
                                                                                                      cGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGCTCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                                 <  2:3304167/150‑1 (MQ=255)
                                                                                                        gacgggggTTGGGTTGCAGCCCTTCGGGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                                               <  1:5663726/146‑1 (MQ=255)
                                                                                                          ccGGGGTTGCGTTGCAGCCCTTCGCTTGCAAGCCCGCTACGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                                             <  2:1941188/150‑1 (MQ=255)
                                                                                                          ccGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                                             <  2:5338887/150‑1 (MQ=255)
                                                                                                           gggggTTGGGTTGCCCCCGTTGGGGTCCACCCCCGCGATGGTTTCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                            <  1:6514349/149‑1 (MQ=255)
                                                                                                           cGGGGTTGCGTTGCAGCCCTTCGCGTCCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                            <  2:3146116/150‑1 (MQ=255)
                                                                                                            ggggTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag                                           <  1:3241440/150‑1 (MQ=255)
                                                                                                                ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                                       >  2:560303/1‑150 (MQ=255)
                                                                                                                ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                                       >  1:3568982/1‑150 (MQ=255)
                                                                                                                   cgttgcCGCCCTTCGCGCGCACCCCCGCGACGGGTGGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                                    <  2:295670/150‑1 (MQ=255)
                                                                                                                    gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCGGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                                   >  2:1422568/1‑150 (MQ=255)
                                                                                                                     gttcccgCCCTCGGCGTGCACGCCCGCTTTGGGTCCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                                   <  1:6016330/144‑1 (MQ=255)
                                                                                                                          cgCCCTGCGGGTGCACCCCGGCGCCGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              <  1:4094376/149‑1 (MQ=255)
                                                                                                                          aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              >  1:812109/3‑150 (MQ=38)
                                                                                                                          aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                              >  1:566494/3‑150 (MQ=38)
                                                                                                                           gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                             >  2:6364199/1‑150 (MQ=255)
                                                                                                                           gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGAGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                             >  2:2122964/1‑150 (MQ=255)
                                                                                                                             ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAATACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                           >  2:2654366/1‑150 (MQ=39)
                                                                                                                             ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                           >  1:4315040/1‑150 (MQ=255)
                                                                                                                             ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCAGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                           >  1:4315068/1‑150 (MQ=39)
                                                                                                                             ccTTCGCGGGCTTGCACGCGAAGGGCTGCAAAGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                           >  1:3447869/1‑150 (MQ=39)
                                                                                                                              cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  2:2119840/1‑150 (MQ=255)
                                                                                                                              cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:6777906/1‑150 (MQ=255)
                                                                                                                              cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:4544556/1‑150 (MQ=255)
                                                                                                                              cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  2:1308051/1‑150 (MQ=255)
                                                                                                                              cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:5815131/1‑150 (MQ=255)
                                                                                                                               ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:3635594/1‑150 (MQ=255)
                                                                                                                               ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAATATCTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:6680997/1‑150 (MQ=255)
                                                                                                                               ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:1429714/1‑150 (MQ=255)
                                                                                                                               ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:5122465/1‑150 (MQ=255)
                                                                                                                               ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:5729516/1‑150 (MQ=255)
                                                                                                                                tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc                        <  2:5541466/150‑1 (MQ=255)
                                                                                                                                 tcgggcgCATGCCCGCGACGGGTTACAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc                        <  2:4593252/144‑1 (MQ=255)
                                                                                                                                 cgcgTGCACGCCCGCGACGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCg                       <  2:1448171/150‑1 (MQ=255)
                                                                                                                                  gcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAATAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg                      >  2:4102243/1‑150 (MQ=255)
                                                                                                                                  gcgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCgg                      >  1:2157172/1‑150 (MQ=255)
                                                                                                                                    cgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGc                     >  1:6279211/1‑150 (MQ=255)
                                                                                                                                    cgTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGc                     >  1:2977007/1‑150 (MQ=255)
                                                                                                                                     gTGCAAGCCCGCTAAGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa                    <  2:4697922/150‑1 (MQ=255)
                                                                                                                                       gCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt                  <  2:5772046/150‑1 (MQ=255)
                                                                                                                                        gcttgcacgcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt                  >  2:5671740/8‑150 (MQ=255)
                                                                                                                                        gcttgcacgcgAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt                  >  1:3082641/8‑150 (MQ=255)
                                                                                                                                        cAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg                 <  2:3899636/150‑1 (MQ=255)
                                                                                                                                         cacgCCCGCGACGGGTTGCAACGCAACCCCAAACAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg                 <  2:6429047/147‑1 (MQ=255)
                                                                                                                                         aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt                >  1:1692757/1‑150 (MQ=255)
                                                                                                                                         aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt                >  2:3504811/1‑150 (MQ=255)
                                                                                                                                          aGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc               <  1:6729223/150‑1 (MQ=255)
                                                                                                                                          aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc               <  1:2781524/150‑1 (MQ=255)
                                                                                                                                          aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc               <  2:5542944/150‑1 (MQ=255)
                                                                                                                                             cccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCgg             >  1:3341911/1‑150 (MQ=255)
                                                                                                                                              cggctatgggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            <  1:2360867/143‑1 (MQ=255)
                                                                                                                                              ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            <  1:5862279/150‑1 (MQ=255)
                                                                                                                                              ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            >  1:5336472/1‑150 (MQ=255)
                                                                                                                                                gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt          >  2:5006309/1‑150 (MQ=255)
                                                                                                                                                gcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt          >  1:3936666/1‑150 (MQ=255)
                                                                                                                                                gcgAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt          >  1:1727804/1‑150 (MQ=255)
                                                                                                                                                  gAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgc        >  2:5726385/1‑150 (MQ=255)
                                                                                                                                                  gAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgc        <  2:3113636/150‑1 (MQ=255)
                                                                                                                                                      ggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa    >  2:6736182/1‑150 (MQ=255)
                                                                                                                                                       gTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGAc   >  1:6591822/1‑150 (MQ=255)
                                                                                                                                                       gTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGAc   >  2:1090732/1‑150 (MQ=255)
                                                                                                                                                        ttGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACa  >  1:4494214/1‑150 (MQ=255)
                                                                                                                                                        ttGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACa  <  2:3016333/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
TGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACA  >  NC_002947/5745340‑5745638

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: