Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,742,452 C→A 27.2% 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,742,4520CA27.2% 195.6 / 89.7 114intergenic (+39/+32)hutF/PP_5037formiminoglutamate deiminase/lipocalin family lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (25/58);  new base A (18/13);  total (43/71)
Fisher's exact test for biased strand distribution p-value = 8.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

Only 100 of 122 total aligned reads displayed.

TGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGT  >  NC_002947/5742304‑5742601
                                                                                                                                                      |                                                                                                                                                     
ggcgcGATGTGAGGGTGAATGGGCAAGGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGccc                                                                                                                                                      <  1:1417528/149‑1 (MQ=255)
  cgcgATGTGATGGTGAAGGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCcgc                                                                                                                                                    <  1:3290729/150‑1 (MQ=255)
     gATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGTCATGCCCAGGAAGATGTAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGaa                                                                                                                                                 <  2:677233/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                                <  1:629242/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                                <  2:1542162/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                                <  1:1420728/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                                <  2:3941962/150‑1 (MQ=255)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAg                                                                                                                                                <  1:1595105/150‑1 (MQ=39)
      aTGTGATGGTGAATGGGCAATGGGTGGTGCACCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGATGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAg                                                                                                                                                <  2:76344/150‑1 (MQ=255)
          gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                            <  2:1832234/150‑3 (MQ=38)
          gATGGTGAAGGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACTGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAgggct                                                                                                                                            <  1:817395/150‑3 (MQ=38)
           ctGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTGTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                           <  2:4457210/149‑1 (MQ=38)
           aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                           <  1:186363/150‑1 (MQ=255)
           aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg                                                                                                                                           <  1:1529658/150‑1 (MQ=255)
           aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg                                                                                                                                           <  2:2050766/150‑1 (MQ=38)
            tGGTGACTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                                          <  1:2661779/150‑1 (MQ=255)
            tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc                                                                                                                                          <  2:1631963/150‑1 (MQ=255)
             ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgca                                                                                                                                         <  2:2843716/150‑1 (MQ=255)
               tGAATGGGCACTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                       <  2:2824006/150‑1 (MQ=38)
               tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGACAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGAGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaac                                                                                                                                       <  2:4499139/150‑1 (MQ=38)
                           gggTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGa                                                                                                                                                  <  1:713682/127‑1 (MQ=38)
                           gggTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGa                                                                                                                                                  >  2:713682/1‑127 (MQ=38)
                                 tgCGCCTGGGGCGGCTTGCCCTGTATGATGAATGCGGGGGGGGTTTTGCCATGTTGTTGCGCCTTGTTTTGGGCTGTGCGTACGGGGTTTGCGTTGCAGCCCTTGGCGTGCAAGCCCGCGATGGTTTGAAACGCAACCCCAAAAAACTCa                                                                                                                     <  2:2142590/150‑1 (MQ=255)
                                 tgCGCCTGGGGCGGCATGCCCTGGATGAGGACAGGGGGGGGGCGTTTGCCATGGTGTGGCGCCTGGTTTTGGGCTGAGGTTACGGGGTTTGCGTTGCAGCCCTTGGCGTGAAAGCCCGCGATGGTTTGAAACGCAACCCCAAAAAACTCa                                                                                                                     <  2:1976364/150‑1 (MQ=255)
                                    ccccggggCGGCTTCCCCAGGACTATCAACGCGGGGGGGCGTGGGCCAGGGTGTTGCGCCAGCTTTTGGCCTGACGGTCCCGGGTTGGCGTTCCCGCCTTTCGGGGGAAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTg                                                                                                                  <  2:3812218/146‑1 (MQ=255)
                                       tggggCGGCTTGCCCAGTAAGATTAACGGGGGGGGGCGTTTGCAATGTTGTTGCGCCAGCTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCAt                                                                                                               <  2:2270914/149‑1 (MQ=255)
                                                  gCCCAGGAAGTTGAACGGGGGCGGGCGTTTGCCATGGTGTTGCGCCAGCTTTTGGGATGAGCGTCCGGGGGTTGGGTTGCCGCCCTCGGGGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTAttt                                                                                                    <  2:4152406/150‑1 (MQ=255)
                                                   tccAGGAAGATGAACGCGGGCGGGCGTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTg                                                                                                   <  2:1040430/149‑1 (MQ=255)
                                                             tGACAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCTTCTTTGGGGCTTAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGACTGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                         <  2:1756725/150‑1 (MQ=255)
                                                             tGAACGGGGGGGGGGTTTTGCCAGGGTGTTGCGCCAGTTTTTGGGGTGAGCGTACCGGGGTTGGGTTGCAGCCCTTCGCGGGCATGCCCGGGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                         <  2:4236771/150‑1 (MQ=255)
                                                             tGAACGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAc                                                                                         <  2:1756706/150‑1 (MQ=255)
                                                              gaccgcggggggGGGTTGGCCAGGGTGTTGCGCCTGCTTTTGGGCGGCGGGTACGGGGGTGGCGTTGCCGCCCTGCGGGTGCAAGCCGGGGATGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcg                                                                                        <  1:721597/146‑1 (MQ=255)
                                                                 aGCGGGCGGGCTTTACCAAAGTTTTTCCCCAACATTTTGGGCTCAGCGAACCGGGGTTGCGTTCCACCCCTTCCCCACCAACCCCCCCAAGGCTTCAAACGCAACCCCAAAAAACTCAGTGCATTATTACTATTTGCTGATCGGAcgccc                                                                                     <  2:4282711/150‑1 (MQ=255)
                                                                     ggcgggcgTTTGCCACGGTGTTGCGCCTGCTTTTGGGCTGAGCGTACCGGGGTTGCGTGGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa                                                                                 <  2:1401108/150‑1 (MQ=255)
                                                                        gggcgTTTGCCATGGTGTTGCGCCAGCTTTTGGGCTGCGCGTCCGGGGGTTGCGTTGCCGCCCTGCGCGGGCACGCCCGGGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGAt                                                                              <  2:2748652/150‑1 (MQ=255)
                                                                           ggTTTGCCATGTTGTTGCGCCAGCTTTTGGGCTGTGCGTAGCGGGGGTGCGTTGCAGCCTTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa                                                                           <  2:3792712/149‑1 (MQ=255)
                                                                              ttGCCAAGGTGTTGCGCCAGCTTTGGGGCGGCGCTTCCGGGGGTTGGGTTCCAGCCCTTCGCGTGCAAGCCCGCCACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCg                                                                        <  2:3058927/150‑1 (MQ=255)
                                                                                  catggTGTTGCGCCAGCTTTTGGGCTGAGCGTCCGGGGTTTGGGTTGCAGCCCTTCGCGTGAAAGCCCGCTACGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTg                                                                    <  1:3674923/147‑1 (MQ=255)
                                                                                   atggTGTTGCGCCAGCTTTTGGGGTGAGCGTACGGGGGTTGGGTTGCCGCCCTTCGCGTGAAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg                                                                   <  2:2712341/148‑1 (MQ=255)
                                                                                             gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTGCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc                                                         <  2:627125/150‑1 (MQ=255)
                                                                                             gcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTGGAGCCCGTTCCGGGCCAACCCCCCGACGGGTGGAAACCAAACCCCAAAAAATCCCTGCCCATTTACTATTTTGAGATAGGACGCCCGCCAGATAAAGCGGCGGGGGTCaaaacc                                                         >  2:4456625/1‑145 (MQ=255)
                                                                                                aacattttGGGCTGAGCGTACCGGTGTTGCGTTGCAGCCCTTCGCGTGAAAGCCCCCCAAGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCct                                                       <  1:358635/146‑1 (MQ=255)
                                                                                                 cgtttttGGGCTGTGGTTACGGGGGTTGGGTTGCAGCCCTTCGCGGCCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCTGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctg                                                      <  2:2267103/147‑1 (MQ=255)
                                                                                                  gcatttGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCctgc                                                     <  2:3605988/147‑1 (MQ=255)
                                                                                                           cTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                            <  2:3760382/150‑1 (MQ=255)
                                                                                                           cTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGcc                                            <  2:4074722/150‑1 (MQ=255)
                                                                                                             gAGCGTACCGGGGTTGCGTTGCCGCCCTGCGCGTGCACGCCCGCGACGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt                                          <  2:3071495/150‑1 (MQ=255)
                                                                                                              cgCGTCCCGGGGTTGCGTCGCCCCCCTCCGCCTGCCCGCCCGCCAAGGCTCCCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTg                                         <  1:4390576/149‑1 (MQ=255)
                                                                                                                gggacgggggTTGGGTGGCCGCCCTTCGCGTGAAAGCCCGCGATGGGTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                       <  2:1558962/144‑1 (MQ=255)
                                                                                                                cttACGGGGTTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                       <  2:3048340/148‑1 (MQ=255)
                                                                                                                cgtagggggtTTGGGTTGCCGCCCTTCGCGGGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGa                                       <  2:767317/144‑1 (MQ=255)
                                                                                                                  tACGGGGGTTGCTTTGCAGCCGTTCGCGGGCACGCCCGCTATGGTTTGAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAca                                     <  1:1035143/150‑1 (MQ=255)
                                                                                                                    ccGGGGTTGCGTTGCCGCCCTTCGCGTGCAACCCCGCTAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGTACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagc                                   <  1:1066556/150‑1 (MQ=255)
                                                                                                                             cgttgcCGCCGTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                          <  1:2613017/150‑1 (MQ=255)
                                                                                                                             cgttgcAGCCCTTCGCGTGCACGCCCGCGACGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                          <  2:611736/150‑1 (MQ=255)
                                                                                                                             cgttgcAGCCCTTCGCGTGCAAGCCCGCGATGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                          <  1:2612914/150‑1 (MQ=255)
                                                                                                                             cgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa                          >  2:1430418/1‑150 (MQ=255)
                                                                                                                              tttgcACCCGTACGGGTGCTAGCCAGCGATGGGTTTAAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                         <  1:2953506/149‑1 (MQ=255)
                                                                                                                              gttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc                         >  2:112673/1‑150 (MQ=38)
                                                                                                                               ttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                        >  1:3576815/1‑150 (MQ=255)
                                                                                                                               ttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACg                        >  2:4302746/1‑150 (MQ=38)
                                                                                                                                 gcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt                      >  2:2617783/1‑150 (MQ=255)
                                                                                                                                  caacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc                     >  1:2072393/4‑150 (MQ=38)
                                                                                                                                  cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc                     >  2:3306598/1‑150 (MQ=255)
                                                                                                                                   aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTGTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  1:2573345/3‑150 (MQ=39)
                                                                                                                                   aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCgg                                                                                           <  1:772865/77‑1 (MQ=18)
                                                                                                                                   aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCgg                                                                                           >  2:772865/3‑79 (MQ=18)
                                                                                                                                   aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  2:2652738/3‑150 (MQ=38)
                                                                                                                                   aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc                    >  1:1244071/3‑150 (MQ=38)
                                                                                                                                   aacccTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGGGCCTTTGTACTATTTGCTGAGCGGACGCCCGCCAGATCACGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGCCAGCAGGTTCTCACGCTcc                    >  2:3228715/3‑150 (MQ=37)
                                                                                                                                    gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACCGCAGGTTCTCACGCTccc                   >  2:3042291/1‑150 (MQ=255)
                                                                                                                                    gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTccc                   >  2:1818092/1‑150 (MQ=255)
                                                                                                                                    gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGGTCTCACGCTccc                   >  2:3042018/1‑150 (MQ=255)
                                                                                                                                     cccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCa                  >  1:76295/1‑150 (MQ=255)
                                                                                                                                      ccTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                 >  1:3525960/1‑150 (MQ=255)
                                                                                                                                      ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                 >  2:2869594/1‑150 (MQ=255)
                                                                                                                                      ccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg                 >  1:2296377/1‑150 (MQ=255)
                                                                                                                                       cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  2:3416090/1‑150 (MQ=255)
                                                                                                                                       cTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg                >  1:1557101/1‑150 (MQ=255)
                                                                                                                                         ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               >  2:789518/1‑150 (MQ=255)
                                                                                                                                         ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               >  1:1520234/1‑150 (MQ=255)
                                                                                                                                         ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               >  2:789352/1‑150 (MQ=255)
                                                                                                                                         ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCGAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               >  1:3180469/1‑150 (MQ=39)
                                                                                                                                         ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               >  1:3081531/1‑150 (MQ=39)
                                                                                                                                         ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc               >  2:4167726/1‑150 (MQ=39)
                                                                                                                                          tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc              >  1:4185346/1‑150 (MQ=255)
                                                                                                                                          tCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc              >  1:1575384/1‑150 (MQ=255)
                                                                                                                                          tCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGcc              <  1:1733582/150‑1 (MQ=39)
                                                                                                                                              gTGCACGCCGGCGAGGGGTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCa          <  1:3663430/150‑1 (MQ=255)
                                                                                                                                               tGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAg         >  1:264488/1‑150 (MQ=255)
                                                                                                                                                gCAAGCCCGCGAAGGTTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgt        <  2:1045429/150‑1 (MQ=255)
                                                                                                                                                   aGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc     >  2:2388941/1‑150 (MQ=255)
                                                                                                                                                    gCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg    >  2:1456416/1‑150 (MQ=255)
                                                                                                                                                    ccccGCGACGGGGAGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg    <  1:963542/149‑1 (MQ=255)
                                                                                                                                                      ccGCGACGGGTTCCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  <  1:4233872/150‑1 (MQ=255)
                                                                                                                                                      ccGCGAAGGGTTGCACCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  <  2:126863/150‑1 (MQ=255)
                                                                                                                                                      ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  >  1:3851438/1‑150 (MQ=255)
                                                                                                                                                      ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  >  2:156195/1‑150 (MQ=255)
                                                                                                                                                      ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTAATATTTGCTGATCGGAAGCACGCCAGATCAAGCGGCTGTTGTCATAAAAATGCTGCAGAGCAATGTACAGCAGGTTCTCACGCTCCAAGGACGGCAGAGTCGgt  >  2:630893/1‑150 (MQ=255)
                                                                                                                                                      ccGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTATACTATTTGCTGATCTGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt  >  2:3314519/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGT  >  NC_002947/5742304‑5742601

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: