Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,440,960 A→C intergenic (‑176/‑98) pfo ← / → ydbJ pyruvate‑flavodoxin oxidoreductase/DUF333 family putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,440,9600AC98.9% 333.9 / ‑4.1 88intergenic (‑176/‑98)pfo/ydbJpyruvate‑flavodoxin oxidoreductase/DUF333 family putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (42/45);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

GGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCG  >  NC_000913/1440813‑1441105
                                                                                                                                                   |                                                                                                                                                 
ggCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATg                                                                                                                                                <  1:1047459/151‑1 (MQ=255)
   gcACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCgg                                                                                                                                             >  2:44447/1‑151 (MQ=255)
    cACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCggg                                                                                                                                            <  2:799606/151‑1 (MQ=255)
     aCGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCgggg                                                                                                                                           <  2:203933/151‑1 (MQ=255)
      cGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAGTGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGc                                                                                                                                          <  2:1087192/151‑1 (MQ=255)
       gAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTCGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGc                                                                                                                                          >  2:478281/1‑150 (MQ=255)
       gAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCa                                                                                                                                         <  2:1203105/151‑1 (MQ=255)
       gAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCa                                                                                                                                         <  2:664464/151‑1 (MQ=255)
        aaaTTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCaa                                                                                                                                        <  2:726041/151‑1 (MQ=255)
          aTTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAAcg                                                                                                                                      <  2:49205/151‑1 (MQ=255)
            tgctgcGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAAcgc                                                                                                                                     <  2:1230877/150‑1 (MQ=255)
            tgctgcGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAAcgc                                                                                                                                     <  2:1202760/150‑1 (MQ=255)
              ctgcGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACgcgc                                                                                                                                   >  2:220858/1‑150 (MQ=255)
                gcgcgcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCa                                                                                                                                  <  1:192619/149‑1 (MQ=255)
                 cgcgcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAgg                                                                                                                                >  2:1193833/1‑150 (MQ=255)
                    gcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGAc                                                                                                                              >  2:445368/1‑149 (MQ=255)
                    gcCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAg                                                                                                                            <  2:226571/151‑1 (MQ=255)
                      ccAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGct                                                                                                                          >  2:1182728/1‑151 (MQ=255)
                      ccAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGct                                                                                                                          >  2:275394/1‑151 (MQ=255)
                          tagtaATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCtctc                                                                                                                       <  1:918219/150‑1 (MQ=255)
                          tagtaATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCtctc                                                                                                                       <  2:1265207/150‑1 (MQ=255)
                              aaTCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGAcg                                                                                                                   <  1:510858/150‑1 (MQ=255)
                                tCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGAcgccg                                                                                                                >  1:644885/1‑151 (MQ=255)
                                 cTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCCGGACAGCTCTCGAcgccgc                                                                                                               >  2:716552/1‑151 (MQ=255)
                                   ttCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGAcgcc                                                                                                                 <  2:562080/147‑1 (MQ=255)
                                      aaTTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCctct                                                                                                          <  2:405664/151‑1 (MQ=255)
                                      aaTTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCctct                                                                                                          <  1:399228/151‑1 (MQ=255)
                                         tttAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCtgtg                                                                                                       >  1:1177705/1‑151 (MQ=255)
                                              cAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGcc                                                                                                  >  1:655692/1‑151 (MQ=255)
                                               aaaTGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCt                                                                                                 >  2:205213/1‑151 (MQ=255)
                                                aaTGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCt                                                                                                 >  1:1145511/1‑150 (MQ=255)
                                                      tttCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTAttattt                                                                                          >  2:717083/1‑151 (MQ=255)
                                                        tcttctGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTAttatttat                                                                                        >  2:125863/1‑151 (MQ=255)
                                                           tctGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATgg                                                                                      >  1:199186/1‑150 (MQ=255)
                                                                aTTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGgtgt                                                                                <  2:552331/151‑1 (MQ=255)
                                                                  tttCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTc                                                                               <  1:77792/150‑1 (MQ=255)
                                                                   ttCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGgtgt                                                                                >  2:237213/1‑148 (MQ=255)
                                                                     cGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGtt                                                                            <  2:86926/150‑1 (MQ=255)
                                                                      gCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGttt                                                                           >  1:166915/1‑150 (MQ=255)
                                                                      gCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTg                                                                          <  1:795312/151‑1 (MQ=255)
                                                                        ttttGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGttttg                                                                         >  1:182622/1‑148 (MQ=255)
                                                                          ttGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCt                                                                       <  2:875740/150‑1 (MQ=255)
                                                                            gtgtCCCCCAAATCAGCGTAATGAATGTTTTGATTAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTg                                                                          <  1:832062/145‑1 (MQ=255)
                                                                                cccccACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGAt                                                                <  2:976142/151‑1 (MQ=255)
                                                                                cccccACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGGGTCGttttggtt                                                                      >  1:34651/1‑140 (MQ=255)
                                                                                  cccACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCa                                                              <  2:224324/151‑1 (MQ=255)
                                                                                      caTCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAgg                                                          <  2:941967/151‑1 (MQ=255)
                                                                                       aTCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCCGTTTGCTTGCAGATCTGAgg                                                          >  2:156751/1‑150 (MQ=255)
                                                                                         cAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTGGGTTGAAGATCAGAgg                                                          >  2:830343/1‑148 (MQ=255)
                                                                                          aGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGGTTGTAGAt                                                                >  1:562238/1‑141 (MQ=255)
                                                                                          aGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAgg                                                          >  2:311258/1‑147 (MQ=255)
                                                                                           gCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGgaagaa                                                     <  1:1038227/151‑1 (MQ=255)
                                                                                             gTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGAGTCGTTTGGGTTGAAGATCCGAGgaagaa                                                     >  1:508948/1‑149 (MQ=255)
                                                                                               aatgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGGAGATCAGAGGAAGaaa                                                    >  2:743147/1‑148 (MQ=255)
                                                                                               aatgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAg                                                  <  2:746981/150‑1 (MQ=255)
                                                                                               aatgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGAGTCGTCTTGCTTCAAGAAGAGGGGAAGaaaa                                                   >  1:381549/1‑149 (MQ=255)
                                                                                                atgaatgTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTCAAGATCAGGGGAAGaaa                                                    >  2:1220096/1‑147 (MQ=255)
                                                                                                      gTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGGCGTTTGGCTTGTAGATCAGAGGAAGGAAAAAGGCGag                                           >  2:266860/1‑150 (MQ=255)
                                                                                                        tttGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGagca                                         <  2:1288894/150‑1 (MQ=255)
                                                                                                           gATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCg                                      <  1:751414/150‑1 (MQ=255)
                                                                                                            aTCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTTAAAATCAAAGGAAGAAAAGATGCggg                                           >  2:501376/1‑142 (MQ=255)
                                                                                                              cAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGtttt                                  <  2:664471/151‑1 (MQ=255)
                                                                                                                aaCAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGGGTCGGTTTGCTTGGAGATCAGAGGGAGAAAAAAAGCGAGCAGCGGTTTgg                                >  2:1057792/1‑151 (MQ=255)
                                                                                                                  cAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGGTTTAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTtgg                               >  2:994574/1‑150 (MQ=255)
                                                                                                                  cAGAGGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGGTTTggg                               >  2:661912/1‑150 (MQ=255)
                                                                                                                   agagGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGGTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGGTTTggg                               >  2:855475/1‑149 (MQ=255)
                                                                                                                    gagGGCAAAAAAACAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGTACAGCTCTCGACGCCGCCCGCTGTGCTGCATATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGtttt                                  <  1:1067677/145‑1 (MQ=255)
                                                                                                                     agGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGGTTTGGGTTgg                           >  2:508630/1‑151 (MQ=255)
                                                                                                                     agGGCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAgg                           <  1:682107/151‑1 (MQ=255)
                                                                                                                      gggCAAAAAAATAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAgg                           <  1:1182656/150‑1 (MQ=255)
                                                                                                                          aaaaaaaTAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGtgtg                      <  1:742369/151‑1 (MQ=255)
                                                                                                                           aaaaaaTAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGc                     <  2:689791/151‑1 (MQ=255)
                                                                                                                               aaTAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGGTTTGGGTAGGGGGTGcc                    >  2:163985/1‑148 (MQ=255)
                                                                                                                                aTAGCCATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCtt                <  1:508616/151‑1 (MQ=255)
                                                                                                                                    ccATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTAtt            <  1:1039180/151‑1 (MQ=255)
                                                                                                                                    ccATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTAtt            <  1:1007492/151‑1 (MQ=255)
                                                                                                                                    ccATAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGCGTCGTTTGGCTTGGAGATCTGGGGGAGAAAAGAAGCGAGCAGCGCGTTGGGGGGGGGGTGCCGCTTTCtt            >  1:1195477/1‑151 (MQ=255)
                                                                                                                                      aTAAAAAAGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGAGGAGTTTTGCTTTTAGATAAGAGGAAGAAAAGAAGCGAGCAGCGTTTTTGGTAGGGGGTGCCGCTTTAt             >  2:143860/1‑148 (MQ=255)
                                                                                                                                        aaaaaaGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAAATGCGAGCAGCGTTTTGGGGGGGGGGTGCCGCTTTAttgtt         >  2:1123399/1‑150 (MQ=255)
                                                                                                                                        aaaaaaGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCGGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAATAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTAttgttg        <  2:789000/151‑1 (MQ=255)
                                                                                                                                         aaaaaGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGGGATGGCTGGCTTTGAGATCAGGGGGGGAAAAGAAGCGAGCAGCGGTTTGGGGAGGGggtgtc                    >  2:965666/1‑136 (MQ=255)
                                                                                                                                          aaaaGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTATTGTTGTc      <  1:17874/151‑1 (MQ=255)
                                                                                                                                           aaaGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTTGCTTTAAGATCAGAGGAAGAAAAAATGCGAGCAGCGGTTTGGGGTGGGGGTGGCGCTTTTttgtt         >  1:380916/1‑147 (MQ=255)
                                                                                                                                             aGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGc   <  1:1234355/151‑1 (MQ=255)
                                                                                                                                             aGTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGc   <  2:672390/151‑1 (MQ=255)
                                                                                                                                              gTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCg  <  1:494011/151‑1 (MQ=255)
                                                                                                                                              gTAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGGGTTGGTTGGCTTGAAGATCAGAGGGAGAAAACATGGGCGCAGCGTTTTGGGGTGGGGGTGGCGCCTTATTGTTGTTgggg  >  1:879331/1‑149 (MQ=255)
                                                                                                                                               tAAACATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGATTCGTTTGGCTTGTAGAGCAGAGGAAGAAAAAATGCGAGGAGCGATTTTGGGAGGGGGGGCCGCTtttttttttttgg    >  2:1175436/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGCACGAAATTGCTGCGCGCCCAGTAGTAATCTTTCAATTTTAGCAAATGGCTTTCTTCTGCATTTTCGCTTTTGTGTCCCCCACATCAGCGTAATGAATGTTTTGATCAAACAGAGGGCAAAAAAATAGCCATAAAAAAGTAAAAATGCGGGGCAACGCGCAGGACAGCTCTCGACGCCGCCCTCTGTGCTGCCTATTATTTATGGCGGTGTCGTTTGGCTTGAAGATCAGAGGAAGAAAAGATGCGAGCAGCGTTTTGGGTAGGGTGTGCCGCTTTATTGTTGTCGGCG  >  NC_000913/1440813‑1441105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑