Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,745,489 T→C 8.2% 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 8.2% 389.5 / 17.9 135intergenic (+49/+22)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (43/78);  new base C (9/2);  total (55/80)
Fisher's exact test for biased strand distribution p-value = 6.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 145 total aligned reads displayed.

TGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACA  >  NC_002947/5745340‑5745638
                                                                                                                                                        |                                                                                                                                                      
tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGTGTGCAAGCCCGCGAAGGGt                                                                                                                                                        <  1:342631/150‑1 (MQ=255)
tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGt                                                                                                                                                        <  1:6333057/150‑1 (MQ=255)
 gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                       <  1:1092254/150‑1 (MQ=255)
 gATGGTGAAGGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAGGCCCGCGAAGGGtt                                                                                                                                                       <  1:502214/150‑1 (MQ=255)
  aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                      <  1:1887122/150‑1 (MQ=255)
  aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                                      <  1:4399839/150‑1 (MQ=255)
     gTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAGGGCTgcaa                                                                                                                                                   <  2:2324198/150‑1 (MQ=38)
      tGACTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  2:187347/150‑1 (MQ=255)
      tGAATGGGCACTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                                 <  1:2469103/150‑1 (MQ=38)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGAGCCCCTAGCGGGCAAGCCCGCCCACGGGTgcaac                                                                                                                                                 >  2:5612096/1‑150 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  2:187521/150‑1 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  1:311425/150‑1 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTAGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  1:311155/150‑1 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACAGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  2:579705/150‑1 (MQ=255)
      tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACAGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaac                                                                                                                                                 <  2:4927344/150‑1 (MQ=255)
           gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                       >  1:4369784/1‑140 (MQ=255)
           gggCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt                                                                                                                                                       <  2:4369784/140‑1 (MQ=255)
             gCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTTGGATTGAGCCCCTAGCCAGCAAGCCCCCGCAGGGGTGGAACGCAAccc                                                                                                                                          >  2:560412/1‑150 (MQ=255)
              cAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTTGTGCCCCACCGCGGGCAGCCCCCCAAAGGGTTGAAACCAAAccc                                                                                                                                          >  2:3079714/1‑150 (MQ=255)
                     tggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGGGGTGGAGCCCATCCCGGGGCAACCCCCCCAGGGGTGGAAACCCCCCCCCaaaaaa                                                                                                                                   >  2:5974901/1‑150 (MQ=255)
                       gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTTGAGCCCCTCGCCGGGAACCCCCCCAAGCGGGTCAACCCCACCCCACAAAAcaa                                                                                                                                >  2:7076473/1‑148 (MQ=255)
                        ggcgcCAGGGGCGGCATGCCCAGGAAGATGAACGCGGGCGGGCTTTGGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCa                                                                                                                               <  2:5071634/149‑1 (MQ=255)
                             ctggggCGGATGGCCCAGGACTATGAATGCGGGCGGGCTTTTGCAATGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACCGGGTTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCa                                                                                                                          <  2:2707035/148‑1 (MQ=255)
                              cggggCGGCATGCCCAGGACGATGAAAGCGGGCGGGCGTTGGCCTTGGTGTTGCGCCAGTTTTTGGGGTGTGGGTACGGGGGTTGCTTGGCCGCCGTTGGGGTGAAAGCCCGCGAGGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAt                                                                                                                         <  2:1234677/149‑1 (MQ=255)
                              aGGGGCGGCATGCCCAGGAAGTTGAACGCGGGCGGGCTTTTGCCAAGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAt                                                                                                                         <  2:4229533/150‑1 (MQ=255)
                                 ggcggcATGCCCAGGACGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAtttt                                                                                                                      <  2:5644151/150‑1 (MQ=255)
                                  gcggaATGCCCAGTAATATGAAAGCGGGCGGGCTTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTCAGCGTACGGGGGTTGCTTTGCAGCCCTTCGCTTCAAAGCCCGCTATGGTTTGAAACGCAACCCCAAAAAACTCAGTGCAttttt                                                                                                                     <  2:7652460/150‑1 (MQ=255)
                                    ggcATGCCCAGTAAGATGAACGGGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTAc                                                                                                                   <  2:5641084/150‑1 (MQ=255)
                                           ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGCGCGTACCGGGGTGGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc                                                                                                            <  2:4363685/150‑1 (MQ=255)
                                               gatgATGAAAGCGGGCGGGCGTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGAGCTTACCGGGGTTGCGTTGCAGCCGTTCGCGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGAt                                                                                                        <  2:5022871/147‑1 (MQ=255)
                                                    tGAACGGGGGCGGGTTTTTTCCATGTTTTTGCGCCTGTTTTTGGGCTGTGCTTACGGGGTTTGGTTTGCAGCCGTTCGCGTGAAAGCCGGGGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                                   <  2:1017383/150‑1 (MQ=255)
                                                          ggggcggggGTTTGCCATGGTGTGGCGCCAGCTTTGGGGCTGCGCGTCCGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgc                                                                                             <  1:4357927/149‑1 (MQ=255)
                                                           ggggggggGTTTGCCACGGTGTGGCGCCTGATTTTGGGCTTCGCGTCCGGGGTTTGGGTTGCCGCCCTGCGCGTGAACGCCCGCGATGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc                                                                                            <  1:4231650/150‑1 (MQ=255)
                                                            ggggggcgTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCATGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                           <  2:8380317/147‑1 (MQ=255)
                                                            ggcgggcgTTTGCCACGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                           <  2:8496426/150‑1 (MQ=255)
                                                            ggcgggcgTTTGCAATGGTGTTGCCCCAGATTTTGGGCTGCGCGTACGGGGGTTGCGTCGCAGCCCTTCGCGTCAAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                           <  2:1842150/150‑1 (MQ=255)
                                                                  cgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                                     <  2:6504539/150‑1 (MQ=255)
                                                                  cgTTGGCCCCGGTGTGGCGCCAGCTTTGGGGCTGCGCGTACGGGGGTGGCGTTGCCGCCCTGCGCGTGCAAGCCGGCGACGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                                     <  2:309226/150‑1 (MQ=255)
                                                                       gccttggTTTTGCGCCTGTTTTTGGGGTGTGGTTCGGGGGTTGGGTTTGCCCCCTTGCGGGTGCACGCCCGGTTTGGTTTCAACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGc                                                                                <  1:6291716/145‑1 (MQ=255)
                                                                         cacggTGTTGCGCCAGCTTTTGGGCTGCGCGTACCGGGGTTGCGTGGCCGCCCTTCGCGGGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAACACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                              <  2:7690590/147‑1 (MQ=255)
                                                                         cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                              <  2:7190033/150‑1 (MQ=255)
                                                                          atggTGTTGCGCCAGCTTTTGGGGTGCGCGTCCGGGGGTTGCGTTGCAGCCCTTCGCTGGCACGCCGGCTACGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                                                             <  1:1436045/148‑1 (MQ=255)
                                                                             gtgtTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtgt                                                                          <  2:3494362/150‑1 (MQ=255)
                                                                              ttttGCCCCAGATTTGGGGCTGTGCGTAGGGGGGTTGCGTTGCCGCCTTGGGGGTGCACGCCCGCTTTGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc                                                                         <  1:6977240/148‑1 (MQ=255)
                                                                               gtTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGTTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCa                                                                        <  2:3460519/150‑1 (MQ=255)
                                                                                tggcgcCTGCTTTTGGGGTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGAAAGCCGGCTATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCAt                                                                       <  2:6686916/148‑1 (MQ=255)
                                                                                     ccAGCTTTTGGGCTGTGCGTACCGGGTTTGCGTTGCAGCCCTTCGCTTGCACGCCCGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAcccc                                                                  <  2:4669697/150‑1 (MQ=255)
                                                                                       aGCTTTTGGGCTGAGCGTACGGGGTTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctg                                                                <  2:6907686/150‑1 (MQ=255)
                                                                                         cTTTGGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCACGCCCGCTACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgct                                                              <  1:5858898/150‑1 (MQ=255)
                                                                                            ttGGGCGGAGGGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTGGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc                                                           <  2:6205128/150‑1 (MQ=255)
                                                                                                gCTGAGCGTACCGGGGTTGCGTTTCAGCCCGTAGCGTGCAAGCCCGCGTAGGGTTGCAAAGCAACCCCAAAAAAATAAGTGCAATTTTTATATTTGGGGTTAGGACCCCCGCCACAACAAGCCGCTGGGGTTATAACCCCGCTTCcgccc                                                       >  2:8042721/1‑148 (MQ=255)
                                                                                                 cTGCGCGTACGGGGGTTGCGTTGCACCCCTTCGCGGGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                                      <  2:8105173/150‑1 (MQ=255)
                                                                                                         aCCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcag                                              <  2:8010340/150‑1 (MQ=255)
                                                                                                           cGGGTTTGCTTTGCACCCCTTCGCGTGAAAGCCCGCTACGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca                                            <  1:2813176/150‑1 (MQ=255)
                                                                                                           cGGGGTTGGCTTTGCAGCCCTTCGCGGGAACGCCCGCTATGGTTTCAACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                                             <  1:4849482/150‑1 (MQ=255)
                                                                                                              ggTTGGGTGGCCGCCCTGCGCGTCCATGCCGGCTATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGt                                         <  1:4884859/150‑1 (MQ=255)
                                                                                                                ttgcgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                                       >  2:1281206/1‑150 (MQ=255)
                                                                                                                 tggcgttgcAGCCCTTCGCGTGCAAGCCCGCGACGGTTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc                                       <  2:164127/148‑1 (MQ=255)
                                                                                                                   gcgtttcAGCCCTTGGGGTGAACGCCGGCTTGGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                                     <  1:5664107/150‑1 (MQ=255)
                                                                                                                   gcgttgcCGCCCTTCGCGTGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                                     <  2:3047355/150‑1 (MQ=255)
                                                                                                                    gctttgcCGCCCTTCGCGTGCACGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtctc                                     <  2:3047338/147‑1 (MQ=255)
                                                                                                                       ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAATCTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                                  >  1:2930515/1‑150 (MQ=255)
                                                                                                                       ttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                                  >  1:6681890/1‑150 (MQ=38)
                                                                                                                         gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAGCTCTGTGCATGTTTACTATTTGCTGATCGGGCGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGCCCGCAGGTTCTCACGCt                                >  2:5727212/1‑150 (MQ=255)
                                                                                                                            gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                             >  1:6368391/1‑150 (MQ=255)
                                                                                                                             cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa                            >  2:3914425/1‑150 (MQ=255)
                                                                                                                             cccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                                               <  2:6668/131‑1 (MQ=255)
                                                                                                                              ccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                           >  2:2787275/1‑150 (MQ=255)
                                                                                                                              ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                           >  1:7940559/1‑150 (MQ=255)
                                                                                                                               cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:6314887/1‑150 (MQ=255)
                                                                                                                               cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  2:2616860/1‑150 (MQ=255)
                                                                                                                               cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:2066046/1‑150 (MQ=255)
                                                                                                                               cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:1949778/1‑150 (MQ=255)
                                                                                                                               cTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  1:759121/1‑150 (MQ=255)
                                                                                                                               cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                          >  2:1348289/1‑150 (MQ=39)
                                                                                                                                ttCGGGTGCACGCCCGCTATGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         <  2:3147068/150‑1 (MQ=255)
                                                                                                                                ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCTAGGc                         >  1:8135948/1‑150 (MQ=255)
                                                                                                                                ttCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:7317262/1‑150 (MQ=255)
                                                                                                                                ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         >  2:4667353/1‑150 (MQ=39)
                                                                                                                                 tcggggTGCACCCCCGCTATGGGTTCAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc                         <  1:7693780/145‑1 (MQ=255)
                                                                                                                                 tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc                        >  1:1021544/1‑150 (MQ=255)
                                                                                                                                 tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc                        >  2:6152905/1‑150 (MQ=255)
                                                                                                                                 tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc                        >  2:3564274/1‑150 (MQ=255)
                                                                                                                                      gTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa                    >  1:5603313/1‑150 (MQ=255)
                                                                                                                                      gTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa                    <  2:8580846/150‑1 (MQ=255)
                                                                                                                                          aaGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtgt                <  1:7967426/150‑1 (MQ=255)
                                                                                                                                            gCCCGCGTTGGGTCCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg              <  1:5496773/150‑1 (MQ=255)
                                                                                                                                            gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg              <  2:6317747/150‑1 (MQ=255)
                                                                                                                                            gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg              >  1:2921381/1‑150 (MQ=255)
                                                                                                                                              ccGCGATGGTTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            <  1:5612096/150‑1 (MQ=255)
                                                                                                                                              ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            >  1:1974385/1‑150 (MQ=255)
                                                                                                                                              ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt            >  1:6189644/1‑150 (MQ=255)
                                                                                                                                               cgcgAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtg           >  1:6721420/1‑150 (MQ=255)
                                                                                                                                                  gAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgc        >  1:6901764/1‑150 (MQ=255)
                                                                                                                                                    ggggTTGCGTTGCAACCCCAAAAAACTCAGGGCATTTTTACTATTTGCTGATCGGACCCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgcgc      >  2:6761121/2‑150 (MQ=255)
                                                                                                                                                      ggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                                    <  2:3359128/118‑1 (MQ=255)
                                                                                                                                                      ggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                                    >  1:3359128/1‑118 (MQ=255)
                                                                                                                                                      ggTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa    >  1:4111350/1‑150 (MQ=255)
                                                                                                                                                       gTTCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGAc   <  1:2423721/150‑1 (MQ=255)
                                                                                                                                                        ttGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGACa  <  2:1347035/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 < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: