Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,750,086 A→C 8.9% intergenic (+26/+29) yfjD → / ← grpE UPF0053 family inner membrane protein/heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,750,0860AC 8.9% 255.8 / 17.3 79intergenic (+26/+29)yfjD/grpEUPF0053 family inner membrane protein/heat shock protein
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base C (4/3);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

GTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCAT  >  NC_000913/2749938‑2750234
                                                                                                                                                    |                                                                                                                                                    
gTCGCAGGCACCCGCGTGCTTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCa                                                                                                                                                      <  4:616326/149‑1 (MQ=255)
gTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCa                                                                                                                                                      <  4:419799/149‑1 (MQ=255)
 tCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCa                                                                                                                                                      <  2:758397/148‑1 (MQ=255)
  cGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAg                                                                                                                                                     <  1:294786/148‑1 (MQ=255)
    cAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAggg                                                                                                                                                   <  1:722785/148‑1 (MQ=255)
     aGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGcc                                                                                                                                                 >  3:92171/1‑149 (MQ=255)
            cgcgTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatt                                                                                                                                          >  3:808166/1‑149 (MQ=255)
             gcgtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatta                                                                                                                                         >  3:411036/1‑149 (MQ=255)
             gcgtgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAatta                                                                                                                                         >  3:602842/1‑149 (MQ=255)
                tgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCgg                                                                                                                                                    <  4:263405/135‑1 (MQ=255)
                tgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCgg                                                                                                                                                    >  3:263405/1‑135 (MQ=255)
                tgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                                      <  4:602842/149‑1 (MQ=255)
                tgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                                      <  4:275434/149‑1 (MQ=255)
                 gcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                                      <  1:8661/148‑1 (MQ=255)
                     aTTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaa                                                                                                                                  >  2:196760/1‑148 (MQ=255)
                       tGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAg                                                                                                                                >  2:424104/1‑148 (MQ=255)
                        ggCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCa                                                                                                                              <  4:486495/149‑1 (MQ=255)
                              tACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAAtt                                                                                                                        >  4:717334/1‑149 (MQ=255)
                                cgatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                                      <  3:683525/149‑1 (MQ=255)
                                cgatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                                      <  4:320166/149‑1 (MQ=255)
                                cgatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                                      <  4:411036/149‑1 (MQ=255)
                                 gatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAg                                                                                                                     <  3:227234/149‑1 (MQ=255)
                                  atatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAg                                                                                                                     >  2:709222/1‑148 (MQ=255)
                                       gatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                                      <  3:176157/126‑1 (MQ=255)
                                       gatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                                      >  4:176157/1‑126 (MQ=255)
                                        atatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCttt                                                                                                                 <  2:98582/146‑1 (MQ=255)
                                        atatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCttt                                                                                                                 >  1:98582/1‑146 (MQ=255)
                                        atatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGc                                                                                                              <  3:717334/149‑1 (MQ=255)
                                          atTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCt                                                                                                             <  1:415035/148‑1 (MQ=255)
                                          atTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCt                                                                                                             >  1:434199/1‑148 (MQ=255)
                                            tctcGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCttt                                                                                                           <  1:424104/148‑1 (MQ=255)
                                                gACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCt                                                                                                       >  1:56237/1‑148 (MQ=255)
                                                   gTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCGGGcc                                                                                                                                                 <  3:181478/103‑1 (MQ=39)
                                                   gTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCGGGcc                                                                                                                                                 >  4:181478/1‑103 (MQ=38)
                                                      cAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTa                                                                                                 <  1:709222/148‑1 (MQ=255)
                                                       aGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAAc                                                                                               >  4:574611/1‑149 (MQ=255)
                                                        ggACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAAt                                                                                                                         >  1:13508/1‑122 (MQ=255)
                                                        ggACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAAt                                                                                                                         <  2:13508/122‑1 (MQ=255)
                                                        ggACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAAcc                                                                                              >  3:315349/1‑149 (MQ=255)
                                                        ggACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAAcc                                                                                              <  4:693343/149‑1 (MQ=255)
                                                         gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                                      <  3:699834/124‑1 (MQ=255)
                                                         gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                                      >  4:699834/1‑124 (MQ=255)
                                                                 gATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATcgccac                                                                                      <  2:56237/148‑3 (MQ=255)
                                                                 gATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCa                                                                                     <  4:385939/149‑1 (MQ=255)
                                                                  aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCAc                                                                                    >  4:598716/1‑149 (MQ=255)
                                                                  aTTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCAc                                                                                    >  4:710198/1‑149 (MQ=255)
                                                                   ttAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCAc                                                                                    >  2:380883/1‑148 (MQ=255)
                                                                       gCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGc                                                                                                        >  3:303239/1‑124 (MQ=255)
                                                                       gCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGc                                                                                                        <  4:303239/124‑1 (MQ=255)
                                                                       gCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATc                                                                               >  3:473305/1‑149 (MQ=255)
                                                                           gTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTAc                                                                           >  4:742463/1‑149 (MQ=255)
                                                                              aaaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGt                                                                                                  >  3:312226/1‑123 (MQ=255)
                                                                              aaaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGt                                                                                                  <  4:312226/123‑1 (MQ=255)
                                                                               aaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTa                                                                                                                       <  1:409434/101‑1 (MQ=255)
                                                                               aaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTa                                                                                                                       >  2:409434/1‑101 (MQ=255)
                                                                                   ttttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCAtt                                                                    >  1:729357/1‑148 (MQ=255)
                                                                                     ttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCACGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAg                                                                 <  4:255966/149‑1 (MQ=255)
                                                                                      tCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAg                                                                 <  1:196760/148‑1 (MQ=255)
                                                                                         tgtgAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTTAGCGt                                                              >  2:566139/1‑148 (MQ=255)
                                                                                         tgtgAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGt                                                              >  2:261485/1‑148 (MQ=255)
                                                                                            gAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGtata                                                           >  1:367441/1‑148 (MQ=255)
                                                                                               aCCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAAccc                                                       >  4:115898/1‑149 (MQ=255)
                                                                                                 cGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCtt                                                     <  4:315349/149‑1 (MQ=255)
                                                                                                      cgcgAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCa                                                <  4:92171/149‑1 (MQ=255)
                                                                                                         gagaGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATc                                                                               <  3:637583/115‑1 (MQ=255)
                                                                                                         gagaGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATc                                                                               >  4:637583/1‑115 (MQ=255)
                                                                                                             gtgtGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGcc                                         <  3:691758/149‑1 (MQ=255)
                                                                                                                tgGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAg                                      <  4:756761/149‑1 (MQ=255)
                                                                                                                   cGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTa                                    >  1:503780/1‑148 (MQ=255)
                                                                                                                           cGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCt                            >  1:363446/1‑148 (MQ=255)
                                                                                                                            gAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTgg                          >  3:691366/1‑149 (MQ=255)
                                                                                                                             ttcaCGGCCCTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgc                         >  4:343114/4‑149 (MQ=255)
                                                                                                                             ttcaCGGCCCTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgc                         <  3:343114/146‑1 (MQ=255)
                                                                                                                                   gCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAAc                    >  1:288567/1‑148 (MQ=255)
                                                                                                                                    cccTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACg                   >  2:848322/1‑148 (MQ=255)
                                                                                                                                    cccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGt                  >  3:426481/1‑149 (MQ=255)
                                                                                                                                     ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtc                 >  3:508957/1‑149 (MQ=255)
                                                                                                                                          aTTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAg            >  4:186488/1‑149 (MQ=255)
                                                                                                                                              gCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTc        <  4:426481/149‑1 (MQ=255)
                                                                                                                                               cATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTc        <  1:380883/148‑1 (MQ=255)
                                                                                                                                                 tGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTccac     <  4:473305/149‑1 (MQ=255)
                                                                                                                                                  gCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTccac     >  2:793896/1‑148 (MQ=255)
                                                                                                                                                    aGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCAt  >  3:532817/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTCATCAGATTCCACCAT  >  NC_000913/2749938‑2750234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: